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load  spectra  and  structural  properties.  Analytical  procedures  and  guidelines 
are  developed  and  described  for  quantifying  the  extent  of  damage  due  to  fatigue 
cracking  in  structural  details  as  a  function  of  service  time  and  design  variables  . 
^  fatigue  quality  model  is  described  £ind  discussed.  This  model  is 

based  on  the  tlme-to-crack-lnitiation  (TTCI)  and  the  equivalent  initial  flaw 
size  (EIFS)  concepts  and  a  deterministic  crack  growth  power  law.  An  experimenta;. 
test  program  and  the  resulting  fractographlc  data  for  cracks  in  clearance-fit 
fastener  holes  are  described  and  discussed.  The  durability  analysis  methodology 
is  demonstrated  for  a  full-scale  fighter  structure  (F-16  lower  wing  skin)  and  for 
a  complex  splice  subjected  to  a  B-1  bomber  spectrum.  Various  durability— related 
studies  are  described  and  discussed:  Initial  fatigue  model  calibrations,  statis¬ 
tical  crack  growth,  scale-up  and  hole  interaction  effects,  development  of 
crack  growth  transfer  functions  for  different  load  spectra  and  parameter  sensi¬ 
tivity  analysis. 

The  durability  analysis  methodology  has  been  demonstrated  for  quantifying  the 
extent  of  damage  due  to  relatively  small  cracks  in  clearance-fit  fastener  holes 
(e.g.,  *  O.IOV).  Further  research  is  needed  to  extend  the  methodology  to 
larger  crack sizes  and  to  verify  the  methodology  for  other  details,  such  as 
ffff®bs,  cutouts  and  lugs.  The -effects  of  fretting,  faying  surface  sealant, 
fastener  clamp-up,  environment,  interference -fit  fasteners,  etc.,  on  the  Initial 
fatigue  quality  needs  to  be  investigated.  A  fundamentally-sound  statistically- 
based  strategy  has  been  developed  for  quantifying  structural  durability  in 
meaningful  terms  for  judging  economic  life  and  optimizing  life-cycle  costs. 

Further  improvements  and  refinements  of  the  procedures  and  durability  design 
data  can  be  developed  within  the  general  framework  established. 
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FOREWORD 


This  report  was  prepared  by  General  Dynamics,  Fort  Worth 
Division,  with  the  support  of  George  Washington  University 
(Dr.  J.  N.  Yang)  and  Modern  Analysis  Inc.  (Dr.  M.  Shinozuka) . 
The  Air  Force  Wright  Aeronautical  Laboratories  (AFWAL/FIBEC) 
sponsored  this  research  under  the  "Durability  Methods  Develop¬ 
ment"  program  (Air  Force  Contract  F33615-77-C-3123) . 

James  L.  Rudd  was  the  Air  Force  Project  Engineer  and  Dr.  Jack 
W.  Lincoln  of  ASD/ENFS  was  a  technical  advisor  for  the  program. 
Dr.  B.  G.  W.  Yee  of  the  General  Dynamics  Materials  Research 
Laboratory  was  the  Program  Manager  and  Dr.  Sherrell  D.  Manning 
was  the  Principal  Investigator.  Dr.  J.  N.  Yang  and  Dr. 

M.  Shinozuka  were  associate  investigators. 

All  tests  were  performed  in  General  Dynamics '  Metallurgy 
Laboratory  by  R.  0.  Nay  under  the  direction  of  F.  C.  Nordquist. 
W.  T.  Kaarlela  was  responsible  for  the  f ractographic  data  ac¬ 
quisition.  Fractographic  readings  were  made  by  D.  E.  Gordon, 

W.  T.  Kaarlela,  A.  Meder,  R.  0.  Nay  and  S.  M.  Speaker. 

S.  M.  Speaker  coordinated  the  testing  and  fractographic  effort 
for  the  program.  J.  W.  Norris  developed  the  computer  software 
for  storing  and  analyzing  the  fractographic  data  acquired  and 
supported  the  initial  fatigue  quality  model  calibration/ 
evaluation  studies.  Dr.  V.  D.  Smith  supported  the  statistical 
analysis  effort.  Peggy  Thomas  and  Ernestine  Bruner  typed  the 
report  and  Ron  Jordan  prepared  many  of  the  illustrations. 

This  report  (Vol.  VII)  summarizes  the  work  performed 
under  Phase  II  of  the  "Durability  Methods  Development"  program. 
The  Phase  II  effort  is  further  documented  in  two  volumes  as 
follows : 

o  Volume  VIII  -  Test  and  Fractographic  Data 

o  Volume  IX  -  Documentation  of  Durability  Analysis  Computer 
Program. 

An  additional  volume  (VI)  was  also  prepared  to  document  the 
computer  programs  developed  under  Phase  I  for  Vol.  IV  - 
Initial  Quality  Representation. 

Work  for  this  report  was  performed  from  June  1979  to 
January  1984. 
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TERMINOLOGY 


Crack  Size  -  is  the  length  of  a  crack  in  a  structural 
detail  in  the  direction  of  crack  propagation. 

Deterministic  Crack  Growth  -  Crack  growth  parameters  are 
treated  as  deterministic  values  resulting  in  a  single 
value  prediction  for  crack  length. 

Durability  -  is  a  quantitative  measure  of  the  airframe's 
resistance  to  fatigue  cracking  under  specified  service 
conditions.  Structural  durability  is  normally  concerned 
with  relatively  small  subcritical  crack  sizes  which 
affect  functional  impairment,  structural  maintenance 
requirements  and  life-cycle-costs.  Such  cracks  may  not 
pose  an  immediate  safety  problem.  However,  if  the 
structural  details  containing  such  cracks  are  not 
repaired,  economical  repairs  cannot  be  made  when  these 
cracks  exceed  a  limiting  crack  size.  The  entire 
population  of  structural  details  in  various  components 
is  susceptible  to  fatigue  cracking  in  service. 
Therefore,  a  statistical  approach  is  essential  to 
quantitatively  assess  the  "durability"  of  a  part,  a 
component,  or  an  airframe. 

Durability  Analysis  -  is  concerned  with  quantifying  the 
extent  of  structural  damage  due  to  fatigue  cracking  for 
structural  details  (e.g..  fastener  hole,  fillet,  cutout, 
lug,  etc.)  as  a  function  of  service  time.  Results  are 
used  to  ensure  design  compliance  with  Air  Force's 
durability  design  requirements. 

Economic  Life  -  is  that  point  in  time  when  an  aircraft 
structure's  damage  state  due  to  fatigue,  accidental 
damage  and/or  environmental  deterioration  reaches  a 
point  where  operational  readiness  goals  cannot  be 
preserved  by  economically  acceptable  maintenance  action. 

Economic  Life  Criteria  -  are  guidelines  and  formats  for 
defining  quantitative  economic  life  requirements  for 
aircraft  structure  to  satisfy  U.S.  Air  Force  durability 
design  requirements.  The  economic  life  criterion 
provides  the  basis  for  analytically  and  experimentally 
ensuring  design  compliance  of  aircraft  structure  with 
durability  design  requirements.  Two  recommended  formats 
for  economic  life  criteria  are: 

o  probability  of  crack  exceedance 

o  cost  ratio:  repair  cost/replacement  cost 
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7. 


Economic  Repair  Limit  -  is  the  maximum  damage  size  that 
can  be  economically  repaired  (e.g.,  repair  0.03"  -  0.05" 
radial  crack  in  fastener  holes  by  reaming  hole  to  next 
size) . 

8.  Equivalent  Initial  Flaw  Size  (EIFS)  -  is  a  hypothetical 
crack  assumed  to  exist  in  the  structure  prior  to 
service.  It  characterizes  the  equivalent  effect  of 
actual  initial  flaws  in  a  structural  detail.  It  is 
determined  by  back-extrapolating  f ractographic  results. 
An  equivalent  initial  flaw  is  ass\amed  to  have  the  same 
flaw  shape  and  origin  as  the  observable  crack  size  at  a 
given  time.  The  EIFS  concept  is  a  convenient 
"mathematical  tool"  for  quantifying  the  IFQ  for 
structural  details  and  the  probability  ^  of  crack 
exceedance  or  extent  of  damage  as  a  function  of  time. 
An  EIFS  is  strictly  a  mathematical  quantity  rather  than 
an  actual  initial  flaw  size.  Within  this  context, 
EIFS's  can  be  positive  or  negative,  depending  on  the 
f ractographic  results  and  the  back  extrapolation  method 
used.  EIFS  values  depend  on  several  factors,  including: 
the  f ractographic  results  used  (and  the  test  variables 
reflected),  the  f ractographic  crack  size  range  used,  the 
form  of  the  crack  growth  equation  used  for  the  back- 
extrapolation,  the  goodness  of  the  curve  fit  to  the 
f ractographic  data,  the  manufacturing  quality  of  the 
structural  details,  fastener  hole  type,  fastener  type 
and  fit,  etc.  EIFSs  for  different  f ractographic  data 
are  not  comparable  unless  the  applicable  IFQ  model 
parameters  are  determined  consistently  (e.g.,  same 
f ractographic  crack  size  range  used,  same  "b  "•  value 
imposed  if  the  Q[a(t)3°  crack  growth  model  is  used, 
same  a  values  imposed  for  comparable  f ractographic  data 
sets,  etc.).  Fractographic  data  pooling  is  essential  to 
quantify  the  IFQ  for  different  fractographic  data  sets 
on  a  common  baseline. 

9.  EIFS  Master  Curve  -  is  a  curve  (e.g.,  equation, 
tabulation  of  a(t)  vs.  t  or  curve  without  prescribed 
functional  form)  used  to  determine  the  EIFS  value  at  t=0 
corresponding  to  a  given  TTCI  value  at  a  specified  crack 
size.  Such  a  curve  is  needed  to  determine  the  IFQ 
distribution  from  the  TTCI  distribution.  The  EIFS 
master  curve  depends  on  several  factors,  such  as  the 
fractographic  data  base,  the  fractographic  crack  size 

range  used,  the  functional  form  of  the  crack  growth 
equation  used  in  the  curve  fit,  etc.  (Ref.  EIFS). 
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10.  Extent  of  Damage  -  is  a  quantitative  measure  of 
structural  durability  at  a  given  service  time.  For 
example,  the  number  of  structural  details  (e.g., 
fastener  holes,  cutouts,  fillets,  etc.)  or  percentage  of 
details  exceeding  specified  crack  size  limits.  Crack 
length  is  the  fundamental  measure  for  structural  damage. 
The  predicted  extent  of  damage  is  compared  with  the 
specified  economic  life  criterion  for  ensuring  design 
compliance  with  U.S.  Air  Force  durability  requirements. 

11.  Generic  EIFS  Distribution  -  An  EIFS  distribution  is 
"generic"  if  it  depends  only  on  the  material  and 
manufacturing/fabrication  processes.  Theoretically,  the 
EIFS  distribution  should  be  independent  of  design 
variables,  such  as  load  spectrum,  stress  level,  percent 
load  transfer,  environment,  etc.  For  "durability 
analysis",  the  EIFS  distribution  for  fastener  holes 
(e.g.,  given  material,  drilling  procedures,  fastener 
types/fit,  etc.)  should  be  justified  for  different 
design  stress  levels  and  load  spectra. 

12.  Initial  Fatigue  Quality  (IFQ)  -  characterizes  the 

initial  manufactured  state  of  a  structural  detail  or 
details  with  respect  to  initial  flaws  in  a  part, 
component,  or  airframe  prior  to  service.  The  IFQ, 

represented  by  an  equivalent  initial  flaw  size  (EIFS) 
distribution,  must  be  defined  using  a  consistent 
f ractographic  data  base.  The  EIFS  distribution  depends 
on  the  f ractographic  crack  size  range  used  and  other 
factors  (Ref.  EIFS  and  EIFS  master  curve).  Whatever 
EIFS  distribution  is  used,  it  should  be  defined 
specifically  for  the  crack  size  range  of  interest  for 
the  structural  details  to  be  used  in  the  durability 
analysis.  A  single  EIFS  distribution  will  not 
necessarily  be  satisfactory  for  a  wide  range  of  crack 
sizes  (e.g.,  0.0005"  -  0.10").  Based  on  current 

understanding,  the  EIFS  distribution  should  be  defined 
for  a  fairly  small  range  of  crack  sizes  (e.g.,  0.020"  - 
0.050"  crack  size  for  fastener  holes).  Further  research 
is  required  to  evaluate  the  effects  and  sensitivity  of 
the  crack  size  range  on  the  EIFS  distribution  and  the 
accuracy  of  the  crack  exceedance  prediction. 

13.  Initial  Fatigue  Quality  Model  -  is  a  "mathematical  tool" 
for  quantifying  the  IFQ  distribution  for  applicable 
structural  details.  Using  the  IFQ  model  and 
f ractographic  results,  an  EIFS  distribution  can  be 

determined  which  is  compatible  with  the  TTCI 
distribution. 
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14.  Probability  of  Crack  Exceedance  (p(i.T))  -  refers  to  the 
probability  of  exceeding  a  specified  crack,  x,  ,  size  at 
a  given  service  time,  t.  It  can  be  determined  from  the 
statistical  distribution  of  crack  sizes  and  can  be  used 
to  quantify  the  extent  of  damage  due  to  fatigue  cracking 
in  fastener  holes,  cutouts,  fillets,  lugs,  etc. 

15.  Reference  Crack  Size  (a.)  -  This  is  the  specified  crack 

size  in  a  detail  used  to  reference  TTCI's.  The  IFQ 

distribution  is  based  on  a  selected  reference  crack 
size. 

16.  Service  Crack  Growth  Master  Curve  (SCGMC)  -  This  curve 
is  used  to  determine  the  EIFS,  yn(r),  corresponding  to 
an  exceedance  crack  size  x,  at  time  T.  The  probability 
of  crack  exceedance,  p(i,T),  can  be  determined  from  the 
EIFS  cumulative  distribution  for  a  given  yii(r).  The 
SCGMC  is  defined  for  the  applicable  design  variables 
(e.g.,  stress  level,  spectrum,  etc.)  and  it  can  be 
determined  using  either  test  data  or  an  analytical  crack 
growth  program.  All  SCGMC' s  must  be  consistent  with  the 
corresponding  EIFS  master  curve  and  the  f ractographic 
data  base.  The  SCGMC  must  be  consistent  with  the  basis 
for  the  IFQ  distribution. 

17.  Structural  Detail  -  is  any  element  in_  a  metallic 
structure  susceptible  to  fatigue  cracking  (e.g., 
fastener  hole,  fillet,  cutout,  lug,  etc.). 

18.  Time-To-Crack-Initiation  (TTCI)  -  is  the  time  or  service 
hours  required  to  initiate  a  specified  (observable) 
fatigue  crack  size,  Sq ,  in  a  structural  detail  (with  no 
initial  flaws  intentionally  introduced). 

J.9.  TTCI  Lower  Bound  Limit  (c)  -  e  is  a  cutoff  value  for 
TTCI's  reflected  in  the  IFQ  model.  It  varies  for_  a 
given  a,  and  it  depends  on  the  EIFS  upper  bound  limit, 
Xu,  and  the  EIFS  master  curve.  TTCI's  for  a  given  crack 
size,  ao,  should  ^  e.  This  Weibull  distribution 
parameter  provides  a  basis  for  quantifying  the  EIFS 
distribution  for  different  TTCI  crack  sizes  on  a  common 
baseline. 

20.  Upper  Bound  EIFS  (x,,  )  -  defines  the  largest  EIFS  in  the 
initial  fatigue  quality  distribution,  x^  ,  specified  by 
the  user,  should  be  consistent  with  e  (TTCI  lower  bound 
limit)  and  the  EIFS  master  curve. 
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SUMMARY 


Objectives  of  the  "Durability  Methods  Development" 
program  were:  (1)  develop  and  demonstrate  a  durability 
analysis  methodology  capable  of  satisfying  the  U.S.  Air 
Force's  durability  design  requirements  for  metallic 
airframes  and  (2)  prepare  a  durability  design  handbook. 

A  statistically-based  durability  analysis  methodology 
was  developed  under  Phase  I  of  this  program  [1-5].  This 
Volume  (VII)  documents  the  Phase  II  effort,  which  was 
concerned  with:  (1)  evaluating,  refining  and  demonstrating 
the  durability  analysis  method  developed  in  phase  I  for 
full-scale  aircraft  structures,  (2)  acquiring  f ractographic 
data  for  fatigue  cracks  in  clearance-fit  fastener  holes,  (3) 
developing  and  evaluating  durability  analysis  procedures  and 
guidelines  and  (4)  performing  related  durability  studies  to 
further  the  understanding.  Raw  test  and  f ractographic  data 
are  reported  in  Volume  VIIlK]  Computer  programs  for 
implementing  the  durability  analysis  method  are  listed  in 
Volume  IX  C7J 

Aircraft  structure  must  be  designed  to  be  durable  to 
minimize  structural  maintenance  requirements,  life-cycle- 
costs  and  functional  impairment  problems  affecting 
operational  readiness.  Durability  is  normally  concerned 
with  relatively  small  subcritical  fatigue  cracks  in 
structural  details  (e.g.,  0.03"-  0.05"  radial  crack  in  a 
fastener  hole).  It  is  a  measure  of  the  structure's 
resistance  to  fatigue  cracking. 

The  population  of  structural  details,  such  as  fastener 
holes,  fillets,  cutouts,  lugs,  etc., in  aircraft  structures 
is  susceptible  to  fatigue  cracking  in  service.  Likewise, 
structural  maintenance  requirements  and  costs  depend  on  the 
number  of  details  requiring  repair.  To  assess  the 
durability  of  the  structure  or  the  extent  of  damage  as  a 
function  of  time,  .the  entire  population  of  structural 
details  must  be  accounted  for.  Therefore,  a  statistical 
approach  is  essential  to  quantify  the  overall  extent  of 
damage  as  a  function  of  time. 

According  to  U.S.  Air  Force  durability  design 
requirements  for  metallic  airframes,  the  economic  life  of 
the  structure  must  exceed  one  design  service  life.  Economic 
life  is  currently  defined  in  qualitative  terms:  "...  the 
occurrence  of  widespread  damage  which  is  uneconomical  to 
repair  and,  if  not  repaired,  could  cause  functional  problems 
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affecting  operational  readiness"  C  8-10],  Acceptable  limits 
for  "widespread  damage"  and  "uneconomical  repair"  must  be 
defined  for  each  aircraft  design  and  such  limits  must  be 
approved  by  the  Air  Force.  A  quantitative  definition  of 
economic  life  has  not  evolved  from  this  program.  However, 
guidelines  and  two  basic  formats  for  quantitative  economic 
life  criteria  have  been  developed:  (1)  probability  of  crack 
exceedance  and  (2)  cost  ratio:  repair  cost/replacement  cost. 
The  analytical  tools  developed  under  this  program  can  be 
used  to  quantify  the  extent  of  damage  and  show  compliance 
with  the  specified  economic  life  criterion. 

The  durability  analysis  method  developed  is  conceptually 
described  in  Fig.  1.  It  accounts  for  initial  fatigue 
quality,  crack  growth  acciamulation  in  a  population  of 
structural  details,  load  spectra  and  structural  properties. 

A  "mathematical  tool"  has  been  developed  for  quantifying  the 
initial  fatigue  quality  of  structural  details  using 
f ractographic  results.  This  model  is  based  on  the  time-to- 
crack-initiat ion  (TTCI)  and  equivalent  initial  flaw  size 
(EIFS)  concepts  and  a  deterministic  crack  growth  power  law. 
Once  the  initial  fatigue  quality  or  EIFS  distribution  has 
been  defined  for  a  given  detail  type,  the  EIFS  population 
for  groups  of  details  with  a  common  stress  history  can  be 
grown  to  a  given  service  time  (for  specified  design 
conditions)  using  a  suitable  analytical  crack  growth 

program.  The  extent  of  damage  is  measured  by  the  nxmiber  of 
details  in  the  population  expected  to  exceed  a  specified 
crack  size  at  a  given  service  time  or  equivalently  by  the 
probability  of  crack  exceedance.  Simple  expressions,  based 
on  the  binomial  distribution,  are  used  to  quantify  the 
overall  extent  of  damage  for  different  details,  oarts, 
components,  or  the  entire  airframe.  Statistical  principles 
can  be  used  to  estimate  upper  and  lower  bound  predictions  for 
the  extent  of  damage.  These  results  provide  a  quantitative 
description  of  the  structure's  resistance  to  cracking  and  a 
basis  for  judging  economic  life. 

The  durability  analysis  method  has  been  evaluated  and 
demonstrated  for  relatively  small  fatigue  cracks  in 
clearance-fit  fastener  holes  (e.g.,  ^  0.10").  However,  the 
basic  strategy  and  concepts  theoretically  apply  to  fatigue 
cracks  in  any  structural  detail  (e.g.,  cutouts,  fillets, 
lugs,  etc.)  as  long  as  the  dominant  crack  in  each  detail  is 
relatively  small  and  such  cracks  are  statistically 
independent  (i.e.,  the  growth  of  the  largest  crack  in  one 
detail  does  not  significantly  affect  the  growth  of  the 
largest  crack  in  neighboring  details  and  vice  versa).  The 
durability  analysis  method  has.  been  successfully 
demonstrated  for  a  full-scale  fighter  aircraft  structure 
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Essential  Elements  of  the  Durability  Analysis  Methodology 
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Fig.  1  Essential  Elements  of  the  Durability  Analysis  Methodology  (Continued) 


(i.e.,  F-16  lover  wing  skin)  and  for  a  complex  splice 

subjected  to  a  B-1  bomber  spectrum. 

The  analytical  tools  and  concepts  developed  under  this 
program  can  be  used  to:  (1)  evaluate  durability  design  trade¬ 
off  options  in  terms  of  structural  design  variables,  (2) 
support  the  full-scale  durability  test  plan/evaluations,  (3) 
evaluate  structural  maintenance  requirements  before  or  after 
the  aircraft  is  committed  to  production  or  service,  and  (4) 
provide  aircraft  user  options  affecting  life-cycle-costs, 
structural  maintenance  requirements  and  operational 
readiness. 

A  large  f ractographic  data  base  was  developed  under  this 
program  for  fatigue  cracks  in  clearance-fit  fastener  holes 
(straight-bore  and  countersunk)  to:  (1)  evaluate  and 
demonstrate  the  durability  analysis  methodology  developed 
(2)  provide  durability  design  data  for  the  handbook. 
Fractographic  results  for  over  800  fatigue  cracks  were 
generated.  The  following  test  parameters  were  considered: 

o  Material  {7475-T7351  Al,  2024-T851  Al, 

D6ac  Steel) 

o  Load  Spectra  (F-16  400-hour,  F-16  500-hour, 

B-1  bomber) 

o  %  Shear  Load  Transfer  Specimens  (0%,  15%,  30%, 

40%) 

o  Clearance-fit  Fasteners  (3/16"  Dia.  and  1/4"  Dia.) 

-  Countersunk  Rivet  (MS  90353) 

-  Countersunk  Bolt  (NAS  1580) 

-  Protruding  Head  Bolt  (NAS  6204) 

o  Multiple  Stress  Levels 

o  Environment  (Lab  Air) 

A  tear-down  inspection  of  the  F-16  lower  wing  skins  from  the 
durability  test  article  was  performed.  Over  1600  fastener 
holes  per  wing  skin  were  inspected  for  cracks.  Results  are 
sxmunarized  in  this  report  and  details  are  given  in  Volume 
VIII . 


Several  studies  were  performed  to  evaluate  the 
durability  analysis  methodology:  initial  fatigue  quality 
model  calibrations-,  statistical  crack  growth,  scale-up  and 
hole  interaction  effects,  crack  growth  transfer  functions 
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for  different  load  spectra  and  parameter  sensitivity 
analyses. 

Further  research  is  required  to:  (1)  extend  the 
durability  analysis  method  to  larger  crack  sizes  (e.g.,  > 
0.10"),  (2)  demonstrate  the  method  developed  for  cutouts, 
fillets,  lugs,  etc.,  and  (3)  increase  confidence  in  the 
extent  of  damage  predictions  for  different  design  conditions 
(e.g.,  materials,  load  spectra,  stress  levels,  etc.).  The 
effects  and  significance  of  factors  such  as  fretting, 
faying  surface  sealant,  environment,  interference  “  fit 
fasteners,  fastener  clamping  force,  fastener  diameter, 
coupon  size,  etc. , on  initial  fatigue  quality  need  to  be 
investigated.  A  better  understanding  of  the  initial  fatigue 
quality  model  parameter  values,  range  of  values  and 
parameter  sensitivities  is  needed  for  different  design 
conditions  and  structural  details.' 
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SECTION  I 


INTRODUCTION 


1.1  OBJECTIVE  AND  SCOPE 

This  report  (Vol.  VII)  documents  and  summarizes  the 
Phase  II  effort.  The  objectives  and  scope  of  this  report 
are; 

1.  Evaluate,  refine  and  demonstrate  the  durability 
analysis  methodology  developed  under  Phase  I  [1-5]  for 
full-scale  fighter  structure  (F-16  lower  wing  skin)  and  for 
a  complex  splice  subjected  to  a  B-1  bomber  spectrum. 

2.  Summarize  the  essential  elements  of  the  refined 
durability  analysis  method  and  describe/illustrate 
implementation  procedures, 

3.  Discuss  important  aspects  of  durability  analysis. 

4.  Describe  the  Phase  II  experimental  program, 
f ractographic  data,  and  evaluate  results.  Raw 
test/f ractography  data  are  reported  in  Volume  VIII. 

5.  Document  various  durability-related  studies  that 
were  performed  to  develop  a  better  understanding  of  the 
durability  analysis  method/procedures.  Key  issues, 
guidelines,  limitations,  overall  conclusions  and  future 
research  needs  are  discussed. 

The  Phase  II  effort  is  dociimented  in  nine  quarterly  pro¬ 
gram  progress  reports  [11-19]  and  two  other  volumes  as  follows 

•  Volume  VIII  -  Test  and  Fractography  Data  [6  ] 

•  Volume  IX  -  Docxamentation  of  Durability  Analysis 

Computer  Program  [  7j 

Volume  VI,  "Doc\amentation  of  Computer  Programs  For  Initial 
Quality  Representation  (Vol.  IV)  [20],  was  also  prepared 
under  Phase  II  to  document  the  computer  programs  developed 
under  Phase  I  for  Volume  IV. 

A  Durability  Design  Handbook  [21]  was  prepared  under 
Phase  III.  This  handbook  has  a  similar  format  to  the  USAF's 
Damage  Tolerance  Handbook  [22].  It  siammarizes  the  Air 
Force's  durability  design  requirements  [s-lO  ]  and  includes 


1 


durability  analysis  methodology,  design  data  and  guidelines 
for  satisfying  these  requirements. 

The  final  understandings  for  this  program,  documented  in 
the  Durability  Design  Handbook  [21),  are  based  bn  the  Phase  I 
and  Phase  II  effort.  Many  studies  were  performed  during 
Phase  II  to  further  the  understanding  of  aircraft  structural 
durability,  durability  analysis  and  related  issues  (e.g., 
initial  fatigue  quality,  generic  EIFS,  IFQ  model  calibration 
procedures,  etc.).  Phase  II  was  definitely  a  "learning 
process."  Therefore,  some  of  the  studies  and  understandings 
from  the  Phase  II  effort  are  documented  in  this  Volxome  for 
information  purposes  only.  In  any  case,  the  Durability 
Design  Handbook  [21]  is  the  governing  document  for 
durability  analysis  applications. 


1.2  Organization 

This  report  is  organized  into  seven  sections  and  three 
appendices.  In  general,  the  Phase  II  effort  is  summarized 
and  discussed  in  the  main  text;  further  details  are 
dociamented  in  applicable  appendices.  Volumes  VIII  and  IX, 
quarterly  reports  and  technical  papers. 

Durability  background  information  is  presented  in 
Section  II;  the  durability  problem,  essential  Air  Force 
durability  design  requirements,  economic  life,  extent  of 
damage,  and  durability  analysis  objectives. 

Essential  elements  of  the  durability  analysis  method 
developed  under  Phase  I  and  applicable  refinements  made 
under  Phase  II  are  summarized  and  discussed  in  Section  III. 

Durability  analysis  details  are  discussed  in  Section  IV, 
including  related  durability  studies.  Procedures  and 
guidelines  are  also  described  for  implementing  the 
durability  analysis  method  described  in  Section  III. 
Initial  fatigue  quality  model  studies  and  the  demonstration 
of  the  durability  analysis  method  for:  (1)  an  f-16  lower  wing 
skin  and  (2)  complex  splice  specimens  (B-l  bomber  spectrum) 
are  discussed. 

The  durability  test  program  is  described  in  Section  V; 
test/f ractographic  results  are  evaluated  and  discussed  in 
Section  VI.  Raw  test/f ractographic  data  are  documented  in 
Volume  VIII  [6  ]. 

The  overall  program  conclusions  and  recommendations  for 
further  research  are  discussed  in  Section  VII.  This  is 

followed  by  three  appendix  sections  which  support  the 
detailed  documentation  of  various  durability  aspects. 
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SECTION  II 


BACKGROUND 


2.1  INTRODUCTION 

The  purpose  of  this  section  is  to:  (1)  describe  the 
aircraft  durability  problem  and  implications,  (2)  describe 
the  essential  elements  of  the  Air  Force  durability  design 
requirements,  (3)  discuss  durability^  analysis  objectives, 
and  (4)  discuss  the  meaning  of  economic  life. 


2.2  AIRCRAFT  STRUCTURAL  DURABILITY  PROBLEM 

Aircraft  structures  have  thousands  of  structural  details 
susceptible  to  fatigue  cracking:  fastener  holes,  fillets, 
cutouts,  etc.  For  example,  the  wing  box  assembly  shown  in 
Fig.  2  has  over  3000  fastener  holes  in  the  wing  skins  alone. 
Aircraft  structural  damage  surveys  conducted  under  this 
program  showed  that  fatigue  cracks  are  likely  to  occur  in 
fastener  holes  under  aircraft  service  conditions  [3]. 
Similar  conclusions  have  also  been  reported  elsewhere  [23- 
27]. 


Durability  is  a  measure  of  the  structure's  resistance  to 
fatigue  cracking.  The  entire  population  of  structural 
details  in  various  components  is  susceptible  to  cracking  in 
service.  Therefore,  to  assess  the  durability  of  the 
structure  or  extent  of  damage  (i.e.,  number  or  percentage  of 
structural  details  in  a  part,  a  structure,  a  component  or 
airframe  exceeding  specified  crack  size  limits  that  cannot 
be  economically  repaired)  as  a  function  of  time,  the  entire 
population  of  structural  details  must  be  accounted  for. 
Thus,  a  statistical  approach  is  essential  to  quantify  the 
extent  of  damage  as  a  function  of  time. 

Structural  durability  is  generally  concerned  with 
relatively  small  subcritical  crack  sizes  which  affect 
functional  impairment,  structural  maintenance  requirements 
and  life-cycle-costs.  Such  cracks  may  not  pose  an  immediate 
safety  problem.  However,  if  the  structural  details 
containing  such  cracks  are  not  repaired,  economical _  repairs 
cannot  be  made  when  these  cracks  exceed  a  limiting  crack 
size.  For  example,  a  0.030-0.050"  radial  crack  in  a 
fastener  hole  can  be  cleaned  up  by  reaming  the  hole  to  the 
next  fastener  size.  The  economical  repair  limit  is  the 
maximiom  crack  size  in  a  detail  that  can  be  cleaned-up 
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without  further  repair  or  part  replacement.  If  structural 
details  are  not  repaired  or  parts  replaced  at  an  opportune 
time,  expensive  repairs  or  part  replacement  may  be  required. 
Also,  unrepaired  cracks  may  reach  sizes  which  could  affect 
structural  safety  during  the  design  life  of  the  aircraft. 

Aircraft  structural  safety  is  ensured  by  satisfying 
damage  tolerance  requirements.  Damage  tolerance _  is 
typically  concerned  with  the  rogue  flaw  size  in  a  single 
detail.  For  example,  in  Fig.  2  the  damage  tolerance 
analysis  of  the  wing  box  is  limited  by  a  few  critical 
structural  details.  On  the  other  hand,  the  durability  of 
the  wing  box  is  concerned  with  the  entire  population  of 
structural  details  and  the  size  of  the  largest  subcritical 
crack  in  each  detail. 

The  conventional  fatigue  analysis  (CFA)  approach  (e.g., 
Palmgren-Miner  rule.  Ref.  28,29)  and  the  deterministic  crack 
growth  approach  (DCGA)  [22]  don't  have  the  proper  format  for 
quantifying  the  "extent  of  damage"  as  a  function  of  service 
time.  The  CFA,  in  its  commonly  used  form,  does  not  quantify 
crack  sizes  for  a  population  of  details  -  an  essential 
requirement  for  any  durability  analysis  method.  A  DCGA  can 
be  used  to  predict  the  growth  of  a  single  crack  in  a  detail 
as  a  function  of  time.  Using  the  DCGA,  details  can  be 
grouped  and  the  "worst-case"  detail  in  the  group  can  be  used 
to  analytically  assure  the  largest  crack  size  in  the  group 
of  details  will  be  =  a  specified  size.  However,  the  DCGA 
does  not  yield  a  quantitative  description  of  the  probable 
crack  sizes  or  ranges  of  crack  sizes  for  the  population  of 
details.  CFA  and  the  DCGA  have  been  evaluated  for  potential 
durability  analysis  applications  [2,30]. 


2.3  AIR  FORCE  DURABILITY  DESIGN  REQUIREMENTS 

Air  Force  durability  design  requirements  for  metallic 
airframes  are  presented  in  three  documents  [  8-10.], 
According  to  these  requirements,  the  airframe  must  be 
designed  to  have  an  economic  life  greater  than  the  design 
service  life.  Furthermore,  the  economic  life  must  be 
verified  by  analysis  and  test.  The  analytical  verification 
of  economic  life  must  account  for  the  effects  of  initial 
quality  variations,  material  property  variations  and  the 
design  loads/environments. 

Durability  analysis  tools  have  been  developed  under  this 
program  for  analytically  quantifying  the  extent  of  damage, 
due  to  fatigue  cracking,  as  a  function  of  time.  These  tools 
can  be  used  to  judge  the  economic  life  of  aircraft  structure 
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for  groups  of  details,  parts,  components,  or  an  entire 
airframe. 


2.4  ECONOMIC  LIFE 

The  economic  life  of  an  aircraft  structure  is  currently 
defined  in  qualitative  terms:  "...the  occurrence  of 
widespread  damage  which  is  uneconomical  to  repair  and,  if 
not  repaired,  could  cause  functional  problems  affecting 
operational  readiness"  [  8-10,].  Acceptable  limits  for 
"widespread  damage"  and  "uneconomical  repair"  must  be 
defined  for  each  aircraft  design  and  such  limits  must  be 
approved  by  the  Air  Force. 

A  quantitative  definition  of  economic  life  has  not 
evolved  from  this  program.  However,  guidelines  have  been 
developed  for  specifying  economic  life  criterion.  Various 
aspects  of  economic  life  and  recommended  formats  for 
economic  life  criterion  are  discussed  in  the  Durability 
Design  Handbook  [21].  In  any  case,  the  economic  life  of  an 
airframe  depends  on  the  user's  acceptable  limits  for  extent 

of  structural  damage,  structural  maintenance  requirements  and 
costs  and  operational  readiness. 


2.5  EXTENT  OF  DAMAGE 

The  extent  of  damage  is  a  quantitative  measure  of  the 
number  of  structural  details  exceeding  specified  crack  size 
limits  as  a  function  of  service  time.  Structural 
maintenance  requirements  and  costs  depend  on  the  number  of 
structural  details  requiring  repair.  The  "durability"  of 
the  structure  depends  on  the  extent  of  damage  .  for  the 
population  of  structural  details  in  a  part,  component,  or 
airframe. 

The  extent  of  damage  can  be  predicted  using  the 
analytical  tools  described  in  this  report  and  the  handbook 
[21].  Extent  of  damage  predictions  provide  the  basis  for 
analytically  ensuring  design  compliance  with  the  governing 
economic  life  criterion. 


2.6  DURABILITY  ANALYSIS  OBJECTIVES 

The  basic  objective  of  the  durability  analysis  tools 
developed  under  this  program  is  to  analytically  assure 
design  compliance  with  Air  Force  durability  design 
requirements  [  8-10  ]  during  the  design  stage  before  the 
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aircraft  structure  is  committed  to  durability  testing, 
production  or  service.  Other  objectives  are: 

1.  Evaluate  durability  design  tradeoffs  for  different 
materials,  stress  levels,  design  spectra,  manufacturing 
variables,  etc.,  which  affect  Air  Force  life-cycle-costs  and 
operational  readiness. 

2.  Support  the  durability  verification  plan  and 
evaluation  of  results  for  full-scale  aircraft  structures. 

3.  Develop  structural  maintenance  policy  before  the 
aircraft  is  committed  to  service. 
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SECTION  III 


DURABILITY  ANALYSIS  METHOD 


3.1  INTRODUCTION 

Essential  elements  and  equations  of  the  durability 
analysis  method,  including  Phase  II  refinements,  are 
presented  in  this  section.  Implementation  procedures  are 
described  and  discussed  in  Section  IV  and  in  the  handbook 
[21].  Further  details  are  given  elsewhere  [1,2,5,  11-19,30- 
36].  Recommended  durability  analysis  procedures  and 
guidelines  are  presented  in  the  Durability  Design  Handbook 
[21 ].  Reference  21  is  the  governing  document  for  durability 
analysis  applications. 


3.2  GENERAL  DESCRIPTION 

The  basic  objective  of  the  durability  analysis 

methodology  is  to  quantify  the  distribution  of  crack  sizes 
as  a  function  of  service  time  for  all  structural  details 
(e.g.,  fastener  holes,  cutouts,  fillets,  lugs,  etc.) 
experiencing  a  common  load  history,  environment,  etc.  The 
extent  of  damage  is  measured  by  the  number  of  details 
expected  to  have  a  crack  size  greater  than  any  specified 
size  at  a  given  service  time  or  equivalently  by  the 
probability  of  crack  exceedance.  Using  the  binomial 
distribution,  the  extent  of  damage  for  different  details, 
parts,  components,  or  an  entire  airframe  can  be  predicted. 
Upper  and  lower  bound  predictions  for  the  "extent  of  damage" 
can  be  estimated  using  statistical  analysis  principles.  The 
extent  of  damage  prediction  provides  a  basis  for  analytically 
assuring  design  compliance  with  the  selected  criterion  for 

economic  life.  ^ 

The  durability  analysis  includes  two  essential  steps: 
(1)  quantify  tl\e  initial  fatigue  quality  (IFQ)  distribution 
for  the  structural  details  to  be  included  in  the  extent  of 
damage  assessment  (Fig.  3),  and  (2)  predict  the  probability 
of  crack  exceedance  using  the  IFQ  distribution  and  the 
applicable  design  conditions  (e.g.,  load  spectrum,  stress 
level,  percent  load  transfer,  etc.)  (Fig.  4). 

3.3  ASSUMPTIONS 

The  durability  analysis  method  is  based  on  the  following 
explicit  or  implied  assumptions: 
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Fig.  3  Initial  Fatigue  Quality  Model 
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Note:  EIFSs  must  be  grown  forward  the 

same  way  the  fractographic  results 
were  grown  backward  to  determine 
the  EIFSs. 


Fig.  4  Growth  of  EIFS  Distribution  as  Function  of  Time 
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1.  Crack  length  is  the  fundamental  measure  of 
durability  damage. 

2.  Each  detail  in  the  population  of  details  in  an 
aircraft  structure  has  a  single  dominant  crack  size;  such 
cracks  are  considered  to  be  random  variables. 

3.  The  largest  crack  in  each  detail  is  relatively  small 
(e.g.,  ^0.05"  corner  crack  in  a  fastener  hole)  and  such 
cracks  are  statistically  independent.  Hence,  the  growth  of 
the  largest  crack  in  one  detail  does  not  significantly 
affect  the  growth  of  the  largest  crack  in  neighboring 
details  and  vice  versa. 

4.  The  initial  fatigue  quality  (IFQ)  of  the  structural 
details  to  be  included  in  the  durability  analysis  can  be 
represented  by  an  equivalent  initial  flaw  size  (EIFS) 
distribution. 

5.  The  population  of  dominant  crack  sizes  for  a  group 
of  structural  details  with  a  common  stress  history  can  be 
grown  from  one  time  to  another  using  a  single  deterministic 
crack  growth  curve. 

6.  Once  the  IFQ  distribution  has  been  defined,  a 
suitable  analytical  crack  growth  program  can  be  used  to 
predict  the  extent  of  damage  as  a  function  of  time  using 
applicable  design  variables  .  [e . g . ,  37-40]. 


3.4  INITIAL  FATIGUE  QUALITY 

Initial  fatigue  quality  (IFQ)  defines  the  initial 
manufactured  state  of  a  structural  detail  or  details  with 
respect  to  initial  flaws  in  a  part,  component  or  airframe 
prior  to  service.  The  objective  of  the  durability  analysis 
is  to  quantify  the  distribution  of  crack  sizes  for 
structural  details  (e.g.,  fastener  holes)  as  a  function  of 
time.  An  IFQ  distribution  provides  a  quantitative 
description  of  the  initial  flaws  in  the  structure  —  an 
essential  starting  point  for  all  durability  analyses. 

The  initial  fatigue  quality  for  a  group  of  replicate 
details  is  represented  by  an  equivalent  initial  flaw  size 
(EIFS)  distribution.  An  equivalent  initial  flaw  is  a 
hypothetical  crack  assumed  to  exist  in  a  detail  prior  to 
service.  An  arbitrary  crack  size,  ag,  is  selected  which  can 
be  readily  detected  or  which  can  be  reliably  observed 
f ractographically  following  testing.  The  time  required  for 
an  initial  defect,  of  whatever  type,  to  become  a  fatigue 
crack  of  size  ag  is  termed  the  time-to-crack-init iation 
(TTCI).  Test  results  of  TTCI  and  crack  growth  rates  using 
coupon  specimens  are  employed  to  define  the  EIFS 
distribution.  A  conceptual  description  of  the  initial 
fatigue  quality  model  is  shown  in  Fig.  3.  Once  the  EIFS 
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distribution  (or  population)  has  been  defined  for  a  group  of 
details,  a  deterministic  crack  growth  analysis  is  used  to 
grow  the  entire  EIFS  population  to  any  service  time  t;  thus 
determining  the  time-varying  crack  size  distribution  in 
service  (Fig.  4). 


3.5  INITIAL  FATIGUE  QUALITY  MODEL 

An  IFQ  model  was  developed  in  Phase  I  for  quantifying 
the  EIFS  distribution  for  structural  details,  such  as 
fastener  holes  [5,32].  This  model  was  evaluated  and  refined 
during  Phase  II.  IFQ  model  refinements  and  EIFS 

distributions  for  two  model  variations  will  now  be  described 

and  discussed.  The  two  EIFS  distributions  are  both  based  on 
the  three-parameter  Weibull  distribution  and  variations  of 
the  same  crack  growth  rate  equation.  The  IFQ  model 

conceptually  described  in  Fig.  3,  is  a  "mathematical  tool" 
for  quantifying  the  EIFS  distribution  for  structural 

details. 

The  time-to-crack-initiation  (TTCI)  distribution  for 
coupon  specimens  is  represented  by  the  three-parameter 
Weibull  distribution,  FT(t). 


F^Ct)  =  P[T<t]  =  1  -  exp 


where  T  =  TTCI,  a  =  shape  parameter,  /3  =  scale  parameter 
and  e  =  lower  bound  TTCI.  The  three  Weibull  parameters  are 
determined  from  fractography  for  coupon  specimens  or  other 
suitable  test  results. 

The  crack  growth  rate  • in  the  small  crack  region  is 
assumed  to  be  of  the  following  form 


(2) 


where  Q  and  b  are  parameters  depending  on  loading  spectra 
structural  and  material  properties,  etc.;  a(t)  is  the  crack 
size  at  t.  Other  crack  growth  laws  could  also  be  used.  Eq. 
2  IS  used  because  of  its  simplicity  and  general  abilitv  to 
fit  crack  growth  data. 


objective  of  the  IFQ  model  is  to  define  an  EIFS 
distribution  which  is  statistically  compatible  with  the  TTCI 
distribution.  A  crack  growth  law,  such  as  Eq.  2,  is  used  to 
obtain  the  EIFS  distribution  through  a  transformation  of  the 
^distribution.  Two  different  EIFS  distributions  are 
obtained  gsing  Eq.  2  for  b  *  1  (case  I)  and  b  =  1  (case 
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II).  The  resulting  EIFS  equations  for  case  I  and  case  II 
are  described  below. 

Case  I  (b  #  1) 

Integrating  Eq.  2  from  t  =  0  to  t  =  T,  the  relation 
between  the  initial  crack  si2e,a(0),  and  the  reference  crack 

size  aQ  =  a(T)  for  TTCI  (flight  hours)  is  obtained  as 
follows : 


EIFS  =  a(0)  =  (a”^+  cQT)"^^^^ 


(3) 


where  c  =  b  -  1. 

The  distribution  of  a{0)  can  be  obtained  from  that  of  T 
given  by  Eq.  1  through  a  transformation  of  Eq.  3  as  follows, 

1 


-[ 


-c  -c  ^ 
^  -  ^0  -  cQe 


cQ3 


a 

] 


;  0<x-x 


u 


(4) 


=  1 


where 


;  x>x 


u 


3  =  3j^(£) 


1/a 


(5) 


u  ^^0^  ‘^Qe) 


0 


-1/c 


(6) 


I n  Eq .  5 , 

/3£  =  Weibull  scale  parameter  for  TTCI 

based  on  the  TTCI's  for  a  given  f ractographic  data 
set  in  which  only  the  fractography  for  the  largest 
fatigue  crack  in  any  "one"  of  fastener  holes 
per  test  specimen  is  used. 

■I  =  Number  of  equally-stressed  fastener 

holes  per  test  specimen  in  which  only  the  largest 
fatigue  crack  in  any  "one"  hole  per  specimen  is 
included  in  the  f ractographic  data  set. 

and, 

a  =  Weibull  shape  parameters  for  a 
given  TTCI  data  set. 


14 


Details^  for  Eq.  5  are  given  in  Section  4.4  and  Ref. 14.  in 
Eq.  4,  if  X  <  0  (i.e.,  negative  EIFS),  refer  to  Section  3.6 
(Step  4)  and  subsection  4.2.3. 

In  Eq.  6,  is  the  upper  bound  EIFS  in  the  initial 
fatigue  quality  ( IF©' distribution.  The  lower  bound  TTCI ,  e, 
for  a  given  reference  crack  size,  an,  is  defined  from  Eg.  6 
as  follows: 


Eq.  4  for  1^®  simplified  by  substituting  the 

expression  for  CQe  from  Eq.  6  into  Eq.  4  as  follows: 


a  > 


In  Eq.  8,  p  and  Xu  are  defined  by  Eqs.  5  and  6, 
respectively. 


Case  II  (b  =  1) 

Integrating _  Eq.  2  from  t  =  0  to  t  =  T  and  considering  b 
=  1,  the  following  relation  between  the  initial  crack  size, 
a(0),  and  the  reference  crack  size  an  ■»  a(T)  for  TTCI 
(flight  hours)  is  obtained: 


EIFS  =  a(0)  =  exp  (-Qx) 


(9) 


The  EIFS  cumulative  distribution,  Fafo')^*),  is  obtained 
by  use  of  Eqs.  1  and  9  as  follows: 


1  -  F 


''a(O) 


=  1 


I 


rin^^ 

L  Q6 


Sq/x)  -  Qe 


^OL\ 


;  0<  x<  X 


u 


x>x 


u 


I 

) 


(10) 


defined  by  Eq,  5  and  Xu  is  the  selected  upper 

(11) 


X  =  a-  exp  (-Qe) 
u  0 
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where  jS  is 
bound  EIFS: 


The  lover_  bound  TTCI  value,  e,  for  a  given  reference  rr«rw 
Size,  aQ,  IS  defined  from  Eq.  11  as  follows: 


;  a  >  X 
0  U 


(12) 


aSd  simplified  by  combining  Eqs.  10 


In  (x^/x)  ]  \ 

.  06  J  )’ 


(13) 


P  and  in  Eq.  13  are  defined  by  Eq.  5  and  6, 


respectively. 


3 . 6  GENERAL  PROCEDURE 

The  general  procedure  for  implement inq  the  durahni^-tr 

euei^e^e  ""Iff  discSIsed"  beloS'^anl 

exiewnere  _  11,5,^1,34-36  3  .  For  final  details  =na 
recomnendations,  refer  to  the  Durability  Design  HandboSk 


regions  ^T.trT  .1^11 

n?g??^^L^en^^?To1e?! 


2.  Use  the  IFQ 
f ractographic  results 
distribiton,  Fa(o)(x). 
Eqs.  3-8.  For  Case  II 


model  shown  in  Fig.  3  and  suitable 
rr  ^®^^ne  the  EIFS  cumulative 
If  Case  I  (i.e,,  b  9^  1)  is  used,  use 
U.e.,  b  =  1),  use  Eqs.  9-13. 


(SCGMC)  service  crack  growth  master  curve 
(SCGMC)  to  determine  ,  for  each  stress  region  tha 
corresponding  EIFS  value,  ^^(r),  which  grows  to®c?ack  si« 
^1  e • ^  service  time  t  as  illustrated  in  Fig  4  rf 
applicable  f ractographic  data  are  available  for 'different 

Sser^thrsCGMc''L^''S''^d^f®^^-®  pooling  procedures  are 
used,  the  SCGMC  can  be  determined  from  Eq.  14  as  follows. 


Two  different  SCGMC s  can  be 
considering  b^  9^  1  and  b^  =  1. 


obtained 


from  Eq.  14 


da(t) 

dt 


Qi 


(14) 
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For  example,  Eq.  14,  for  h±*  1,  can  be  integrated  from  a(0) 
*  yii(r)  to  a(r)  =  xi  to  obtain  Eq.  15. 


y^jCt)  -  [  (15) 


In  Eq.  15,  =  bj^  -  1  and  is  given  by  Eq.  16. 


Qi  =  (16) 


Qi  in  Eq.  16  is  assumed  to  be  a  power  function  of  the 
maximum  applied  stress  a,  and  ^  and  7  are  constants  to  be 

determined  from  available  f ractographic  data. 

» 

Equations,  for  bi  =  1,  can  be  integrated  from  a(0)  =  yi,-(T) 
to  a(r)  =  xi  to  obtain  y2^^(T)  in  Eq.  17. 


yii(T)  =  exp  (-Q^t) 


Qi  in  Eq.  17  is  given  by  Eq.  16. 

If  applicable  f ractographic  results  are  not  available 
for  the  desired  design  conditions  (e.g.,  load  spectra,  % 
l^olt  load  transfer,  stress  level,  etc.),  an  analytical  crack 
growth  program  [e.g.,  37]  can  be  used  to  generate  the  SCGMC 
needed  to  determine  yii(r)  for  a  given  xi  and  t  [37].  When 
an  analytical  crack  growth  program  is  used,  yii(T)  must  be 
determined  in  a  manner  which  is  consistent  with  that  used  to 
determine  the  EIFS  distribution  from  the  f ractographic  data. 


4. 

stress 
Eq.  18 


Compute  the  probability  of  crack  exceedance  for  each 
region,  i.e.,  p (i ,  t)=P  [a (t)  >  x  ]  =  1-F  [yT(T)].  Use 

and  19  for  b  1  and  for  b  =  l}  respec^i^ely^ 


For  b  ^  1 
p(i,T)  =  1-exp 
=  0 


I'l 

cQ3 

J  j;0<yii(x) 

;yii(T^)=x 

(18) 
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For  b  -  1 


P(1.T>  -  1 


exp 


a 


Q6 


0<yii(T)<x^ 


P(1,t)  «  0 


;  yu(T)»x^ 


(19) 


In  Eq .  18  and  19,  3  is  given  by  Eq.  5  and  is  given  by 
Eqs.  6  and  11.  yii{T)  in  Eqs .  18  and  19  is  given  by 
Eqs.  15  and  17,  respectively,  or  by  a  suitable  analytical 
SCGMC . 


A  flow  diagram  for  determining  y,  .(r)  and  p(i,r)  is 
shown  in  Fig.  5  for  yii{r)  >  0  and  ^for  <  0. 

Equations  20  and  21  are  an  extension  of  Eqs.'^^lS  and  18* 
respectively.  The  negative  EIFS  issue  is  discussed 
subsection  4.2.3. 


in 


5.  The_  average_  number  of  details  N(i,r)  and  the 
standard  deviation  a(i,r)  in  the  ith  stress  region  with  a 
crack  size  greater  than  at  service  time  r  are  determined 
using  the  binomial  distribution  and  are  expressed  as 
follows: 


NCi,T)  =  p(i,T) 

(22) 

% 

aj^(i,T)  =  p(i,T)  [l-p(i,T)]| 

(23) 

in  which  Nj[  denotes  the  total  number  of  details  in  the  ith 
stress  region.  The  average  number  of  details  with  a  crack 
size  exceeding  x]^  at  the  service  time  t  for  m  stress 
regions,  X  (r),  and  its  standard  deviation,  aT(r),  can  be 
computed  using  Eqs.  24  and  25. 


L(t) 


S  N(l,t) 
i=l 


(24) 


(25) 
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Fig,  5  Flow  Diagram  and  Equations  for  Computing  yii(T)  and  p(i. 


Equations  24  and  25  can  be  used  to  quantify  the  extent  of 
damage  for  a  single  detail,  a  group  of  details,  a  part,  a 
component,  or  an  airframe.  Upper  and  lower  bounds  for  the 
prediction  can  be  estimated  using  L(r)  iZa^U),  where  Z  is 
the  number  of  standard  deviations,  aj(r),  from  the  mean, 
17(7).  Equations  22  through  25  are  valid  if  cracks  in  each 
detail  are  relatively  small  and  the  growth  of  the  largest 
crack  in  each  detail  is  not  -affected  by  cracks  in 
neighboring  details.  Hence,  the  crack  growth  accimiulation 
for  each  detail  is  statistically  independent  (Refs.  1,  5, 
21 , 34-36) . 
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SECTION  IV 


DURABILITY  ANALYSIS  DETAILS 


4.1  INTRODUCTION 

Durability  analysis  concepts,  procedures,  and  final 
understandings  evolved  from  the  Phase  II  effort.  This  was 
a  "learning  process."  The  purpose  of  this  section  is  to 
document  the  evolving  understandings  from  the  Phase  II 
effort  and  to  support  the  Durability  Design  Handbook  121]. 
In  any  case.  Ref.  21  is  the  governing  document  for 
durability  analysis. 


4.2  INITIAL  FATIGUE  QUALITY 

Initial  fatigue  quality  (IFQ),  or  EIFS  distribution,  is 
one  of  the  most  important  elements  of  the  durability 
analysis  method  developed.  Four  aspects  of  IFQ  are 
discussed  in  this  section:  (1)  crack  growth  power  law 
variations,  (2)  generic  EIFS,  (3)  negative  EIFS  and  (4) 
details  other  than  fastener  holes.  IFQ  is  also  discussed  in 
Section  3.4  herein  and  elsewhere  [1,2,5,21,31,32,34-36]. 


4.2.1  Crack  Growth  Power  Law  Variations 

Two  variations  of  the  crack  growth  power  law,  Eq.  2, 
were  used  to  develop  equations  for  the  EIFS  cumulative 
distribution,  Fa(0)[x)*  For  example,  F^/q\(x)  based  on  h  * 
1  and  b  *=  1  are  discussed  in  this  section.  ' 

Either  Eq.  8  or  13  can  be  used  to  define  the  IFQ  or  EIFS 
cumulative  distribution.  However,  if  Eq.  8  is  used  and  b  < 
1,  negative  EIFS  values  are  theoretically  possible  (Ref. 
Fig.  5).  The  negative  EIFS  issue  is  discussed  in  subsection 
4.2.3.  In  Eq.  8,  if  b  >  1  all  EIF  values  will  be  greater 
than  zero.  In  Eq.  13,  for  b  =  1  (Case  II)  all  EIFS  values 
will  be  positive. 

Intuitively,  the  two-parameter  form  of  Eq.  2  (i.e.,  Q,  b 
*  1;  Case  I)  should  provide  a  better  overall  fit  to  the 
f ractographic  data  than  the  one-parameter  form  (i.e.,  Q,  b  = 
1;  Case  II).  However,  the  one-parameter  form  is  simpler  and 
it  worked  very  well  for  the  final  demonstration  of  the 
durability  analysis  methodology  [35,36]. 
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As  long  as  consistent  procedures  are  used  to  determine 
the  crack  growth  parameters  in  Eq.  2  for  a  given 
f ractographic  data  set(s),  either  Eq.  8  or  Eq.  13  can  be 
used  to  define  the  IFQ.  If  the  service  crack  growth  master 
curve  (SCGMC)  is  made  compatible  with  the  applicable  EIFS 
master  curve(s),  the  resulting  p{i,T)  predictions  will  be 
virtually  the  same  for  either  variation  of  Eq.  2  (Ref. 
Subsection  4.5.3  and  Table  6). 


4.2.2  Generic  EIFS  Distribution 

An  EIFS  distribution  is  "generic"  if  it  depends  only  on 
the  material  and  manufacturing/fabrication  processes. 
Theoretically,  the  EIFS  distribution  should  be  independent 
of  design  variables,  such  as  load  spectrum,  stress  level, 
percent  bolt  load  transfer,  environment,  etc. 

The  IFQ  model  parameters  in  Eq.  8  and  13  can  be 
determined  for  fastener  holes  using  applicable  f ractographic 
results  for  one  or  more  data  sets.  A  f ractographic  data  set 
refers  to  the  f ractographic  results  (i.e.,  a(t)  versus  t 
values)  for  replicate  fatigue  tests  (e.g.,  same:  material, 
specimen  geometry,  drilling  procedure,  load  spectrum,  stress 
level,  %  bolt  load  transfer,  fastener  type/fit,  etc.). 
Since  the  fractography  reflects  specific  fatigue  test 
variables,  the  derived  EIFS  cumulative  distribution, 
F„/o')(x),  will  not  be  strictly  "generic".  This  situation  is 
further  complicated  by  the  effects  of  scatter  in  the  fatigue 
test  results  on  the  IFQ  model  parameters. 

The  IFQ  model,  described  in  Fig.  3,  was  derived  for  a 
single  TTCI  data  set.  However,  the  model  can  be  applied  to 
one  or  more  f ractographic  data  sets  using  the  data  pooling 
procedures  described  in  the  Durability  Design  Handbook  [21]. 
This  is  a  practical  approach  for  determining  a  suitable  EIFS 
distribution  for  durability  analysis. 

Suppose  two  sets  of  replicate  specimens  are  fatigue 
tested  using  the  same  loading  spectrum  but  different  stress 
levels.  Using  the  f ractographic  results  for  each  data  set, 
the  respective  TTCI's  for  a  given  a  can  be  determined  for 
each  data  set  (Ref.  Fig.  6).  For  the  EIFS  cumulative 
distribution,  Fa(03(x),  to  be  "generic",  the  TTCI 
distribution,  F^Ct)^,  for  data  set  1  and  2  should  transform 
into  the  same  EIFS  distribution  as  shown  in  Fig.  6. 

Theoretically,  a  "generic"  EIFS  distribution  can  be 
obtained  by  forcing  the  IFQ  model  parameters,  3^  ,  af  t>  and 
Q^,  to  be  constants  for  one  or  more  f ractographic  data  sets. 
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Fig.  6  Illustration  of  Generic  EIFS  Condition 
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Using  this  approach,  compatible  IFQ  model  parameters  can  be 
obtained  for  different  f ractographic  data  sets.  This  is 
important  in  order  to  obtain  compatible  EIFS  values  for  each 
data  set  or  sets  and  a  viable  EIFS  distribution. 

It^  appears  that  cx  and  Q/3  are  material  constants 
depending  on  the  fastener  hole  drilling  procedures  (and 
quality)  and  type  of  fastener  hole.  Further  research  is 
needed  to  evaluate  and  compare  QjS  values  for  different 
f ractographic  data  sets  [e.g.,  6  ].  IFQ  model  parameter 
studies  are  documented  in  Appendix  A  and  B  in  this  volume. 
Recommended  procedures  for  determining  the  IFQ  model 
parameters  are  illustrated  in  Ref. 21  . 


4.2.3  Negative  EIFS 

Negative  EIFSs  are  not  physically  possible.  However, 
in  some  cases,  when  f ractographic  results  are  back- 
extrapolated,  negative  EIFSs  can  result.  The  purpose  of 
this  section  is  to:  (1)  explain  the  conditions  under  which 
negative  EIFSs  can  occur  and  (2)  discuss  the  meaning  and 
significance  of  negative  EIFSs  in  terms  of  the  IFQ  model. 

The  IFQ  model  is  a  "mathematical  tool"  for  quantifying 
the  initial  manufactured  state  of  a  structural  detail  in 
terms  of  an  EIFS.  This  "tool"  provides  a  mathematical  basis 
for  making  crack  exceedance  predictions  for  design 
conditions  different  than  the  f ractographic  data  base. 

The  EIFS  master  curve  is  obtained  by  curve  fitting  a 
suitable  crack  growth  law  (e.g.,  Eq.  2)  to  f ractographic 
results.  Using  the  same  f ractographic  results,  different 
EIFS  master  curves  can  be  obtained  for  each  crack  growth  law 
and  curve  fit  criterion  used.  The  EIFS  master  curve  can 
also  be  numerically  defined  using  a  series  of  crack  size¬ 
time  points  rather  than  a  deterministic  equation  [5]. 

The  EIFS  cumulative  distribution  Fa(o)(x)  is  based  on  a 
transformation  of  the  TTCI  distribution,  Eq.  1,  using  the 
applicable  EIFS  master  curve.  A  different  Fa('o')(x)  equation 
is  obtained  for  each  EIFS  master  curve  used.  ''Significantly, 
each  equation  will  be  statistically  compatible  with 
the  TTCI  distribution:  i.e.,  Fo(o)(x)  =  1  -  F-T.(t). 
Whatever  Fa(0)(x)  equation  is  used,  the  resulting  "service 
crack  growth  master  curve"  used  should  be  compatible  with 
the  EIFS  master  curve  for  making  p(i,T)  predictions  (Ref. 
Sections  4.5.2  and  4.5.3). 
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Negative  EIFSs  V7ill  now  be  explained  two  ways  usinp- 
(1)  an  untruncated  EIFS  distribution  and  (2)  a  truncated 
EIFS  distribution.  The  following  discussion  provides  ad¬ 
ditional  insight  into  the  features  and  flexibility  of  the  IFQ 
model  for  durability  analyses. 

4 . 2 . 3 . 1  Untruncated  EIFS  Distribution 

Suppose  an  EIFS  distribution  with  both  positive  and 
negative  EIFSs  is  obtained  as  illustrated  in  Fig.  7.  Assume 
the  resulting  Fg(Q)(x)  distribution  was  determined  following 
the  transformation  procedure  described  for  the  IFQ  model. 

Since  negative  flaws  are  not  physically  possible,  how 
can  the  IFQ  model  be  justified  for  both  positive  and  negative 
EIFSs? ^  The  main  purpose  of  the  IFQ  model  is  to  quantify  an 
EIFS  distribution  that  will  track  the  observable  distribution 
of  crack  sizes  when  grown  outward.  The  fact  that  a  portion 
of  the  EIFS  population  is  negative  may  create  an  emotional 
problem  (i.e.,  negative  flaws  are  not  physically  possible). 
However,  either  positive  or  negative  EIFSs  can  be  handled 
mathematically  by  the  IFQ  model. 

In  Fig.  7,  a  positive,  zero,  and  negative  EIFS  are 
shown.  Using  the  applicable  EIFS  master  curve,  there  is  a 
corresponding  TTCI  for  each  EIFS  value.  Note  that  "negative 
EIFSs"  will  eventually  grow  to  a  crack  size  aQ.  Keep  in 
mind  that  the  IFQ  model  is  simply  a  "mathematical  tool"  for 
making  crack  exceedance  predictions. 

The  positive  EIFSs  support  the  fracture  mechanics 
theory  that  initial  falws  exist  in  a  structural  detail  and 
such  flaws  grow  according  to  fracture  mechanics  theories. 

On  the  other  hand,  the  negative  EIFSs  suggest  that  a  portion 
the  EIFS  population  may  behave  like  an  "initiation  process" 
where  no  initial  flaws  are  considered  present.  In  any  case, 
the  IFQ  model  developed  can  be  used  to  predict  crack  exceedance, 
P(i/T),  as  a  function  of  service  time  whether-or-not  the  EIFS 
distribution  includes  negative  EIFSs  (Ref.  Subsection  4.5.3). 

The  upper  tail  portion  of  the  EIFS  distribution  is  of 
most  interest  for  p(i,T)  predictions  because  the  extent  of 
damage  is  governed  by  the  larger  initial  flaw  sizes.  Since 
the  negative  EIFSs  will  be  in  the  lower  tail  portion  of  the 
^a{0)(^)  distribution,  in  most  cases  negative  EIFSs  would 
not  have  to  be  considered  anyway. 
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Fig.  7  IFQ  Distribution  With  Both  Positive  and  Negative  EIFSs 


4. 2. 3. 2  Truncated  EIFS  Distribution 


The  concept  of  a  truncated  EIFS  distribution  is  used  to 
explain  cases  where  a  portion  of  the  EIFS  population  may  be 
negative.  For  example,  in  Fig.  8  the  EIFS  population 
includes  both  positive  and  negative  values.  The  cross- 
hatched  area,  A^,  under  the  TTCI  distribution  (Fig.  8) 
denotes  the  portion  of  the  TTCI  population  with  TTCIs  s  r, 
or  the  probability  that  T  >  Tj.  Note  that  TTCIs  >  Tj 
result  in  EIFSs  <  0.  The  EIFS  distribution  is 
statistically  compatible  with  the  TTCI  distribution. 
Therefore,  the  area  under  the  EIFS  probability  density 
function,  “  dF^(Q\ (x)/dx,  is  also  equal  to  the 
corresponding  area  A^  under  the  TTCI  distribution  (Fig.  8). 

Assume  that  EIFSs^  <  0  are  inadmissible  in  Fig.  8.  The 
area  A^  under  the  EIFS  probability  density  function  can  be 
replaced  by  a  bar  of  finite  width  with  area  Ai.  The  "bar", 
located  at  EIFS  =  0,  can  be  acknowledged  mathematically 
using  the  dirac  8  function. 

The  EIFS  distribution  shown  in  Fig.  8  is  called  a 
"truncated  EIFS  distribution"  because  the  EIFS's  <  0  are 
treated  as  inadmissible  quantities.  The  truncated  EIFS 
distribution  can  also  be  used  to  represent  the  IFQ  of 
replicate  structural  details. 

In  Fig.  8,  notice  that  the  distribution  of  crack  sizes 
is  truncated  until  Tj  is  reached.  Also  note  that  the  bar 
representing  the  population  of  EIFSs  <  0  decreases  in  area. 
For  example,  the  bar  area  is  A^  at  t  =  0,  A,  at  t  =  Ti,  and 
zero  at  t  =  Tj.  The  percentile  corresponding  to  Tj  can  be 
arbitrarily  selected  by  making  A^  a  selected  finite  quantity 
(e.g.,  0.001). 

The  EIFS  cumulative  distribution,  F  ® 
truncated  EIFS  distribution  is  shown  in  Fig.  9.  Note  that 
F^(0) starts  at  "Ai"  instead  of  0  when  a  truncated  EIFS 
distribution  is  used. 

EIFS  and  TTCI  rankings  are  noted  in  Fig.  10  for  a 
truncated  EIFS  distribution.  Note  that  the  truncated  EIFS 
distribution  is  interpreted  in  the  same  manner  as  an 
untruncated  EIFS  distribution.  Therefore,  the  derived 
equations  for  the  IFQ  model  apply  to  the  truncated  EIFS 
distribution. 
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Fig.  8  Truncated  IFQ  Distribution 
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Fig.  9  Truncated  EIFS  Cumulative  Distribution 
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Fxg.  10  Ranked  EIFS  Crack  Growth  for  Truncated  Distribution 


30 


4.2.4  Details  Other  Than  Fastener  Holes 


The  durability  analysis  method  described  in  Section  III 
theoretically  applies  to  different  types  of  structural 
details  susceptible  to  fatigue  cracking  (e.g.,  fastener 
holes,  lugs,  cutouts,  fillets,  etc.).  Although  the  method 
has  been  demonstrated  for  clearance-fit  fastener  holes, 
further  rsearch  is  required  to  verify  the  method  for  other 
detail  types  and  for  different  combinations  of  details. 

Theoretically,  the  IFQ  for  a  given  detail  type  can  be 
quantified  using  applicable  f ractographic  results  for  the 
desired  detail.  Suitable  specimens  need  to  be  designed  for 
acquiring  fatigue  cracking  data  for  lugs,  fillets,  cutouts, 
etc.  Crack  initiation  data  should  be  generated  for 
different  detail  types  for  different:  •  materials,  stress 
levels,  load  spectra,  manufacturing  techniques,  etc.  These 
data  should  be  used  to  verify  and  refine,  if  required,  the 
durability  analysis  method  for  those  details  which  may  have  a 
significant  effect  on  the  structural  maintenance 
requirements  and  economic  life. 


4.3  DETERMINATION  OF  IFQ  MODEL  PARAMETERS 


An  IFQ  model  was  developed  to  determine  the  key  input 
for  the  durability  analysis  method  developed  -  the  EIFS 


cumulative  distribution 


Because  of  the 


importance  of  F^/qn(x),  the  IFQ  model  was  evaluated 
extensively,  using^  available  f ractographic  results  for 
fastener  holes,  over  the  course  of  the  program.  There  are 
many  ways  that  the  IFQ  model  parameters  can  be  determined 
and  several  numerical  procedures  and  variations  thereof  were 
investigated  [e.g.,  1,5,6,11-19  and  31-26  . 

♦  The  recommended  procedures  and  guidelines  for 
calibrating  the  IFQ  model  parameters,  documented  in  the 
Durability  Design  Handbook  [21],  evolved  from  the  total 
effort  of  this  program.  The  Phase  II  evaluation  of  the  IFQ 
model  was  definitely  a  "learning  process".  The  purpose  of 
this  section  is  to  document,  for  reference  purposes,  the  IFQ 
model  calibration  procedures  used  for  two  Phase  II  studies 
(Ref.  Appendices  A  &  B),  to  discuss  f ractographic  data 
pooling  concepts  and  to  support  Ref.  21.  In  any  case.  Ref. 
21  is  the  governing  document  for  durability  analysis. 


4.3.1  Fractographic  Data  Pooling 


The  IFQ  for  fastener  holes  can  be  quantified  using  the 
IFQ  model  of  Section  III  and  suitable  fractographic  results. 
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IFQ  model  parameters  must  be  consistently  quantified  to  make 
the  parameters  directly  comparable  for  different 
fractographi^c  data  sets.  This  is  very  important. 

A  "fractographic  data  set"  refers  to  the  set  of 
f ractographic  data  obtained  for  replicate  test  specimens 
(e.g.,  same  material,  same  geometry,  etc.),  and  fatigue  test 
conditions  (e.g.,  same  stress  level,  same  spectrum,  etc.). 
Fractographic  data  pooling  techniques  have  been  developed 
and  evaluated  for  quantifying  the  IFQ  model  parameters  on  a 
common  baseline  (Ref.  Section  4.8.1  and  Appendices  A  and  B). 

Pooling  the  fractographic  results  for  different  data 
sets  is  a  practical  approach  for  calibrating  the  IFQ  model 
parameters.  Pooling  the  fractographic  results  effectively 
increases  the  sample  size.  This  not  only  increases 
confidence  in  the  IFQ  model  parameter  values  but  it  also 
ensures  that  the  resulting  values  will  be  "representative 
of  the  overall  fractographic  data  base  used.  Furthermore, 
since  the  IFQ  model  parameters  are  on  the  same  baseline, 
parameter  values  for  different  fractographic  data  sets  can 
be  directly  compared  and  the  effects  of  different  variables 
on  the  parameter  values  can  be  more  readily  evaluated. 
Finally,  fractographic  data  pooling  is  essential  to  justify 
using  the  resulting  IFQ  distribution  for  different  stress 
levels,  load  spectra,  etc.,  than  reflected  in  the 
fractographic  data  base. 


4.3.2  EIFS  Master  Curve  Parameters  (Q^  and  b*) 

EIFS  master  curves,  based  on  two  variations  of  Eq.  2, 
are  given  by  Eqs.  3  and  9  for  a  single  TTCI  data  set.  To 
distinguish  the  crack  growth  parameters  Q  and  b  in  Eq.  2  tor 
each  fractographic  data  set,  Eq.  26  is  used. 


da^tl  ^ 
dt 


(26) 


In  Eq.  26,  Q?  and  b»  are  crack  growth  constants  for  the  ith 
fractographic  data  set.  To  make  the  Qt  and  b*  parameters 
"compatible"  for  different  fractographic  data  sets,  bf  =  bj 

s . h*  *=  constant  must  be  ^imposed.  ^  if  Eq.  26  is  used  for 

a  single  fatigue  crack,  use  Q.  ,  and  b.  .  to  distinguish  values 
for  the  i^^  fractographic  data’ set  and’-^jth  fatigue  crack 
from  the  i^^  data  set. 


Using  the  notations  of  Eq.  26,  the  EIFS  master  curve  for 
>  1  and  b*  =  1  is  given  by  Eqs.  27  and  28,  respectively. 
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(27) 

(28) 


EIFS  =  a(0)  =  (Sq  ^  ^ 

EIFS  =  a(0)  =  a^  exp  (-qJt) 


in  Eq.  27,  -  1. 


The  procedures  described  in  this  section  3re 
in  Appendices  Aand  B  for  different  f ractographic 


illustrated 
data  sets. 


4. 3. 2.1  Case  I  (0^-  b^  ^  1) 


Q*  and  b*  are 
for  a  given  b*  there 
[33].  Therefore,  b*  must  be  a 
EIFS  c\amulative  distributions 
data  sets. 


This  means 
vice  versa 
constant  to  obtain  compatible 
for  different  f ractographic 


cross-correlated  parameters 
is  a  corresponding  and 


and  b’^f  can  be  determined  as  follows; 

^2.  1- 

1.  Transform  Eq.  26  into  a  linear  form  for  least  square 
fit  as  follows; 


In  =  In  Q*  +  b*  In  a(t) 


(29) 


2 .  Compute  b* 
f ractographic  results 
fit  criterion) ; 


for  a  single  crack  or  for 
using  Eq.  30  (based  on  least 


pooled 

square 


b .  =  NExy  -  (^x) (ly) 

1  9  2 
NEx  -  (Ex) 


(30) 


where;  x  =  In  a(t) 


y  =  In  [da(t)/dt] 

N  =  total  nxamber  of  x,y,  pairs  pooled 

If  fractographic  results  for  multiple  data  sets  are  used  in 
Eq.  30,  b  =  b*.  The  subscript  "i"  is  used  to  dis¬ 
tinguish  values  for  individual  data  sets. 


3. 

quantify 
imens . 


Once  bi  has  been  determined,  Eq.  31  can  be  used  to 
for  a  single  specimen  or  for  a  group  of  spec- 


* 

=  exp 


■k 

Ey  -  b. 


N 


.Ex\ 
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(31) 


The  terms  shown  in  Eq.  31  are  defined  as  noted  in  Eq.  30. 

A  conceptual  description  of  the  calibration  procedure  is 
shown  in  Fig.  11.  A  computer  program  has  been  developed  for 
determining  b’^  and  values  for  individual  f ractographic 
data  sets  I  7],  This  program  has  been  tailored  for  pooling 
f ractographic  data  for  multiple  data  sets. 


4. 3. 2. 2  Case  II  (01.  b^  =  1) 

Q*  for  Case  II  can  be  determined  using  Eq.  31  and  b^  = 
1.  can  be  defined  using  f ractographic  results  for  one  or 
more  fatigue  cracks  in  a  f ractographic  data  set. 


4.3.3  TTCI  Parameters 

The  purpose  of  this  section  is  to  document  the  way  the 
Weibull  parameters  (a  and  p)  in  Eq.  1  were  determined  for 
two  Phase  II  modeling  studies  (Appendices  A  and  bK  For 
these  studies,  e  =  0  and  =  bq  were  considered. 

Procedures  are  described  for  determining  a  and  P  for  a 
single  data  set  and  for  pooled  f ractographic  data  sets. 

The  notations  for  a,  p  and  e  in  Eq.  1  are  modified  to 
distinguish  values  for  different  f ractographic  data  sets 
when  the  data  pooling  procedure  is  used.  For  example,  the 
following  notations  have  been  adopted  for  the  Durability 
Design  Handbook  [21]. 

a  =  Weibull  shape  parameter  used  in  Fa(o<x)  based  on  one 
or  more  TTCI  data  sets  and  the  TTCI  data  pooling 
procedure 

=  Weibull  shape  parameter  for  the  ith  TTCI  data  set. 

Pj_  -  Weibull  scale  parameter  for  the  ith  TTCI  data  set, 
including  "statistical  scaling  of  p”  if  applicable 
(Ref.  Section  4.4) . 

Pn  =  Weibull  scale  parameter  for  the  ith  TTCI  data  set 
with  no-, /Statist ical  scaling  of  P” .  Note:  P^  = 
Pz.i  L)  ,  where  =  number  of  fastener  holes 

peV  specimen  where  only  the  largest  fatigue  crack 
in  1  of  ^^holes  per  specimen  is  used  to  define  the 
TTCI  data  set. 
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In  da(t)/dt 


SINGLE  SPECIMEN 
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SINGLE  DATA  SET 


#2 


In  da(t)/dt 


POOLED  DATA  SETS 


Fig.  11  Calibration  Concepts  for  b* 
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ci  »  Weibull  lower  bound  TTCI  for  the  ith  TTCI  data  set 
(Note:  ei  «  0  considered  for  modeling  studies 

documented  in  Appendices  A  and  B). 


4. 3. 3.1  For  One  TTCI  Data  Set 

and  for  a  given  TTCI  data  set  can  be  determined  as 
follows: 

1.  netermine  the  TTCI  values  for  a  selected  reference 
crack  size  (e.g.,  a  =  0.03")  using  the  applicable 
f ractographic  data.  ^ 

2.  Rank  the  TTCI  values  in  ascending  order  and  estimate 
the  cumulative  distribution  of  TTCI,  F^(t),  using  Eq.  32. 

F  (t)  =  ^  (32) 

n+1 


Where  r  =  rank  of  TTCI  value,  t  (1,  2,...,n) 

n  =  number  of  TTCI  values  in  the  data  set  for  given  . 

3.  Transform  Eq.  1  into  the  following  least  squares  fit 
form: 

Z  =  a^Y  +  U  (33) 


where  Z  =  In  {-ln[l-I^  ( t )  ] } 

Y  =  In  t 

U  =  -  a.  In  /3, 

1  t 

P-  =  Pn  (£.)  (Ref.  Subsection  4.4) 


a , 


and  p.  can  be  determined  using  Eqs.  34  and35. 


r^pectively^  In  Eq.  34, 


NSYZ-(EY)  (E^) 

NEY^-(EY)^ 


(34) 
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N  =  total  number  of  Y,Z  pairs; 


Y  and  Z  are  defined  in  Eq.  33. 


(35) 


4. 3. 3. 2  For  Pooled  TTCI  Data  Sets 

A  conceptual  descript  ion  of  the  pooling  procedures  used  to 
determine  a  and  in  the  studies  of  Appendices  A  and  B  is 
shown  in  Figure  12.  Essential  steps  areas  follows; 

1,  Determine  TTCI  values  for  a  selected  reference  crack 
size  (e.g.,  =  0.03")  for  each  f ractographic  data  set. 

2^  Compute  the  average  TTCI  value  for  each  TTCI  data 
set  (X^).  Then,  normalize  the  TTCI  value  in  each  data  set 
by  dividing  each  TTCI  by  the  corresponding  average  TTCI  for 
the  data  set. 

3.  Normalize  TTCI  values  for  each  data  set  to  be 
pooled  and  rank  the  pooled  values  in  ascending  order.  Then, 
estimate  F^(t)  using  Eq.  32. 

4.  The  a  value  for  the  pooled  data  sets  (Fig.  12,  Frame 
B)  can  be  determined  using  Eq.  34  and  the  applicable  Y  and  Z 
values  from  Eq.  33. 

5.  Given  a  (based  on  pooled  TTCI  results),  the 
corresponding  Pi  for  each  data  set  can  be  determined  using 
Eq.  35  (where  a  =  ai)  and  the  TTCI  results  for  each  data  set 
separately.  This  procedure  is  conceptually  illustrated  in 
Frame  C  of  Fig.  12  for  two  TTCI  data  sets. 


4.4  STATISTICAL  SCALING  OF  p 

Laboratory  tests  are  both  expensive  and  time  consuming. 
In  order  to  save  the  cost  and  duration  of  experimental 
tests,  several  fastener  holes,  say  /  holes,  can  be  used  in 
one  coupon  specimen.  The  coupon  test  can  be  stopped  when  a 
visible  crack  appears  in  any  hole.  Consequently,  such  a 
crack  is  the  largest  crack  in  /  holes.  The  time  to  crack 
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Fig,  12  Conceptual  Description  of  Procedures  for  Determining 
^and/^,  for  Pooled  TTCI  Data  Sets 


initiation  (TTCI)  estimated  from  the  largest  crack 
represents  the  smallest  TTCI  in  /  holes,  i.e.,  the  shortest 
TTCI  of  a  sample  of  size  -f. 


As  a  result,  one  has  to  determine  the  parameter  values 
a,P  and  e  appearing  in  the  three  parameter  Weibull 
distribution  of  TTCI,  Eq.  1,  from  a  set  of  coupon  data  which 
represents  the  smallest  value  of  TTCI  of  a  sample  of  size 
In  this  connection,  the  theory  of  minimum  value  of  a  sample 
of  size  /  can  be  used. 

Let  T£  be  the  minimum  TTCI  of  a  sample  of  size  Then, 
the  cumulative  distribution  function  of  Tg  can  be  derived 
from  the  three-parameter  Weibull  distribution  of  TTCI  given 
by  Eq.  1  as  follows: 


(t) 


PlT^_<t]  =  1  -  PlTj^>t] 


1  -  ■{PlT>t3}^  =  1  -  ■  {l-PlT<t3}^ 


(36) 


1  -  Il-F,j.(t)]^ 


Substituting  Eq.  1  into  Eq.  36,  one  obtains  the 
cumulative  distribution  for  Tjj^  as 


(37) 


or 


(38) 


in  which  is  the  scale  parameter  of  the  distribution  of 

the  smallest  TTCI  of  a  sample  of  size  /  given  by 


B 


or 


B  = 


39 


(39) 


It  is  important  to  observe  from  Eq.  38  that  the 
distribution  of  the  smallest  TTCI  of  a  sample  of  size  /  is 
also  three-parameter  Weibull.  Furthermore,  the  shape 
parameter  a  and  the  lower  bound  c  are  identical  to  those  of 
the  TTCI  while  the  scale  parameter  is  related  to  0 
through  Eq.  39. 

Thus,  the  test  data  of  the  smallest  TTCI  can  be  best 
fitted  by  the  three-parameter  Weibull  distribution  [see  the 
procedures  described  in  subsection  4.3.3]  to  estimate  the 
parameter  values  a  and  0  .  Then,  the  scale  parameter,  0,  of 
TTCI  is  estimated  from  0^  using  Eq.  39. 

Equation,  39  can  be  used  to  estimate  0.  from  ,  (Note: 
^  =  /3£.  (-^ r' .  This  important  concept  is  further 
^scussed^  and  illustrated  in  the  Durability  Design  Handbook 
[21] . 


4.5  PROBABILITY  OF  CRACK  EXCEEDANCE 

The  objectives  of  this  section  are;  (1)  show  that 
comparable  p{i,r)  predictions  can  be  obtained  using  either 
variation  of  the  crack  growth  law  of  Eq.  2  (i.e.,  b  =  1  or  b 
#  1),  (2)  illustrate  how  to  determine  a  SCGMC  that  is 
compatible  with  the  EIFS  master  curve,  and  (3)  show  that 
valid  p(i,T)  predictions  can  be  made  using  either  ±  Yi;? 
(EIFS)  values. 


4.5.1  IFQ  Based  on  AFXMR4  Data  Set 

For  illustration  purposes,  the  IFQ  model  parameters  were 
determined  for  the  AFXMR4  data  set.  A  description  of  the 
AFXMR4  data  set  is  given  in  Table  1.  Fractographic  results 
from  Vol.  VIII  [  6]  were  used  and  included  a  crack  size 
range  of  0.030"  -  0.050".  IFQ  model  parameters  were 
determined  for  two  different  variations  of  Eq.  26  and  the 
results  are  referred  to  as  Case  I(b^  *  1)  and  Case  II  (b*  = 


Using  an  upper  bound  EIFS  =  x„  =  0.0275",  the 

corresponding  lower  bound  TTCI's  (e)  for  five  different  Bq's 
(i.e.,  0.030V  0.035V  O.OAO','  0.045V.&0 . 050" )  were  determined.  The 
TTCI’s  for  five  different  crack  sizes  were  normalized  using 
the  average  TTCI  for  each  sq.  Normalized  results  were  then 
pooled  and  a  was  determined  using  the  procedures  described 
in  Subsection  4.3.3.  Using  a,  the  corresponding  /I,  and  the 
corresponding  Q^0ii^  values  were  determined  for  ^%he  five 
different  crack  sizes.  IFQ  model  parameter  results  for  Cases 
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Table  1  Description  of  AFXMR4  Fractographic 

Data  Set 


Material:  7475-T7351  A1 

Type  Specimen:  15%  Load  Transfer 

(Reverse  double  dog-bone) 

Fasteners:  MS  90353  (%”  dia.)  Csk  Rivets 

Hole  Drill:  Modified  Winslow  Spacematic 

Gross  Stress: 

Load  Spectrum: 


<^Q  =  34  ksi 


F-16  400 -hour 

(1  service  life  =  8000  hours) 


Number  of  Specimens; 
Fractography  Basis: 


Largest  crack  for  1  of  4  holes 
per  specimen  (  £  =  4) 
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I  and  II  are  summarized  in  Table  2.  These  results  will  be 
used  in  Sections  4.5.2  and  4.5.3. 


4.5.2  Compatible  SCGMC 

Procedures  are  described  in  Section  3.6  for  making  the 
SCGMC  compatible  with  the  EIFS  master  curve.  These 
procedures  are  discussed  and  illustrated  in  this  section. 

Analytical  predictions  for  crack  growth  in  a  straight -bore  fastener 
hole  with  15%  bolt  load  transfer  in  7475-T7351  aluminum  subjected  to  the 
F-16  500-hour  spectrum  are  shown  in  Table  3.  These  results  [12]  were  de¬ 
veloped  using  an  analytical  crack  growth  program  [37]  .  The  crack  growth 
from  an  initial  flaw  size  of  0.03"  (comer  flaw)  to  0.05"  was  determined 
as  follows.  The  Vlheeler  retardation  model  parameter  was  estimated  using 
applicable  spectrum  crack  growth  data.  A  irultiplying  factor  was  used  to 
adjust  the  stress  intensity  for  the  15%  bolt  load  transfer.  This  factor 
combines  the  effects  of  tension  only  (no  bolt  load  transfer  or  zero  bear¬ 
ing  stress)  and  the  ratio  of  the  bolt  bearing  stress  to  tension  stress. 
Then,  applicable  da/dN  versus  AK  data  were  used  with  a  cycle-by-cycle 
crack  growth  program  and  cycle  counting  scheme  to  make  the  crack  growth 
predictions . 

To  make  the  SCGMC  compatible  with  the  EIFS  master  curve, 
the  same  "bj^"  value  reflected  in  the  EIFS  master  curve  is 
imposed.  For  example,  b^  =  b^  =  0.7372  (Case  I)  and  b  .  =  b* 
=  1.0  (Case  II)  are  used,  along  with  the  a(t)  versus  t  .  data 
from  Table  3,  to  determine  the  Q?  values  for  the 
corresponding  SCGMC.  The  resulting  b.  's^nd  Q.  's  for  the 
two  different  service  crack  growth  Jnaster  curVes  are  shown 
in  Table  4. 

Predicted  a(t)  values  versus  the  analytical  values  are 

summarized  in  Table  5  for  a  =  22.9  ksi  and  o=  28.3  ksi.  The 
SCGMC  parameters  for  the  two  variations  of  the  crack  growth 
equation  (Eq.  26)  are  indicated.  The  predicted  a(t)  values 
were  determined  using  a(t,)  =  0.050"  and  Eqs.  40  and  41. 

Equations  40  and  41  are  generalized  crack  size-time 
relationships  based  on  Eq.  26  for  ^  1  and  =  1, 

respectively  (see  Ref.  21  for  details).  In  Eqs.  40  ana  41, 


(40) 


a(tp  -  a.(.t^y  exp[-Q^(t2-tp]  (4^) 
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Table  2  IFQ  Model  Parameters  for  AFXMR4 
Data  Set 


Case 

’^i 

O^xio'^ 

a 

Ave. 

a 

X 

u 

% 

I 

0.7372 

0.788 

0.617 

1.389 

0.0275 

4 

II 

1.0 

1.854 

1.451 

1.389 

0.0275 

4 

Note:  c*  =  b*  -1 
1  1 


Table  3  Analytical  Crack  Growth  Predictions 
for  7475-T7351  Aluminum  and  the  F-16 
500  Hour-Spectrum 


0 

t 

(KSI) 

a(t) 

(Fit .Hrs . ) 

22 

.9 

0.0300 

0 

0.0325 

321 

0.0364 

821 

0.0404 

1321 

0.0444 

1821 

0.0484 

2321 

0.0500 

2518 

28 

.3 

0.0300 

0 

0.0365 

313 

0.0469 

813 

0.0500 

959 

Ref.  12 
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Table  4 


Summary  of  SCGMC  Parameters 


a 

(KSI) 

Case  I 

Case  II 

b. 

1 

Q.xlO^ 

b. 

1 

Q. xlO^ 

22.9 

0.7372 

0.8586 

1.0 

2.0052 

28.3 

0.7372 

2.2464 

1.0 

5.2333 

Note :  C-»  =  b .  -1 
1  1 
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Table  5 


Comparison  of  a(t)  Predictions  (SCGMC) 
Versus  Actual  a(t)  Values 


cr 

(KSI) 

t 

(Fit,  Hrs,) 

a(t) 

Analytical 

Predicted 

Table  3 

Eq.  41 

Eq,  40 

2: 

>.9 

0 

0,0300 

0.0302 

0.0301 

321 

0,0325 

0.0322 

0.0322 

821 

0,0364 

Cf.0356 

0.0358 

1321 

0,0404 

0.0393 

0.0396 

1821 

0,0444 

0.0435 

0.0437 

2321 

0.0484 

0.0481 

0.0482 

2518 

0.0500 

0.0500* 

0.0500* 

1*0 

0.7372 

-bi 

2.0052 

0.8586 

a'xio^ 

28 

.3 

0 

0.0300 

0.0303 

0.0302 

313 

0.0365 

0.0357 

0.0358 

813 

0.0469 

0.0463 

0.0465 

959 

0.0500 

0.0500* 

0.0500* 

1.0 

0.7372 

1 

5.2333 

2.2464 

Note 

^(^2)  "  0,0500”  (starting  point  for  predictions) 
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a(ti)  and  a(tj)  denote  the  crack  size  at  time  t^  and  t,, 
respectively;  c^  =  -  1  and  is  a  crack  growth 
parameter. 

In  this  case,  there  is  excellent  agreement  between  the 

analytical  and  predicted 'a  (t)'  values .  Since  the  "b"  value  for 

the  EIFS  master  curve  is  imposed  on  the  SCGMC,  the  goodness- 
of-fit  of  Eq.  26  depends  on  the  results  of  the  analytical 
crack  growth  program  (e.g.,  Table  3).  Further  work  is 
required  to  evaluate  "compatible  SCGMC s",  based  on  the 
procedure  described,  for  different  durability  design 
conditions . 


4.5.3  Crack  Exceedance  Predictions 


It  will  be  shown  in  this  section  that  if  the  SCGMC  is 
compatible  with  the  applicable  EIFS  master  curve:  (1)  it 
doesn't  make  a  lot  ^f  difference  which  variation  of  Eq.  14 
is  used  for  p(i,T)  predictions  (i.e.,  Eq.  18  or  19),  and  (2) 
comparable  p{i,T)  predictions  can  be  obtained  for  both 
positive  and  negative  EIFS  or  values. 

Crack  exceedance  analysis  results  are  presented  in  Table 
6  for  two  different  stress  levels.  Results  are  presented 
for  two  different  IFQ  distributions  for  the  AFXMR4  data  set 
and  for  the  corresponding  service  crack  growth  master 
curves.  The  crack  exceedance  predictions  were  based  on  the 
IFQ  model  parameters  of  Table  2  and  the  SCGMC  parameters  of 
Table  4.  The  equations  used  to  compute  yi^Cr)  and  p(i,T) 
have  been  appropriately  modified  for  yi^Cr)  <  0  (see  Fig.  5). 

For  illxostration  pxjrposes,  different  and  rvalues  were  asstjmed  to  make 
p(i,r)  predictions. 


In  Table  6,  the  EIFS's  for  a  given  x^  and  t 
for  two  different  cases. 


for  two  different  cases.  Although  the  yxi(r)  values  ar^ 
different  for  Cases  I  and  II,  virtually  the  same  p(i,T)  i: 


are  shown 
are 
is 

obtained  for  a  given  Xj,r.  yii(r)  values  <  0  resulted  for 
Case  I  in  two  instances  for  a  =  22.9  ksi  and  28.3  ksi.  For 
example,  for  Xi  =  0.030",  r  =  20000  flight  hours  and  a  = 
22.9  ksi,  y-|  i  ( t )  =  -0.0000144"  and  0.000544"  were  obtained 
for  Cases  I  an^II,  respectively.  As  shown  in  Table  6, 
comparable  p(i,T)  values  were  obtained  (e.g.,  0.635  versus 
0.630)  for  two  different  EIFS's  although  the  EIFS  was 
negative  for  Case  I  and  positive  for  Case  II.  If  p(i,T) 
predictions  are  made  using  an  IFQ  distribution  and  a 
compatible  SCGMC,  both  ±  EIFS  can  be  mathematically  handled. 


Further  work  is  required  to  evaluate  the  accuracy  of 
p(i,r)  predictions  for  different  durability  design 
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Table  6  Crack  Exceedance  Analysis  Based  on 
Ccrapatible  SCGSyiC 


cr 

Case 

T 

KFS  -  y^j_(r) 

p(i,r) 

(KSI) 

® ® 

(Inch) 

(Fit,  Hra.) 

(Inch)  (2) 

(3) 

22 

!.9 

1 

0.030 

i 

20000 

.  J 

-0.0000144 

0.635 

II 

0.000544 

0.63O 

I 

0.040 

■  4000 

0.016288 

0.048 

II 

i 

0.017936 

0.045 

I 

0.050 

i 

6000 

1 

0.013045 

0.074 

II 

0.015013 

0.072 

28 

.3 

^  I 

0.030 

4 

15000 

i 

-0.065027 

0.926 

II 

0,0000117 

0.923 

I 

0.040 

i 

4000 

4 

0.0019125 

0.278 

II 

0.0049311 

0.271 

I 

0.050 

1 

6000 

1 

0.0001619 

0.427 

II 

i 

i 

0.0021642 

0.420 

Notes ; 
®  Ref. 
®  Ref. 
®  Ref. 


Table  2  for  IFQ  Model  Parameters 

Table  4  for  SCGMC  Parameters 

Fig,  5  for  Appropriate  Equations  Used 


/ 
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conditions  and  for  different  crack  size  ranges.  Also,  the 
one  and  two  parameter  forms  of  Eq.  14  should  be  investigated 
further  to  justify  using  the  b.  «=  1  form  for  different 
durability  analysis  applications.  ^ 


4.6  FORMATS  FOR  PRESENTING  DURABILITY  ANALYSIS  RESULTS 

The  basic  objective  of  the  durability  analysis  is  to 
analytically  assure  design  compliance  with  the  U.S.  Air 
Force's  durability  requirements  [8-lo].  Durability 
analysis  results  can  be  presented  in  four  basic  formats  as 
follows;  (1)  probability  of  crack  exceedance,  (2)  extent  of 
structural  damage,  (3)  structural  damage  rate  and  (4)  cost 
ratio;  repair  cost/replacement  cost.  The  probability  of 
crack  exceedance  is  the  fundamental  output  of  the  durability 
analysis. 

Durability  analysis  results  can  be  plotted  in  various 
forms  as  shown  in  Fig.  13  .  Extent  of  damage  illustrations 
are  presented  in  the  Durability  Design  Handbook  [21  ]  and 
elsewhere  [1,5,34-36]. 


4.7  DEMONSTRATION  OF  DURABILITY  ANALYSIS  METHOD 

The  durability  analysis  method  developed  has  been 
demonstrated  for  a  full-scale  fighter  structure  (F-16  lower 
wing  skin)  and  for  a  complex  splice  subjected  to  a  B-1 
bomber  spectrum.  Results  are  siammarized  in  the  Durability 
Design  Handbook  [21 1  and  details  are  given  elsewhere  [34- 
36].  The  durability  analysis  demonstrations  for  both 
fighter  and  bomber  load  spectra  are  discussed  in  this 
section. 


4.7.1  F-16  Lower  Wing  Skin 

A  preliminary  and  final  demonstration  of  the  durability 
analysis  methodology  is  documented  in  Ref.  34  and  36, 
respectively.  Both  analyses  were  based  on  coupon 
f ractographic  results  for  fatigue  cracking  in  fastener 
holes.  No  scale-up  of  coupon  data  to  full  size  structure 
was  considered. 

The  preliminary  analysis  reflected  the  following; 

o  IFQ  of  the  countersunk  fastener  holes  was  based 
on  f ractographic  results  for  straight  bore 
fastener  holes  [41]. 
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Crack  Length, 


Crack  Length, 


Crack  Length,  x^ 


Flight  Hours,  T 


"  Service  Life 

cr  ^  Stress  Level  (Fixed  Value) 


(T"  “  Stress  level 
T*  =  Service  Life  (Fixed  Value) 


P  Exceedance  Probability 
T  &  cT are  fixed  values 


For  Given  xj^ 


Average 
Percentage 
Of  Crack 
Exceedance 


Average 
Percentage 
Of  Crack 
Exceedance 


or 


or 


Ave,  No,  Of 
Details  With 


Ave,  No,  Of 
Details  With 


Percentage 
Of  Crack 
Exceedance 


Average 
Percentage 
of  Crack 
Exceedance 


or 

No*,  Of 
Details 
With 
a  ^  X- 


or 


Ave.  No.  Of 

Details  With 

> 

a  *  x^ 


Fig.  13  Formats  for  Presenting  Extent  of  Damage  Results 
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Fig.  13  Formats  for  Presenting  Extent  of  Damage  Results 
(Continued) 
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o  IFQ  model  parameters  were  based  on  a  sinqle 
f ractographic  data  set  (XWPF)  from  Ref.  4i 
and  e  =  0  was  used. 

o  A  crack  growth  power  law,  Eq.  2,  for  the 
EIFS  master  curve  was- used  where  b  1. 

o  "Service  crack  growth  master  curves"  were  used 

based  on  Ref.  37. 

o  The  F-16  lower  wing  skin  was  divided  into  three 

stress  regions. 

The  final  demonstration  of  the  durability  analysis 
method  reflected  the  following: 

o  Fractographic  results  for  countersunk  fastener 
holes  and  data  pooling  procedures  were  used  to 
define  the  IFQ  model  parameters  for  the  counter¬ 
sunk  fastener  holes  in  the  F-16  lower  wing  skin, 
c  9*  0  was  considered. 

o  The  F-16  lower  wing  skin  was  divided  into  ten 
stress  regions  (see  Fig.  14). 

o  A  crack  growth  power  law,  Eq.  2  (with  b  =  1)^ 
was  used  for  the  EIFS  master  curve. 

o  Average  a  and  Q+^.  values  were  used  for  crack 
exceedance  predictions. 

o  The  service  crack  growth  master  curve  was  based 
on  pooled  fractographic  results  and  the  applicable 
stress  level  for  each  stress  region. 

Tear-down  inspection  results  from  the  F-16  wing  durability 
test  article  (after  r  =  16000  hours  of  testing  using  the  F- 
16  500-hour  spectrum)  were  compared  with  extent  of  damage 
predictions.  Excellent  correlations  and  trends  were 
obtained. 


4.7.2  Complex  Splice  (B-1  Bomber  Spectruir) 

A  durability  analysis  of  the  complex  splice  shown  in 
Fig.  15  was  performed  to  demonstrate  the  methodology  for  a 
bomber  load  spectrum.  Details  of  this  analysis  are 
documented  in  Ref.  35.  The  durability  analysis  and  basic 
conclusions  are  briefly  discussed  in  this  section. 
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Fig .  14 


Stress  Zones  for  F-16 


Lower  Wing  Skin 
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Figo  15  Complex-Splice  Specimen 
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Several  complex  splice  specimens  were  fatigue  tested 
using  a  B-1  bomber  load  spectrum,  A  tear-down  inspection 
and  f ractographic  evaluation  were  made  to  determine  the 
number  of  fastener  holes  in  the  outer  row  of  the  splice  with 
a  crack  size  ^  0.050"  at  t  =  13500  flight  hours.  The  extent 
of  damage  was  predicted  for  the  outer  row  of  fastener  holes 
and  the  results  were  compared  with  test  results. 

The  durability  analysis  reflected  the  following: 

o  A  simplified  stress  analysis  was  used  to  estimate  the 
gross  stress  on  the  cross  section  at  the  outer  row 
of  fastener  holes  at  the  interfaces. 

o  The  IFQ  was  based  on  f ractographic  results  for 
countersunk  fastener  holes  in  reverse  dog-bone 
specimens  (15%,  30%  and  40%  load  transfer). 

o  Nine  f ractographic  data  sets  were  pooled  to  calibrate 

the  IFQ  model  parameters.  The  f ractographic  crack  size 
range  used  was  0.005"  -  0.10". 

o  an  =  X  =  0.050" 
u  u 

o  c  =  0;  /  =  4 

o  The  EIFS  master  curve  was  based  on  Eg.  9  (b  =  1) 

o  The  SCGMC  was  based  on  Eqs .  14,  16  and  17  with  b  =  1 

In  general,  the  predicted  results  and  trends  compared 
reasonably  well  with  test  results.  Further  details  are 
given  elsewhere  [21,35]. 


4.8  ADDITIONAL  STUDIES 

Many  studies  were  performed  during  Phase  II  to  evaluate 
the  durability  analysis  methodology  developed  in  Phase  I. 
The  most  important  studies  are  discussed  in  this  section. 
Miscellaneous  durability-related  topics  are  also  discussed. 


4.8.1  IFQ  Model  Studies 

Initial  fatigue  quality  (IFQ)  is  one  of  the  most 
important  elements  of  the  durability  analysis  methodology 
developed.  For  this  reason,  several  studies  were  made  in 
Phase  II  to  evaluate  the  IFQ  model  using  the  wealth  of 
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f ractographic  data  available  for  both  straight-bore  [41]  and 
countersunk  fastener  holes  [6  ], 

The  main  objectives  of  these  studies  were  to: 

1.  Develop,  evaluate,  and  refine  the  procedures  for 
quantifying  the  IFQ  model  parameters  (i.e.,  a  ,  /3^,  e,  Q* , 
bl,  x^)  on  a  comparable  baseline  for  different  f ractographic 
data  sets. 

2.  Evaluate  the  IFQ  model  and  goodness-of-f it  using 
f ractographic  results  for  different  design  and  test 
variables . 


3.  Determine  the  IFQ  model  parameter  values,  including 
the  range  and  variance  of  values,  for  selected  f ractographic 
data  sets. 

4.  Evaluate  the  sensitivity  of  the  IFQ  model  parameter 
values  for  different  design  and  test  variables. 

5.  Compare  the  IFQ  model  parameter  values  for  different 
fractographic  data  sets  for  consistency  and  trends. 

6.  Justify  using  the  same  IFQ  distribution  for  a  given 
material,  fastener  hole  type,  drilling  procedure  and 
fastener  type/fit  for  different  stress  levels,  load  spectra 
and  percentage  of  fastener  load  transfer. 

7.  Develop  confidence  in  the  use  of  the  IFQ  model  for 
different  variables  and  applications. 

The  final  conclusions  and  recommendations  for  IFQ,  given 
in  the  Durability  Design  Handbook  [21  ],  were  based  on  the 
"lessons  learned"  over  the  course  of  the  program. 
Therefore,  the  IFQ  model  procedures  and  results  for  some 
studies  are  not  in  complete  accord  with  the  final 
understandings  reached. 

Two  IFQ  model  studies  are  docximented  in  Appendices  A  and 
B  for  straight  bore  and  countersunk  fastener  holes, 
respectively.  These  studies  were  based  on  fractographic 
results  from  Ref.  41  and  6  . 

The  studies  documented  in  Appendices  A  and  B  reflected 
the  following: 

1.  Equation  1  was  used  for  the  TTCI  distribution  with  e 

=  0. 
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2.  The  crack  growth  law  of  Eq.  2  was  used  with  b  1 
(Case  I). 

3.  <  =  0  and  aiQ  =  were  used. 

Based  on  the  final  understanding  for  the  program,  the 
following  was  recommended:  (1)  use  the  three-parameter 
Weibull  distribution  for  TTCI  with  e  2:  0,  (2)  define  the 
upper  bound  EIFS  value,  in  conjunction  with  e,  (3)  use 
the  two-parameter  form  of  Eq.  2  (i.e.,  Q  and  b)  for  b  >  1; 
otherwise,  use  the  one-parameter  form  (i.e.,  Q  and  b  =  1). 
Further  work  is  required  to  evaluate  the  one-parameter  form 
for  different  design  variables. 

The  following  conclusions,  based  on  the  studies 
documented  in  Appendices  A  and  B,  are  consistent  with  the 
final  understandings  of  the  program  [21]; 

1.  IFQ  model  parameters  should  be  determined  using  the 
f ractographic  data  pooling  procedure  in  order  to  make  the 
resulting  parameters  compatible  for  different  f ractographic 
data  sets. 

2.  a  and  b|  appear  to  be  material  constants  independent 
of  load  spectra,  stress  level  and  percentage  of  shear  load 
transfer. 

3.  a  and  b’^  depend  on  the  type  of  fastener  hole.  For 
example,  typical o?"  values  were  approximately  4.4  and  2.5  for 
straight  bore  holes  and  countersunk  fastener  holes, 
respectively,  b*  values  were  approximately  0.88  and  0.61 
for  straight  bore  and  countersunk  fastener  holes, 
respectively.  Note;  The  one-parameter  form  of  Eq.  2  (b  = 
1)  is  recommended  for  quantifying  the  IFQ  if  bj  <  1. 

4.  EIFS  values  <0  are  theoretically  possible  using  the 
IFQ  model  when  b  <1.0.  This  will  not  create  a  durability 
analysis  problem  as  long  as  the  IFQ  model  is  considered  as  a 
"mathematical  tool"  and  the  SCGMC  is  compatible  with  the 
EIFS  master  curve  (Ref.  Section  4.2.3). 

5.  QIP±  appears  to  depend  on  the  material  and  fastener 
hole  type.  However,  Qi/^i  seems  to  be  insensitive  to  stress 
level,  load  spectra  and  %  shear  load  transfer. 

6.  The  following  empirical  relationships  look  very 
promising  for  durability  analysis  applications: 


(42) 
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7.  and  a  provide  the  common  denominator  for  linking 
different  I ractographic  data  sets  on  a  comparable  basis  and 
the  keys  for  generic  EIFS. 

8.  and  b*  a^e  cross  correlated  constants.  For 
example,  tor  a  given  b.  value  there  is  a  corresponding  QA 
value  and  vice  versa.  Careful  attention  must  be  given  to 
the  and  b?'  values  used  to  define  IFQ. 

9.  In  general,  f ractographic  results  for  straight  bore 
fastener  holes  exhibited  less  scatter  than  those  for 
countersunk  fastener  holes  [  6,41  ]. 


4.8.2  Statistical  Crack  Growth  Studies 

Statistical  crack  growth  studies,  performed  in  Phase  II 
of  the  program,  are  documented  in  Ref.  33  .  These  studies 
were  performed  to:  (1)  better  understand  IFQ  model 
parameters,  and  statistical  variability,  (2)  investigate 
possible  relationships  between  parameters  Q  and  b,  and  (3) 
study  crack  exceedance  predictions  based  on  probabilistic 
and  deterministic  considerations. 

Two  f ractographic  data  sets  from  the  FHQ  program  [41] 
wereused:  WPF  and  XWPF.  Both  data  sets  vere  for  7475-T7351 
aluminum  and  straight  bore  fastener  holes.  The  WPF 
specimens  were  designed  for  no  load  transfer  and  these 
specimens  were  subjected  to  the  F-16  400-hour  spectrum.  The 
XWPF  specimens  were  designed  for  15%  shear  load  transfer 
through  the  fasteners  and  these  specimens  were  subjected  to 
the  F-16  400-hour  spectrum.  These  studies  were  based  on  Eq. 
2  with  b  1. 

Essential  conclusions  from  these  studies  are: 

1.  Crack  growth  rate  parameters  Q  and  b  can  be 
characterized  by  a  log  normal  distribution  reasonably  well. 

2.  Since  the  economic  repair  crack  size  limit  is  usually 
small,  the  influence  of  the  crack  growth  variability  on  the 
extent  of  damage  predictions  and  economical  life  may  not  be 
very  significant. 

3.  Since  the  EIFS  values  are  determined  using  the  TTCI 
data,  such  values  indirectly  account  for  the  crack  growth 
variability.  Hence,  the  initial  flaws  can  be  grown 
deterministically  to  determine  the  probability  of  crack 
exceedance. 
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4.  Similar  studies  should  be  made  using  Eq.  2  and  b  »  1. 


4.8.3  Related  Durability  Studies 
The  following  studies  were  also  performed  during  Phase 


1.  Investigate  effect  of  scale-up  and  hole  interaction 

on  crack  growth. 

2.  Evaluate  method  for  "scaling"  the  crack  growth  for 

different  load  spectra  (development  of  transfer 
functions) . 

3.  Parameter  sensitivity  analysis. 

These  studies  were  exploratory  and  limited  in  scope. 
Details  are  documented  in  Appendix  C. 
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SECTION  V 


DURABILITY  TEST  PROGRAM 


5.1  INTRODUCTION 

An  overview  of  the  durability  test  program  is  given  in 
this  section.  Essential  elements  of  the  test  program  are 
described  and  discussed, including  the  test  matrix,  specimen 
details,  experiment!  methods  and  data  acquisition.  Details 
of  the  test  program,  raw  f ractographic  data  and  other 
specifics  are  documented  in  Volume  VIII  [  6]. 

A  total  of  701  specimens  were  manufactured  and  tested. 
Variables  included  in  the  test  program  were  material, 
spectrum,  stress  level,  fastener  diameter/type,  load 
transfer  and  specimen  geometry  (width  and  thickness). 
Fractographic  data  were  obtained  for  over  800  fatigue  cracks 
and  a  tear-down  inspection  of  the  F-16  lower  wing  skin  from 
the  durability  test  article  was  performed. 

The  basic  objective  of  the  test  program  and  the 
fractographic  data  acquired  was  to:  (1)  provide  the  data 

base  needed  to  evaluate,  refine  and  demonstrate  the  durability  analysis 
methodology  developed  and  (2)  provide  fractographic  data  for  quantifying 
the  initial  fatigue  quality  of  fastener  holes. 


5.2  TEST  MATRIX 

The  durability  program  test  matrix  is  presented  in  Table 
7  and  is  further  detailed  in  Table  8.  Different  specimen 
configurations  were  used  as  noted  in  Figure  16.  Data  set 
designations  and  specimen  details  are  described  in  Sections 

5.3  and  5.4,  respectively.  Refer  to  Volume  VIII  for  further 
details. 


5.3  DATA  SET  DESIGNATIONS 

_  Test  specimens  for  a  given  test  series  were  identified 
using  the  coding  scheme  shown  in  Table  9.  Data  sets  were 
coded  to  identify  the  key  parameters  for  a  given  test 
series.  Each  data  set  is  identified  by  six  characters 
where  each  character  represents  a  different  test  parameter' 
as  shown  in  Table  9.  The  first  symbol  identifies  the 
material  used,  the  second  identifies  the  spectrum,  etc. 
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Table  7  Durability  Test  Matrix 
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Table  8  Teat  Matrix  Detalla 
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Table  8  Test  Matrix  Details  (Continued) 
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Table  9  Data  Set  Designations 


ORDER  OF 
SYMBOL 

TEST  1 

PARAMETER  | 

POSSIBLE  1 

SYMBOLS  1 

1 

MATERIAL  I 

A:  7475-T7351  Al.  | 

T:  2024-T851  Al.  I 

S :  D6hc  steel  I 

2 

SPECTRUM 

F:  F-16  400  HR  BLOCK  I 

B:  B-1  BOMBER  | 

D:  F-16  500  HR  RANDOM  I 

3 

LOAD  1 

TRANSFER 

:  NO  LOAD  TRANSFER  I 

X:  157.  LOAD  TRANSFER  | 
Y:  307.  LOAD  TRANSFER  I 
Z:  407.  LOAD  TRANSFER  I 

4 

STRESS 

LEVEL 

1  L:  LOW  1 

M:  MEDIUM 

H:  HIGH  J 

5 

FASTENER 

TYPE 

1  R:  MS-90353  RIVET  I 

C:  NAS  1580  C’SUNK  BOLT 
P:  NAS  6204  P.H.  BOLT 

6 

FASTENER 

3:  3/16  ! 

DIAMETER 

4:  4/16”  or  1/4” 

NOTES: 

1, 


2. 


Lower  case  letters  for  stress  level  synibols  represent 
a  stress  level  that  is  lower,  but  in  the  same  range 
as  the  upper  case  symbol. 

CBSL:  Complex  Splice/Lov7  Stress  Level 
CBSH:  Complex  Splice/High  Stress  Level 


EXAMPLE: 


ABXHR3 

Material: 

A: 

7475-T7351  Al. 

Spectrum; 

B: 

B-1  BOMBER 

Load  Transfer; 

X: 

157. 

Stress  Level: 

H: 

HIGH  C40.8  KSI) 

Fastener  Type; 

R: 

ms-90353  rivet 

Fastener  Dia.: 

3: 

3/16” 
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Lower  esse  letters  for  the  stress  level  indicste  3  stress 
level  that  is  lower  than  an  upper  case  letter,  but  is  in  the 
same  general  range.  Letters  in  parentheses  at  the  end  of  a 
designation  indicate  one  of  the  specimen  configurations 
shown  in  Figure  16.  Other  peculiarities  of  the  naming 
system  are  indicated  in  Table  9. 


5.4  Load  Spectra 

Three  load  spectra  were  used;  the  F-16  400-hour  wing 
root  bending  moment  spectrum,  the  F-IS  500-hour  wing  root 
bending  moment  durability  spectrum  and  the  B-1  bomber  lower 
carry-through  box  bending  moment  spectrum.  The  F-16  400- 
hour  block  spectrum  used  was  the  same  spectrum  used  for  the 
"Fastener  Hole  Quality"  program  [4l].  Several  specimens 
were  tested  using  the  F-16  500-hour  block  spectrum,  which  is 
the  same  spectrum  used  to  test  the  F-16  wing  component  for 
the  full-scale  durability  demonstration. 


5.5  SPECIMEN  DETAILS 

Five  specimen  designs  were  used  for  testing  the  7475- 
T7351  aluminum  material,  consisting  of  four  different 
amounts  of  load  transfer  and  a  complex  splice.  The  details 
for  the  0%,  15%,  30%  and  40%  load  transfer.  and  complex 
splice  specimens  are  shovm  in  Figs.  17.-20  and  15,  respective  y. 


Most  of  the  0%  load  transfer  specimens  tested  were 
0.375"  thick  as  shown  in  Figure  17.  Specimens  with  a  0.250" 
tfhickness  were  also  tested.  These  specimens  were  identified 
with  a  (B)  after  the  data  set  designation  (Ref.  Table  8). 

Two  different  widths  were  used  for  the  15%  load  transfer 
specimen.  Most  of  the  15%  load  transfer  specimens  tested 
had  a  width  of  1.500"  as  shown  in  Figure  18.  However,  ten 
specimens  with  a  2.000"  width  were  also  tested.  These 
specimens  were  designated  as  the  "A3XMR4(B)"  data  set. 

The  "complex-splice”  specimen,  shown  in  Fig.  1-) ,  sim¬ 
ulates  a  fuselage  splice.  This  specimen  was  made  from  7475- 
'J'7351  aluminum  plate  (1/2"  thickness)  and  included  blind, 
pull-through  rivets  (MS90353-06  (3/16"  dia.)).  No  faying 
surface  sealant  was  used. 

The  07o  load  transfer  specimen  design  shown  in  Figure  21 
was  used  for  testing  the  D6ac  steel  material. 
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Note; 


Two  different  spacings  used:  S  -  1.00" 

S  -  3.00" 


T 

0J50'* 


Fig.  17  No-Load  Transfer  Specimen  Design 
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Fig.  18  157„  Load  Transfer  Specimen  Design 
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Fig,  19  307o  Load  Transfer  Specimen  Design 
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002  TfP 


Fig.  20  407o  Load  Transfer  Specimen  Design 


Fig,  21  D6ac  Steel  Specimen  Design 


TYP. 


70 


All  fastener  holes  for  the  test  specimens  were  drilled 
using  the  modified  Winslow  Spacematic  drill.  Fasteners  were 
installed  using  typical  aircraft  production  procedures.  All 
load  transfer  specimens  were  match-drilled  in  pairs. 

The  load  transfer  specimens  (i.e.,  15%,  30%,  and  40%) 
were  designed  assuming  equal  longitudinal  deformation  for 
the  upper  and  lower  dog  bones  between  lug  pin  and  fastener. 
Pins  and  fasteners  were  assumed  to  be  rigid  and  to 
"perfectly  fit"  their  mating  holes.  Due  to  normal  hole  size 
variations  and  clearance  fits,  the  actual  percentage  load 
transfer  for  each  specimen  varies  depending  on  the  load 
level  required  to  take  up  the  slop  between  the  hole  and  the 
fastener.  A  free  body  diagram  is  shown  in  Fig.  22 . 


5.6  TEST  SETUP  AND  PROCEDURES 

Spectrum  fatigue  tests  were  performed  at  the  GD/FWD  on 
servo-controlled,  hydraulically-actuated  load  frames.  A 
special  hardware  interface  was  used  to  monitor  each  load 
frame  and  to  assure  proper  load  control.  This  system 
provided  a  permanent  record  of  test  events. 

Specimen  dimensions  were  recorded.  Then,  the  required 
load  for  the  desired  stress  level  was  determined.  The  hole 
diameters  and  countersink  depths  were  also  measured  and 
recorded.  Results  are  reported  in  Volume  VIII  [6]. 

Each  load  transfer  specimen  contained  four  fastener 
holes  in  the  test  section.  The  load  transfer  through  the 
fastener  affects  the  stress  on  the  cross-section  of  each 
dogbone.  Therefore,  careful  records  were  kept  on  hole 
positions  in  the  specimens  and  the  size  of  the  largest  crack 
in  each  hole.  Results  are  documented  in  Volume  VIII  [  6  ]. 

Specimens  were  then  loaded  in  the  test  frames  and  left 
to  run.  When  the  specimens  had  fatigued  for  the  desired 
length  of  time,  the  unbroken  dogbones  were  pulled  apart  and 
the  failure  load  was  recorded.  If  the  specimen  failed  in 
the  load  frame,  no  failure  load  was  recorded.  The  broken 
dogbones  were  then  cut  up  and  the  fatigue  cracks  in  each 
hole  were  exposed  and  numbered.  The  D6ac  steel  specimens 
were  sprayed  with  a  clear  lacquer  to  protect  the  fracture 
surface  from  rust.  The  final  sizes  of  all  fatigue  cracks 
were  measured  and  the  largest  crack  for  each  specimen  was 
read  f ractographically .  Many  specimens  had  more  than  the 
largest  crack  read  and  these  data  are  also  presented  in 
Volume  VITI. 
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Fig.  22  Load  Transfer  Specimen  Free  Body  Diagram 


72 


Details  of  the  fractography  set-up  and  procedures  are 
detailed  in  Section  VI  and  in  Volume  VIII. 

The  complex  splice  specimen,  shown  in  Figure  15,  was 
tested  without  lateral  support  using  the  same  test  frames 
used  for  the  other  tests.  This  same  specimen  design,  for  a 
different  aluminum  alloy  (2124-T851),  had  been  previously 
tested  with  and  without  lateral  support  [42].  Results  of 
the  strain  survey  indicated  negligible  effect  of  the  lateral 
support  on  the  resulting  strains  [42]. 
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SECTION  VI 


EVALUATION  OF  TEST  AND  FRACTOGRAPHIC  RESULTS 
6 . 1  INTRODUCTION 


The  basic  objective  of  this  section  is  to  present,  compare 
and  evaluate  selected  test  and  fractographic  results  from  Vol¬ 
ume  VIII  [6  ].  A  comprehensive  evaluation  of  "all  the  data"  in 
Volume  VIII  was  beyond  the  scope  of  this  program.  The  evalua¬ 
tion  of  the  test  and  fractographic  results  herein  is  limited 
to  selected  data  sets.  Further  evaluation  of  the  data  in  Volume 
VIII,  including  a  statistical  analysis,  is  recommended  to  better 
understand  the  effects  of  the  test  and  specimen  variables  on 
the  time-to-crack-initiation  CTTCI)  and  the  time-to-failure  C.TTF)  . 

The  main  objectives  of  the  Phase  II  test  program  were  to: 

1.  provide  fractographic  data  for  evaluating  the  IFQ 
model  for  different  conditions  Ce.g.,  material,  stress  level, 
spectrum,  %  load  transfer,  fastener  type/fit,  etc.), 

2.  develop  test  results  and  data  needed  to  demonstrate 
the  durability  analysis  methodology  developed,  and 

3.  provide  a  "fractographic  data  base"  that  could  be 
used  in  the  future  to  determine  the  IFQ  for  clearance-fit 
fastener  holes. 

Volume  VIII  documents  the  raw  test  results  and  the  fracto¬ 
graphic  data  and  provides  useful  information  (e.g.,  specimen 
dimensions,  hole  diameters,  photographs,  etc.)  for  future  studies. 

Since  the  "durability"  of  fastener  holes  is  primarily  con¬ 
cerned  with  relatively  small  cracks  (e.g.,  <0.10"),  the  evaluation 
in  this  section  emphasizes  the  TTCI  aspects  more  than  the  TTF. 

TTCI  and  TTF  results  for  selected  fractographic  data  sets  (7475- 
T7 351  aluminum)  are  compared  for  dog-bone  specimens.  A  comprehen¬ 
sive  evaluation  of  crack  initiation  origins  and  trends  was  made 
during  Phase  II;  observations  are  discussed  in  this  section.  Test 
results  for  the  complex  splice  specimen  for  the  F-16  tear-down 
inspection  are  also  presented.  Overall  conclusions  and  recommen¬ 
dations  for  the  Phase  II  test  program  are  presented  in  Section  6.6. 

Several  variables  included  in  the  test  program  were  material, 
load  spectrxom,  stress  level,  fastener  diameter /type,  percent  bolt 
load  transfer,  and  specimen  geometry  (width  and  thickness).  Many 
of  the  factors  which  can  influence  fatigue  life  have  been  studied 
[e.g.,  43-46]. 
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6.2  DOG- BONE  SPECIMENS 


6.2.1  TTCI  Comparisons 


TTCI  comparisons  are  presented  in  Tables  10  and  11  for 
selected  fighter  and  bomber  fractographic  data  sets,  respec¬ 
tively.  The  fractographic  data  sets  listed  in  Tables  10  and 
11  were  the  key  data  sets  used  in  Phase  II  to  evaluate  the 
IFQ  model  for  countersunk,  clearance-fit  fastener  holes.  The 
minimxim,  average  and  maximum  TTCI  values  are  given  for  each 
data  set,  along  with  the  sample  size  used.  In  some  cases, 
the  number  of  TTCI  samples  for  a  data  set  may  be  less  than  the 
total  nximher  of  specimens  tested, 

The  TTCI  values  for  the  selected  data  sets  were  determined 
by  interpolation  from  the  fractographic  data  given  in  Ref,  6 
With  a  few  exceptions,  only  those  specimens  with  fractographic 
data  encompassing  aQ=0,03''  were  used  to  determine  the  minimum, 
average  and  maximum  TTCI  values.  Some  results  were  used  when 
the  last  fractographic  crack  size  was  near  0.03".  Also,  only 
the  fractographic  results  for  the  largest  fatigue  crack  per 
specimen  were  used. 

For  comparison,  the.  minimum,  average  and  maximum  TTCI  re¬ 
sults  from  Table  10  are  plotted  in  Figs.  23-25  for  the  "fighter" 
data  sets  (F-16  400-hour  spectrum).  In  a  similar  manner,  the 
TTCI  results  from  Table  11  are  plotted  in  Figs.  26-28  for  the 
"bomber"  data  sets  (B-1  spectrum).  In  Figs.  23-28,  the  average 
TTCI  values  are  indicated  by  the  open  circles  and  the  sample 
size  used  is  shown  in  parentheses. 


6.2.2  Fatigue  Life  Comparisons 


Average  fatigue  life  is  plotted  as  a  function  of  the 
maximum  baseline  stress  for  selected  data  sets  in  Fig.  29  and  30 
for  the  F-16  400  hour  spectrum  and  B-1  bomber  spectrum,  respec¬ 
tively.  These  plots  are  based  on  results  from  Ref,  6. 

In  general  the  results  plotted  in  Fig.  29  and  30  appear  to 
be  fairly  consistent.  For  example,  the  average  fatigue  life 
decreases  as  the  stress  level  increases.  Also  the  results  tend 
to  show  that  for  a  given  stress  level  the  fatigue  life  decreases 
as  the  percentage  of  load  transfer  increases.  In  Fig.  29,  for 
example,  the  average  fatigue  lives  for  the  407o,  30%,  15%,  and  07o 
load  transfer  data  sets  are  in  the  order  expected.  Exceptions 
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TABLE  10  Comparison  of  TTCI  Results  for  Key  Fighter  Data  Sets  (7475-T7351  Aluminum) 
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MS90353  =  Countersunk,  blind,  pull-through  rivet 

NAS 1580  =  Countersunk  bolt 

cr  =  Gross  stress 

TTCIs  from  Ref.  6 

Ref.  Fig.  16  for  specimen  types 


TABLE  11  Comparison  of  TTCI  Results  for  Key  Bomber  Data  Sets  (7475-T7351  Aluminum) 
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Fig.  23  TTCI  Comparisons  for  Selected  No-Load  Transfer  Data  Sets 
(7475-T7351  Aluminum;  F-16  400-Hour  Spectrum;  a«  -  0.03") 
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Fig,  24  TTCI  Comparisons  for  Selected  15%  Load  Transfer  Data  Sets 
(7475-T7351  Aluminum;  F-16  400-Hour  Spectrum;  aQ  *=  0.03") 
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(7475-T7351  Aluminum;  F-16  400-Hour  Spectrum;  “  0,03”) 
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Fig.  26  TTCI  Comparisons  for  Selected  No-Load  Transfer  Data  Sets 
(7475-T7351  Aluminum;  B-1  Bomber  Spectrum;  aQ  «  0,03") 
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MAX  BASELINE  STRESS  (ESI) 


7A75-T7351  ALUMINUM 
400  HR  BLOCK  F-16  FIGHTER  SPECTRUM 
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Fig.  29  Fatigue  Life  as  a  Function  of  Stress  Level  for 
F-16  400-Hour  Spectrum 
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MAX  BASELINE  STRESS  (KSI) 


7475-T7351  ALUMINUM 
B-1  BOMBER  SPECTRUM 
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Fig.  30  Fatigue  Life  as  a  Function  of  Stress  Level  for 
B-1  Bomber  Spectrum 
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are  also  shown  in  Fig.  29  at  stress  levels  of  32  and  34  Vs-? 


set  (AFXMR4)  at  the  34  ksi  stress  level.  Further  analyses  a: 
required  in  this  case  to  evaluate  the  effects  of  scatter  and 
percent  load  transfer  on  fatigue  life. 


6.2.3  Test  Results  for  D6ac  Steel 


Seventeen  D6ac  no-load  transfer  specimens  (Fig.  21)  were 
tested  using  the  F-16  400-hour  block  spectriom.  All  holes  were 
1/4- inch  diameter  with  NAS- 6204  straight  shank  bolts  installed. 
Test  results  are  summarized  in  Table  12  and  further  details, 
including  fractographic  results,  are  given  in  Ref.  6.  The 
size  of  the  largest  fatigue  crack  in  a  fastener  hole  and  the 
corresponding  number  of  flight  hours  to  failure  are  shown  for 
each  specimen  in  Table  12.  Fatigue  crack  initiation  sites  are 
also  indicated.  The  fatigue  life  as  a  function  of  stress 
level  is  shown  in  Fig.  31. 


6.2.4  Discussion 


Test  results  for  the  dog-bone  specimens  are  discussed  in 
this  section  in  terms  of  the  specimens  and  test  variables.  This 
discussion  is  limited  to  general  observations  and  understandings 
reached  during  Phase  II.  Further  evaluation  of  the  test  results 
[6]  is  recommended. 

The  "load- transfer"  specimens  (see  Figs.  18-20)  were  de¬ 
signed  for  157o,  30%,  and  40%  load  transfer  through  the  fasteners. 
These  specimens  were  designed  by  equalizing  the  axial  deformation 
for  the  upper  and  lower  dog  bones  between  the  lug  pin  and  fastener. 
Pins  and  fasteners  were  also  assumed  to  be  rigid  and  to  "perfectly 
fit"  their  mating  holes.  Due  to  normal  hole  size  variations  and 
clearance  fits,  the  actual  percentage  of  load  transfer  obtained 
varied  from  specimen- to- specimen  and  depended  on  the  magnitude 
of  the  applied  load.  A  typical  free-body  diagram  is  shown  in 
Fig  22. 
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Table  12  Test  Results  for  D6ac  Steel  Specimens 


a 

o 


•H  0) 

M  -H  M 

w  ^ 

a 

CJ 


(U 

u 

o 

PQ 


u  u 

0)  Q) 


c 

o 


0) 

a) 

a> 

0) 

<u 

M 

Pi 

u 

u 

O 

o 

o 

Q 

o 

pq 

pq 

pq 

m 

pq 

o 


M 

Pi  Pi 

Q) 

0)  0) 

C  0) 

p:  c  cu  (u  0) 

P  ^ 

M  Pi  Pi  Pi  P( 

o  o 

O  O  O  O  0 

u  pq 

U  a  pq  pq  pq 

CJ  /-N 

CO  c 

O 

Pi  M 

m 

U 

o 

« 

rH  (U 

c 

CO  N 
a  •H 

V 

•H  CO 

flH 

CO  vo  00  r>*  o 

00  CO  o  lO  o 

'd'  00  o  CO  ^ 

CM  CM  rH  CO  CNJ  CM 


00  00  CM  ^  CM 

lio  vo  ^  00  m 

UO  O  00  1-H 

CM  ^  CM  1-H  CM 


m  ^  CM  CO  vo 

CO  CO  *— »  sO  LO 

0^  vO  o  Mt 

^  CM  CM 


o  o  o  o  o  o 


O  O  O  O  O 


O  O  O  O  O 


no 

0) 


0) 


a  Q) 
B  U 
O  P 

U  tH 
•H 
W  CO 
U  Pm 


3 

o 


rC  <U 

00  pq 

•H 


O 

O 

O 

\D 

A 


vD  vO  UO  C30  0^  00 

O  "d-  CO  ^  m 

CO  o  CM 

^  ^  0% 

00  VO  CM  00  O  o 


a\ 

VO 

C7^ 

uo 

vO 

00 

\o 

CM 

Ov 

CO 

vO 

o 

CO 

CO 

o 

o 

CJs 

m 

vO 

<r 

00 

oo 

00 

oo 

00 

ON 

vO 

0s 

0s 

m 

\0 

00 

<r 

VO 

rH 

CM 

CM 

CO 

0) 

PI 

■H 

•H 

rM 

CO 

Q) 

CO 

CO 

pq 

CO 

CO 

Q) 

» 

u 

X 

•M 

cO 

CO 

0 


LO 

00 


o 

o 


uo 

CM 


OJ 

Pi 

B 

OJ 

•H 

rO 

U 

H 

0) 

P) 

P. 

7^ 

CO 

00  OV  O  ^  CM  CO 
^0  ^0  1—1 


Lo  'O  QO  ov 

•— <  1— I  f— I  1— (  1— j 

r>*  r>. 


O  »H  CM  CO 
O  CM  CM  CM  CM 


I  ■s© 


PLh 

CO 


Ph 

♦J 

CO 


<r 

w 


CO 


CO 

«  0  0  0 


88 


Ref.  21  pp.  B“125  -  B-133  for  fractography 
Gross  section;  F-16  400  hour  spectrum 
Ref.  Fig.  34 


MAX  BASELINE  STRESS  (KSI) 


Fig,  31  Fatigue  Life  as  a  Function  of  stress  Level  for 
D6ac  Steel  Specimens  (P-16  400-Hour  Spectrum) 
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s:z.:“  :?!3S  sr;ri‘srs.i.;  rs  ;r% 

specimens  hed  to  be  "bulkier"  then  those  foi"  ?he  m  lilt  lllZftl 

result,  for  the  same  fastener  fits."  the  30%  and 

157  load  specimens  required  a  larger  applied  load  than  the 

5/c  load  transfer  specimens  before  the  bolts  could  start  transferring 

on^the  load^  the  magnitude  of  the  specimen  loads  varied  depending 
on  the  load  spectrum  and  stress  level,  the  actual  percentaL  of  load 
transferred  through  the  bolts  also  varied.  ^  ^  ^  ^  ^ 

for  transfer  specimens  seemed  to  perform  very  well 

d"a  the  30%  and  40%  load  trLsfer 

specimens  did  not  reach  the  design  percent  load  transfer  consiL 

tran<?f*  d^^  reverse,  double  dog-bone  type  specimen  for  bolt  load 
win  ^  not  assure  the  "design"  percent  bolt  load  traLLr 

will  be  obtained  for  each  load  in  the  spectrum.  To  assure  the 

to  tr^nsf^^^^^^d^d^^  transfer,  the  specimen  should  be  designed 

r:ialJ;;rL^for:ftL^:?^^^  -^^^ectly  by 

also  varies  depending^on  the^fastener^tv^^/f aircraft  structure 
load  levels.  In  order  to  quantitati val applied 
load  transfer  on  the  TTCI  and  TTF  t-bo  ^  assess  the  effects  of 
must  be  controlled.  F^thjf  tests  and^'a^a?"" 

develop  a  better  understanding  of  the  ef fects^of ^bolt^l^^d^'^ 
on  the  TTCI  and  TTF.  exrects  of  bolt  load  transfer 

the  test  program.  ^ThirincluderthrMS-9035rblinr  '"^11 

fit  bolt  because  the  sLnk  i  ®  dearance- 

a  typical  "bucked"  rlve^  it  Slf ''f  like 

and  hole,  for  the  nomlnai  1/trn™  d^L^t^rMriostrr'' 
the  diameter  can  ranse  from  n  9?s  •  f  ^  rivet, 

diameter  according  to  General  n  0.260-inch  while  the  hole 

inch  to  0.263-S  TSrallowrri"?  ^-198  is  0.260- 

0.005-inch  clearance  For  Jha  tolerance  ranging  from  0  to 

bolt,  the  fastener  diamarfr^^  nominal  1/4- inch  diameter  NAS-1580 

the  NAS- 1580  bolt°''^n\lumfnum^all°^^l^''/^''  obtained  with 

decrease  in  fa^igie  Ufe  wirS  ?^  transfer  specimens,  a 

n  xatigue  life  with  increase  in  clearance  fit  has  been 
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observed  in  clearance  fit  holes  [47].  In  general,  specimens  with 
the  larger  diameter  MS-90353  rivet  had  slightly  longer  fatigue 
lives-  than  the  smaller  diameter  rivet.  Clamp-up  normally  increases 

an  accompanying  increase  in  fatigue  life. 
Urzi  [44]  has  studied  the  effects  of  clamp-up  on  fatigue  life  and 
concluded  that  the  fatigue  life  decreases  as  the  stacking-thickness 
to  fastener  diameter  ratio  increases, 

1-h.  1^6]  that  as  the  countersink  depth  approaches 

the  thickness  of  the  material,  the  stress  concentration  factor,  K<- , 
increases  For  countersink  depth/sheet  thickness  ratios  greater 
than  0.5,  Kt  increases  from  a  range  of  2-4  to  a  range  of  8-10.  For 
our  no  load  transfer  specimens,  this  ratio  was  -0.29  for  the  0  375- 
inch  coupons,  and  -0.44  for  the  0.250-inch  thickness  specimens; 
therefore  little  increase  in  stress  concentration  factor  would  be 
expected  Countersink  depths  were  approximately  the  same  for  each 
type  of  1/4- inch  diameter  fastener  (0.105-inch  with  the  MS90353 
rivet  and  O.lll-inchwith  the  NAS-1580  bolt).  In  the  load- transfer 
coupon,  the  countersink  depth/sheet  thickness  ratio  was  -0  58  for 
the  1/4- inch  diameter  fasteners,  therefore,  a  slight  increase  in 
Kt  would  be  expected  in  this  area. 

For  the  same  maximum  baseline  stresc;  i-i, 

was  longer  for  specimens  tested  using  thrB-rb™be^  specir™ 

sneet^°’^  tested  using  the  F-16  400 -hour  block  fights 

spectra.  For  example,  at  a  maximum  baseline  stress  of  14 
ksi,  the  average  fatigue  life  for  the  fighter  Spec^ru^  te^s 

b2i  ^  tested  und^  thf 

omber  spectrum  survived  approximately  3  lives. 

random^Jpectiu^'’®''^h?2®,-''^22  500-hour 

F-16  durability 'test  article!  anris^rii^ir!"^ 

fe!  ^7475^al'^°''’'  spectrum.  Both  07,  and  iS  l^I^^rLs- 

fa:;eL'd'w!trtrMS?9%%^?rs%::rri^^^ 

Specimens  with  two  different  widths  a-nA 
fatigue  tested.  No  conclusionrwLrreaShef on 
specimen  thickness  on  TTCI  and  TTF.  However  th^l  5"^wlde 
specimens  generally  had  slightly  longer  fatigue  livef  than  the 
2  wide  specimens.  The  effects  of  scatter  nee^d  to  be  investigated. 

7475-Tn5°l1nf 202l?8”“™ 

generated  [48]  I^o^dS^L'^bl^''''^^  growth  data  have  been 

for  the  "F^tener  LjfqLutr  program  ft 

yuctj.xcy  program  [41]^  more  testing  was  done 
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using  7475-T7351  material  than  2024-T851  aluminum 
spectrum  fatigue  data  on  2024- T851  has  been  obtained 
other  programs  [47,49]. 


Considerable 
on  two 


Previous  tests  with  an  F-16  spectrum  on  2124-T851  showed 
fatigue  life  than  similar  testing  with  7475-T7651  material 
I  J .  However,  it  has  been  shown  that  overaging  the  7475-T651 
alloy  to  the  T7351  condition  improves  crack  growth  resistance 
under  spectrum  loading  considerably  at  the  expense  of  reduced 
yield  strength  [51]. 


6 . 3  COMPLEX  SPLICE 


testerusine^a“B  specimens  (Fig.  15)  were  fatigue 

testea  using  a  B-1  bomber  load  spectrum.  Specimens  were  tested 

without  intemediate  lateral  support  between  the  lug  ends  using 

servo-controlled,  hydraulically-actuated  load  frame!  11,1^ 

ame  specimen  design,  for  a  different  aluminum  alloy  (2124-T851) 

lateral  support  [42  ’ 

Results  of  the  strain  survey  from  Ref.  42  showed  a  negligibir 
effect  of  the  lateral  support  on  the  resulting  strains® 

Table'^l’T^'^F  splice  test  results  from  Ref.  6  are  summarized  in 

faaif  results  for  fastener  hole  fatigue  cracks 

fastener  hole  dimensions,  final  crack  sizes  for  the  outer  row  ’ 
astener  holes  and  other  testing  details  are  given  in  Ref.  6. 

fatigue  tested  to  failure  except  one 
(CBSi;  22  ksi;  ref.  Table  13).  All  failures  occured  in  th^ 

731^’'(TaMe^l3?^®\-"b''f®-r^  except  for  specimen  number 

(  13)  ,  which  failed  in  the  inner  row  of  holes  Frettine 

appeared  to  be  the  major  cause  of  failure  in  specimen  731.  ^ 

Following  the  fatigue  test,  all  fastener  holes  in  the  outer 

and  bv  th  inspected  for  fatigue  cracks  visually 

and  by  the  eddy  current  technique.  The  size  of  the  largest 

fatipe  crack  in  each  outer  row  fastener  hole  (10  holes  per 

f  [  6  ]  .  A  fractographic  evaluation  was 

largest  fatigue  crack  in  any  one  of  ten 

holes  per  specimen. 

in  Table  13  the  size  of  the  largest  fatigue  crack  per  specimen 

T  ^^^"^^^^^hole  and  the  corresponding  numLr  of 
shown.  The  location  of  the  fastener  hole  with  the 

inrSr  fatigue  crack  initiation  sites  are  also 

oted.  Crack  initiation  sites  were  difficult  to  determine  in 
those  cases  where  cracks  from  adjoining  holes  overlapped.  Also 
locations  of  the  fatigue  crack  origins  varied,  ranging  from  the’ 
countersink-bore  intersection  area  to  the  faying  surface. 


92 


table  13  Complex  Splice  Test  Results 


Q> 

BASEL INI 

!  NO,  FLT3, 

BAT/ 

SPECIMEN 

1  STRESS 

COMPLETED 

SET 

NO. 

© (ksl) 

©  1 

CBBj 

72^ 

22 

3840 

CBSL 

725 

25 

2273 

72  J 

1319 

728 

1769 

729 

1979 

730 

2079 

7J1 

'  2304 

732 

2159 

733 

1954  1 

734 

2589 

1 

' 

735 

1 

1539 

1 

\  , 

71& 

f 

1839 

CBSH 

724 

30 

1079 

737 

759 

738 

769 

739 

899 

740 

3162 

741 

827 

742 

779 

743 

1059 

744  1 

939 

1 

745 

729 

746 

J 

769 

! 

747 

f 

979 

final  crack 

SIZE 

(^>04) 


0, 0276 


0.4142 
0,2711 
0.  2B75 
0.3087 
0. 2698 
0. 1690 
0,3434 
0.3600 
0.3253 
0.3166 
0.2919 


0.2070 

0.2283 

D.2S49 

0.2S72 

0.2955 

0.3010 

0.3183 

0, 3037 

0,2613 

0.2854 

0,2724 

0.2570 


© 

HOLE  NO. 
(OUTER  ROW) 


2 

2 

4 

4 

A 

3 

A 

2 

A 

2 

2 


3 
2 

4 

3 
A 
2 
2 

4 

A 

2 

4 

4 


fatigue  crack  0 

ORIGIN  LOCATION 


(C.S.-B)  Intersection 

(C.S,-B)  " 

Countersink  Area 

Corner 

Bore 

Bore 

Bore 


Comer 

(C.S.-B)  Intersection 


Countersink  Area 
Bore 

(C.S.-B)  Intersection 

Multi:  Bore 

Comer 

(C.S.-B)  Location 
(C.S.-B) 

(C.S.-B)  Location 
Comer,  (C.S.-B) 

Bore 

Bore 


© 

© 

© 

© 

© 

© 


Ref.  6  ,  pp  C-175-C-186  for  fractography 
Gross  section  stress  at  outer  row  of  fasteners 

1  Bomber  Spectrum;  1280  fllghts»13500  fH^hr  h  i 
Largest  fatigue  crack  in  onto,-  /  ^  hours=l  service  life 

iTz:t  ™  II  zrizzr^- 
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Crack  growth  rates  were  much  faster  at  the  stress 


stress  levels  tested 
transfer  coupon  data. 


fer  specimens , 


6.4  CRACK  INITIATION  ORIGINS  AND  TRENDS 


Fatigue  crack  origins  and  trends  were  discussed  in  this 
section  for  the  no-load  (0%),  15%,  30%.  and  40%  load  transfer 
specimens  and  for  the  complex  splice  specimens.  These  obser¬ 
vations  are  based  on  the  comprehensive  fractographic  evalua¬ 
tions  performed  during  Phase  II  [  6] .  All  fracture  surfaces 
considered  were  examined  to  determine,  if  possible,  the  origin 
of  the  dominant  fatigue  crack.  In  many  cases,  fatigue  crack 
origins  were  difficult  to  determine  due  to  the  coalescence  of 
multiple  fatigue  cracks  (Fig.  33). 

A  coding  system,  shown  in  Fig.  34,  was  developed  to  iden- 
tify  fatigue  crack  origins  for  straight-bore  and  countersunk 
fastener  holes.  This  system  applies  to  test  specimens  with 
and  without  a  bolt  load  transfer.  For  example,  in  Fig  34  the 
coding  for  the  no-load  transfer  specimens  applies  to  only 
one  section  and  the  one  with  the  applicable  type  of  fastener 
hole  (i.e.,  countersunk  or  straight-bore). 


6.4.1  No-Load  Transfer  Specimens 


Most  of  the  observed  fatigue  crack  origins  were  randomly 
distributed  in  the  bore  of  the  hole.  In  general,  origin 
locations  in  specimens  tested  with  the  fighter  spectrum  were 
consistent  with  those  tested  using  the  B-1  bomber  spectriam. 
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mx  BASELINE  STRESS  (KSI) 


COMPLEX  SPLICE  SPECIMEN 
B-1  BOMBER  SPECTRUM 


Fig.  32  Fatigue  Life  as  a  Function  of  Stress  Level 
Spectr'l)^''  Specimens  (B-1  Bomber 
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HfrnwF 


FASTENER  HOLE 


Fig.  33  Coalescence  of  Multiple  Cracks  to  Form  a 
Single  Large  Crack 
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CS  -  Countersink  Region 

CS  B  -  Countersink-Bore  Intersection 

B  -  Bore  of  the  Hole 

C  -  Corner 

F  -  Faying  Surface 


Fig.  34  Fatigue  Crack  Origin  Location  Designations 
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Fatigue  cracking  was  distributed  between  the  countersink  and 
straight  bore  portions  of  the  hole.  Typical  crack  origins 
are  shown  in  Fig.  35.  Similar  trends  were  observed  in  Ref. 

41  for  straight-bore  holes. 

Early  fatigue  failures  appeared  to  be  due  to:  (1)  multi¬ 
initiation  sites  coalescing  to  cause  rapid  crack  propagation, 
(2)  large  initiation  sites  were  crack  growth  is  unstable  until 
an  a/2c  ratio  of  0.5  is  obtained,  and  (3)  large  fatigue  cracks 
on  opposite  sides  of  a  fastener  hole. 


6.4.2  157o  Load  Transfer  Specimens 

Several  conclusions  can  be  made  about  crack  origins  and 
failure  modes  in  the  154  load  transfer  specimens. 

1.  Most  fatigue  crack  initiation  sites  were  located  at 
or  near  the  fastener  hole-faying  surface  intersection.  Typical 
crack  origins  are  shown  in  Fig.  36. 

2*  Failures  frequently  occurred  in  the  hole  section  con¬ 
taining  the  countersink.  The  countersink  depth/section  thick¬ 
ness  ratio  was  greater  than  0.5.  Therefore,  a  slight  increase 
in  stress  concentration  would  be  expected  in  this  area, 

3.  Consistent  with  results  obtained  from  no  load  transfer 
specimens,  load  transfer  specimens  with  multiple  initiation 
sites  also  had  shorter  fatigue  lives.  In  general,  the  multi¬ 
origins  caused  a  broader  crack  front  at  an  earlier  stage  of 
propagation.  An  example  of  a  corner  crack  converging  with  a 
crack  in  the  bore  of  the  hole  is  shown  in  Fig.  37.  Again,  as 
was  found  in  no-load  transfer  specimens,  faster  crack  growth 
rate  and  earlier  fatigue  failure  was  observed  in  holes  with 
large  cracks  on  opposite  sides  of  the  hole. 

4.  Similar  crack  origins  and  trends  were  observed  for 
the  specimens  tested  in  Ref.  41. 

5.  Since  a  surface  treatment  was  not  applied  to  the 
faying  surfaces,  some  fatigue  damage  due  to  fretting  occurred. 
Fretting  is  normally  a  more  serious  problem  at  low  fatigue 
stresses,  since  at  higher  stresses  cracks  initiate  rapidly 
even  without  fretting.  No  significant  difference  was  observed 
in  the  failure  modes  of  the  15%  load  transfer  specimens  at 
two  different  stress  levels  (i.e.,  32  ksi  and  40.8  ksi) . 
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F-16  400  HOUR  SPECTRUM 
NO  LOAD  TRANSFER 


Fig.  35 


Typical  Fatigue  Crack  Origins  in  No-Load 
Transfer  Specimen 
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232TA 


232HB 


Fig.  36  Typical  Fracture  Surface  for  15%  Load  Transfer  Specimen 
(ABXMR4(A) ; (r=  38ksi;  B-1  Bomber  Spectrum) 


Fig,  37  Convergence  of  Crack  Fronts  Causing  Rapid  Crack  Growth 
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6.4.3  307o  and  407o  Load  Transfer  Specimens 


Fatigue  crack  initiation  sites  varied  in  these  specimens. 
Many  cracks  originated  at  or  near  the  fastener  hole-faying 
surface  intersection.  A  large  number  were  also  found  at  or 
near  the  countersink-bore  intersection.  Several  early  failures 
were  due  to  broad  origins  around  the  countersink-bore  inter¬ 
sections  . 


In  a  few  of  the  407o  load  transfer  specimens  (AFZLR4,  AFZMR4, 
ABZLR4,  and  ABZMR4) ,  cracks  initiated  on  the  faying  surface  rather 
than  in  the  fastener  hole. 


6.4.4  Complex  Splice 

Fatigue  crack  initiation  sites  were  difficult  to  determine 
when  cracks  from  adjacent  holes  overlapped.  Also  fatigue  crack 
origins  varied,  ranging  from  the  countersink-bore  intersection 
area  to  the  faying  surface  (see  Table  13) . 

Most  of  the  largest  fatigue  cracks  in  the  outer  row  of 
fastener  holes  occurred  in  the  three  holes  near  the  center 
(i.e. ,  hole  number  2,  3,  4,  in  Fig.  38).  In  20  out  of  24 
specimens , the  largest  fatigue  crack  in  the  outer  row  of  fastener 
holes  occurred  in  either  hole  number  2  or  4  (ref.  Fig.  38  and 
Table  13). 

Considerable  fretting  was  observed  in  most  specimens.  In 
two  cases ,  surface  cracks  were  initiated  on  the  faying  surface 
away  from  the  bore  of  the  hole. 


6.5  F-16  LOWER  WING  SKIN  TEAR-DOWN  INSPECTION 


The  full-scale  F-16  wing  box  structure  in' Fig.  39  was 
fatigue  tested  under  the  F-16  durability  certification  program 
[52,53  ].  Both  the  right  and  left  hand  wing  boxes  were  fatigue 
tested  the  same  way  (i.e.,  lab  air,  same  load  spectrum,  stress 
levels,  etc.)  using  the  F-16  500-hour  spectrum. 
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Fig,  38  Identification  of  Fastener  Holes  in  Outer  Row 
of  Complex  Splice  Specimen 


102 


Upper  Attach  Fltt 


16W101  Up 
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16V140  N6.  10 
16V139  Ho.  9  Spai 
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16V134  No 
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16W136  No.  a  Spar 


16W137  No.  7  Spar 


16V13&  Ko.  6  Spar 

! 

16W131  No.  1  Sp«r. 
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I 

16V133  No.  3  Spa 

1 
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Lower  Attach  Fitting 
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16V263,  Sta  325.8 
16V264.  Sta  341.8 
16W265.  Sta  357.8 


Fig.  39  F-16  Basic  Wing  Box  Structure 
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The  F-16  durability  test  article  (wings)  was  fatigue  tested 
to  two  service  lives  Ci.e.,  16000  flight  hours)  using  the  F-16 
500  hour  spectrum.  After  successfully  completing  the  F-16 
durability  certification  requirements,  the  durability  test  art¬ 
icle  (wings)  was  dedicated  to  the  "Durability  Methods  Develop¬ 
ment"  (DMD)  program. 

A  tear-down  inspection  of  the  F-16  durability  test  article 
(wings)  was  performed  under  Phase  II  of  the  DMD  program.  All 
fastener  holes  and  cutouts  in  both  wing  lower  skins  (7475-T7351 
aluminum)  were  inspected  (over  3000  holes)  using  eddy  current 
techniques.  Also,  approximately  forty  fastener  holes  per  spar 
(20  upper  and  20  lower  at  the  inboard  end)  were  inspected. 


All  fastener  holes  in  the  lower  wing  skins  with  an 
indicated  crack  were  broken  open  and  then  evaluated  fracto- 
graphically.  Fastener  holes  with  an  indicated  crack  after 
16000  hours  of  fatigue  testing  (F-16  500  hour  spectrum)  are 
identified  in  Figs.  40  and  41  for  the  right  hand  and  left 
hand  lower  wing  skin,  respectively.  Fractographic  results 
are  summarized  in  Tables  14  through  17.  Results  of  the 
F-16  wing  tear-down  inspection  were  used  to  demonstrate 
the  durability  analysis  methodology  [21  ,  34,  36] .  The 
tear-down  inspection  is  also  documented  in  Volume  VIII  and 
elsewhere  [11,15,18]. 


6.6  CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  and  recommendations  are  based 
on  the  Phase  II  test  program  and  results  [  6 ] • 

1.  A  wealth  of  fatigue  crack  initiation  and  crack  pro¬ 
pagation  data  for  clearance-fit  fastener  holes  has  been  dev¬ 
eloped  and  documented  under  this  program  [‘6  ] .  These  results 
can  be  used  to  determine  the  IFQ  for  clearance-fit  fastener 
holes  for  future  durability  analyses. 

2.  The  "ton"  of  data  in  Volume  VIII  needs  to  be  further 
evaluated  to  better  understand  the  effects  and  significance  of 
various  test  and  specimen  variables  on  the  TTCI  and  TTF  (e.g., 
7o  bolt  load  transfer,  fastener  type/size/fit ,  specimen  size, 
clamp-up,  etc.).  A  statistical  analysis  of  the  data  is  needed 
to  evaluate  the  effects  of  scatter  and  to  determine  trends . 
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Fig,  41  Identification  of  Fastener  Holes  in  the  F-16  Durability 
Component  Lower  Left  Hand  Wing  Skin  1th  Observed  Cracks 
CAfter  16000  Hrs,  of  Testing  with  F**16  500-Hour  Spectrina.) 
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Table  15  for  F-16  I,o,«r  King  Skin  (Right  Hand; 

Durability  Test  Article)  -  Stress  Zone  II 
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P-16  500-Hour  Block  Spectrum  ***Perimeter  dimensions  along  direction 

Vaoue  growth/crack  measurements 
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Table  17  Fractographic  Results  for  F-16  Lower  Wing  Skin  (Left  Hand; 
Durability  Test  Article) 
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*** 

-16  500 -Hour  Block  Spectrum  Perimeter  dimension  along  direction 

of  crack  growth/ crack  measurements 


3.  Other  fractographic  data  sources  for  fastener  holes 
[e.g.,  47,  49,  54,  55]  should  also  be  evaluated  and  compared 
with  the  results  in  Volume  VIII  [6] . 

4.  Variable  percentages  of  bolt  load  transfer  were 
obtained  using  the  double,  reverse  dog-bone  type  specimens 
(see  Figs.  18-20).  The  actual  amount  of  bolt  load  transfer 
for  each  specimen  depends  on  factors  such  as  the  fastener 
type  and  fit,  load  spectrum,  stress  level,  etc.  Since  the 
actual  percentage  of  bolt  load  transfer  varies  for  clearance- 
fit  fastener  holes  in  aircraft  components,  test  specimens  with 
variable  percentages  of  bolt  load  transfer  may  be  satisfactory 
for  generating  the  IFQ  for  fastener  holes.  The  effects  of  bolt 
load  transfer  on  the  TTCI  and  IFQ  need  to  be  investigated  fur¬ 
ther,  using  specimens  with  a  controlled  load  applied  directly 
to  the  bolt. 

5.  Specimens  with  countersunk  fastener  holes  typically 
exhibited  more  "scatter"  in  the  fractographic  results  than 
those  for  straight  bore  holes  [e.g.,  41].  The  increased 
scatter  for  the  countersunk  fastener  holes  is  probably  due 
to  differences  in  :  (1)  fit  for  countersunk  and  protruding 
head  fasteners  and  (2)  crack  initiation  origins. 

6.  The  following  crack  initiation  origins  and  trends 
were  observed: 

(1)  No-load  transfer  Specimens  -  Most  origins  were 
in  the  bore  of  the  hole.  Multiple  fatigue 
cracks  converging  into  a  single  crack  front 
were  also  observed. 

(2)  15%  Load  Transfer  Specimens  -  Fatigue  cracks 
generally  originated  at  the  faying  surfaces  - 
bore  intersection.  Signs  of  fretting  were 
also  observed  in  several  specimens. 

(3)  30%  and  40%  Load  Transfer  Specimens  -  Origins 
of  failure  varied  considerably.  Many  origins 
were  at  or  near  the  fastener  hole- faying  surface 
intersection  and  at  or  near  the  countersink-bore 
intersection.  Signs  of  fretting  were  observed 
but  the  problem  appeared  to  be  no  worse  than  that 
for  the  15%  load  transfer  specimens. 
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(4)  Complex  Splice  Specimens  -  Fatigue  crack  origins 
varied  considerably,  ranging  from  the  countersink- 
bore  intersection  area  to  the  faying  surface  (see 
Table  13) .  Initiation  sites  were  particularly 
difficult  to  determine  when  cracks  from  adjacent 
holes  overlapped.  Most  of  the  largest  fatigue 
cracks  occurred  in  the  outer  row  of  fastener 
holes  near  the  centerline.  Considerable  evi¬ 
dence  of  fretting  was  observed. 

7.  Early  fatigue  failures  were  observed  when  cracks  of 
approximately  equal  size  occurred  on  both  sides  of  the  hole. 
Fractographic  results  were  obtained  for  the  largest  fatigue 
crack  per  specimen  and  for  selected  secondary  cracks.  A 
small  fatigue  crack  (e.g.S  0.10")  for  one  side  of  a  fastener 
hole  does  not  appear  to  have  a  significant  effect  on  the. 
crack  initiation  for  a  crack  on  the  opposite  side  of  the  hole 
and  vice  versa.  The  effects  of  cracks  on  both  sides  of  the 

hole  and  crack  geometry  on  TTCI  and  TTF  need  to  be  investigated 
further. 

8.  The  wealth  of  available  fractographic  data  for  fastener 
holes  [e .g. ,  6 , 24, 41 , 47 , 49 , 54, 55]  should  be  used  to  determine 
suitable  IFQ  model  parameters  for  durability  analysis.  IFQ  Model 
parameters  should  be  determined  for  pooled  fractographic  data 
sets  and  the  accuracy  of  the  extent  of  damage  predictions  evalu¬ 
ated  to  increase  confidence  in  the  durability  analysis  methodology. 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

7.1  CONCLUSIONS 


The  durability  analysis  methodology  developed  in  Phase  I 
of  this  program  has  been  evaluated  and  refined  in  Phase  II. 
This  methodology  has  been  demonstrated  for  coupon  data, 
full-scale  fighter  aircraft  structure  {i.e.,  F-16  lower  wing 
skin)  and  a  complex  splice  subjected  to  a  B-1  bomber 
spectrxim.  The  analytical  tools  developed  under  this  program 
can  be  used  to  analytically  assure  design  compliance  with 
the  U.S.  Air  Force's  durability  design  requirements. 
Although  the  methodology  has  been  demonstrated  only  for 
fastener  hole  cracking,  the  basic  durability  strategy  and 
concepts  apply  to  cracking  in  any  structural  detail  (e.g., 
cutouts,  lugs,  fillets,  etc.).  The  methodology  has  been 
demonstrated  for  quantifying  the  extent  of  damage  in  grouped 
fastener  holes  due  to  relatively  small  sub-critical  cracks 
(e.g.  ^  0.10").  Further  work  is  required  to  demonstrate  the 
methodology  for  larger  crack  sizes. 

A  large  f ractographic  data  base  has  been  developed  for 
cracking  in  fastener  holes.  Much  has  been  learned  about  IFQ 
characterization  and  the  associated  parameters  involved.  An 
effective  IFQ  model  has  been  developed,  evaluated  and 
refined  for  quantifying  the  IFQ  for  structural  details 
susceptible  to  fatigue  cracking.  The  IFQ  model  is  a  very 
promising  tool  for  quantifying  the  initial  manufactured 
state  of  structural  details  in  terms  of  equivalent  initial 
flaw  sizes. 

The  following  -  conclusions  and  observations  are  based  on 
the  results  and  investigations  of  this  program: 

1.  Several  useful  applications  of  the  durability 
analysis  tools  developed  are:  (1)  evaluate  durability  design 
trade-offs  in  terms  of  structural  design  variables,  (2) 
evaluate  structural  maintenance  requirements  before  or  after 
aircraft  is  committed  to  service,  and  (3)  provide  aircraft 
user  options  affecting  life-cycle-costs,  structural 
maintenance  requirements,  and  operational  readiness. 

2.  The  IFQ  model  is  simply  a  "mathematical  tool"  for 
quantifying  the  IFQ  of  structural  details.  Therefore,  the 
resulting  EIFS's  must  be  considered  in  the  context  of  the 
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IFQ  model  and  the  f ractographic  results  used  to  calibrate 
the  model  parameters.  EIFS's  should  be  considered  as 
hypothetical  cracks  used  to  make  crack  exceedance 
predictions  rather  than  actual  initial  cracks  per  se. 

3.  Back  extrapolations  of  f ractographic  data  must  be 
done  consistently  to  put  the  EIFS's  on  a  common  baseline  for 
different^  data  sets.  _  Inconsistent  EIFS  results  will  be 
obtained  if  the  EIFS  distribution  is  determined  by  back 
extrapolating  the  fractographic  '  results  for  individual 
specimens  and  then  fitting  a  statistical  distribution  to  the 
EIFS  results  for  different  data  sets.  Two  problems  result 
if  this  approach  is  used:  (1)  the  EIFS's  are  not  on  a  common 
baseline  for  different  data  sets,  and  (2)  the  resulting  EIFS 
distribution  is  not  statistically  compatible  with  the  TTCI 
distribution  and  the  fatigue  wear  out  process.  The 
result-ing  EIFS  distribution  should  be  statistically 
compatible  vith  the  TTCI  distribution.  The  IFQ  model 
developed  satisfies  this  requirement. 

4.  The  IFQ  model  developed  can  be  used  to  define  the 
EIFS  cumulative  distribution  using  suitable  fractographic 
results.  The  par^eters  a  and  Q*^.  provide  the  basis  for 
putting  fractographic  results  on  ^a^  common  baseline  for 
quantifying  the  initial  fatigue  quality,  a-  and  Q*/J.  should 
be  constants  for  different  fractographic  data  siti  (same 
material,  fastener  type/fit,  and  drilling  technique  but 
different  loading  spectra,  stress  level  and  percent  .load 
transfer) . 


5.  For  generic  EIFS,  a  and  qJ^-  should  be  constants. 
Encouraging  results  have  been  obtained  to  justify  the  use  of 
the  same  EIFS  cumulative  distribution  for  making  crack 
exceedance  predictions  for  different  design  conditions. 
Further  research  is  required  to  confirm  the  IFQ 
distributions  for  different  materials,  load  spectra,  stress 
levels,  fastener  types/diameters/fit,  %  load  transfer,  etc. 
A  considerable  amount  of  fractographic  results  exist  which 

need  to  be  evaluated  using  the  IFQ  model  [e . g.  , :6 , 41 , 47 , 49 , 

54^5  5  jj  * 


6.  Theoretically,  the  IFQ  model  can  be  used  to  quantify 
the  EIFS  cumulative  distribution  for  various  structural 
details  so  long  as  fractographic  results  are  available  for 
the  details  to  be  included  in  the  durability  analysis.  The 
IFQ  model  has  been  evaluated  using  fractographic  results  for 
fastener  holes.  Suitable  specimens  and  guidelines  need  to 
be  developed  for  generating  crack  initiation  and  crack 
growth  results  for  details  such  as,  cutouts,  fillets,  lugs, 
etc.  Fractographic  results  should  be  developed  and 
evaluated  for  such  details  so  that  the  durability  analysis 
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methods  described  can  be  efficiently  applied  to  different 
types  of  structural  details  in  typical  aircraft  structures. 

7.  The  accuracy  of  crack  exceedance  predictions,  based 
on  the  same  EIFS  cumulative  distribution,  needs  to  be 
evaluated  for  different  design  conditions.  Also,  IFQ  model 
parameter  sensitivity  studies  need  to  be  performed  to  better 
understand  the  average  parameter  values  and  variances  and 
the  impact  of  these  parameters  on  the  IFQ  for  different 
f ractographic  data  sets. 

8.  The  durability  analysis  method  was  developed  for 
making  crack  exceedance  predictions  for  relatively  small 
crack  sizes  (e.g.,  ^  0.10")  in  structural  details.  The 
largest  crack  in  each  detail  was  assvimed  to  be  statistically 
independent  to  justify  using  the  binomial  distribution  for 
combining  crack  exceedance  predictions  for  structural 
details.  If  the  largest  crack  in  a  given  detail  doesn't 
significantly  affect  the  growth  of  cracks  in  neighboring 
details,  perhaps  the  proposed  durability  analysis  method  can 
be  extended  to  crack  sizes  >  0.10".  The  simplistic  crack 
growth  rate  equation  (Eq.  2)  used  in  the  IFQ  model  may  not 
be  suitable  for  defining  the  EIFS  cumulative  distribution 
for  crack  sizes  >  0.10".  Other  functional  forms  for  crack 
growth  rate  may  be  required  to  justify  the  same  EIFS 
distribution  for  both  small  and  large  crack  sizes.  Eq.  2 
may  be  acceptable  for  making  p(i,T)  predictions  for  larger 
crack  sizes  if  the  "EIFS  master  curve"  is  curve  fit  to  the 
larger  crack  sizes.  This  would  result  in  a  different  EIFS 
cumulative  distribution  for  the  small  and  large  crack  size 
range.  Further  research  is  required  to  extend  the 
probabilistic  fracture  mechanics  approach  developed  to 
larger  crack  sizes. 

9.  Two  different  Fa(0)^*)  equations  (i.e.,  Eqs.  8  and 
13)were  presented  for  representing  the  IFQ.  Either  equation 
is  valid  but  Eq.  13  is  recommended  for  two  reasons:  (1)  it 
assures  all  EIFS's  in  the  IFQ  distribution  will  be  >  0,  and 
(2)'  the  crack  growth  rate  parameter.Q  or  Q^,  can  be  easily 
determined  from  the  f ractographic  results  and  the  resulting 
qJ  values  for  different  data  sets  will  be  directly 
comparable.  If  Eq.  8  is  used,  a  common  b  parameter  (Eq.  2) 
must  be  imposed  for  different  f ractographic  data  sets  to  put 
the  Q±P±  values  on  a  comparable  baseline.  As  long  as  b  > 
1,  all  EIFS's  in  Eq.  8  will  be  ^  0. 

10.  The  EIFS  cumulative  distribution  ^a(O) 
independent  of  the  TTCI  reference  crack  size,  a,,  used  (ref. 
Eqs.  8  and  13).  By  using  consistent  fractographic  data 
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pooling  procedures,  results  for  different  data  sets  can  be 
put  on  a  common  baseline. 

11.  Fractographic  results  for  cracks  in  given 
structural  details  are  required  to  calibrate  the  IFQ  model 
parameters.  The  IFQ  distribution  is  based  on  the  observed 
fractographic  results  without  direct  consideration  of  small 

cracky  kinetics.  This  means  that  the  IFQ,  or  EIFS 
distribution,  can  be  determined  without  addressing  directly 
the  question  of  small  crack  growth  behavior.  The  effects  of 
small  crack  growth,  if  desired, .  could  be  accounted  for  in 
the  EIFS  master  curve  used  to  define  the  IFQ  distribution. 
Due  to  the  generality  of  the  IFQ  model  concepts,  IFQ  can  be 
determined  using  state-of-the-art  understandings  of  the 
crack  growth  behavior. 

12.  The  EIFS.  concept  is  convenient  for  quantifying  IFQ. 
Care  and  judgment  must  be  exercised  when  defining  the  IFQ 
distribution  to  assure  EIFS  results  are  on  a  comparable 
basis  for  oifferent  data  sets.  One  of  the  most  important 
lessons  learned  from  this  program  is  that  EIFS's  must  be 
defined  using  pooled  fractographic  results  to  make  the 
EIFS's  for  one  data  set  comparable  with  those  for  another 
data  set.  If  a  constant  a  can  be  justified  for  the  pooled 
fractographic  data  sets,  this  provides  a  simple  way  to  put 
the  TTCI  results  on  a  compatible  basis.  In  order  to  obtain 
compatible  EIFS  results  for  different  fractographic data 
sets,  the  EIFS  master  curve  must  also  be  defined  in  a 
consistent  manner.  This  is  very  important  because  the 
resulting  IFQ  distribution  depends  on  the  transformation  of 

the  TTCI  distribution  using  the  applicable  EIFS  master 
curve. 


13.  ^  The  EIFS  master  curve  depends  on  the  fractographic 
crack  size  range  used  and  the  acceptable  curve  fit 

master  curve  changes,  the  resulting 
EIFS  s  also  change.  Therefore,  different  EIFS's  can  be 
obtained  for  the  same  set  of  TTCI's.  Use  consistent 
procedures  and  fractography  data  pooling  for  defining  EIFS's 
to  make  results  directly  comparable  —  otherwise  EIFS 
comparisons  will  be  meaningless. 

consistent fractographic  crack  size 
ranges  for  different fractographic  data  sets  when  defining 
the  corresponding  EIFS  master  curves.  For  example,  a  crack 
size  range  from  0.005"  to  0.100"  appears  reasonable,  for 
quantifying  the  IFQ.  Consistent  crack  size  ranges  must  be 
used  to  obtain  compatible  EIFS  master  curves  for  different 
fractographic  data  sets. 
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15.  The  durability  analysis  methodology  has  been 
demonstrated  for  the  F-16  lower  wing  skin  and  for  a  complex 
splice  specimen  subjected  to  a  B-1  bomber  spectrum.  In  both 
cases,  f ractographic  results  for  load  transfer  coupon 
specimens  were  used  directly  to  predict  the  extent  of  damage 
for  the  corresponding  full-scale  structures.  Excellent 
correlations  were  obtained  without  considering  scale-up 
effects  (i.e.,  coupon  structure  to  full-scale  structure). 
However,  such  effects  should  be  investigated  further  to 
determine  their  significance  for  durability  analyses. 

16.  Fractographic  data  are  available  for  countersunk 
fasteners  for  two  different  diameters  (Vol.  VIII).  The 
effects  of  fastener  diameter  on  the  IFQ  distribution  have 
been  studied  but  further  research  is  needed  to  reach  a 
better  understanding. 

17.  IFQ  model  calibration  studies  (see  Appendices  A  and 
B)  were  made  using  fractographic  results  for  straight  bore 
[41]  and  countersunk  fastener  holes  fG/]  .  These 
studies  were  based  on  the  Q[a(t)P  model  for  the  EIFS  master 
curve.  The  following  observations  are  based  on  these 
studies  for  7475-T7351  aluminum: 

o  cc  and  b  appear  to  be  material  constants 

independent  of  load  spectra  and  stress  level. 

o  Typical  ft  values  ranged  from  approximately 

4.4  for  straight  bore  holes  to  2.5  for  counter¬ 
sunk  holes.  More  scatter  in  the  results  were 
observed  for  countersunk  holes  than  for  the 
straight  bore  holes. 

o  Typical  b  values  of  0.88  and  0.61  were 

observed  for  straight  bore  and  countersunk 
holes,  respectively. 

o  Average  values  of  2.8  and  0.62  were 

observed  ifer  straight  bore  and  countersunk 
bore  holes,  respectively  (aQ*0.03"). 

o  appeared  to  be  independent  of  load 

spfectra,  stress  level,  and  percent  load 
transfer. 

Further  studies  are  needed  to  evaluate  the  IFQ  distribution 
based  on  the  Q[a(t)]  model  for  the  EIFS  master  curve.  A 
better ^  .  understanding  of  the  IFQ  model  parameter 
sensitivities  to  design _  variables  is  needed.  Further 
studies  should  be  made  using  available  fractographic  data 
,[e.g.  ,  6,  41,  47,  49,  54,55]  . 
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18.  Variable  shear  load  transfer  through  the  fasteners 
was  obtained  for  the  various  %  load  transfer  fractographic 
data  sets.  This  variability  was  apparently  due  to  the 
following; 

o  Variable  clearance  between  fasteners  and  mating 
holes  (test  section  and  lugs) . 

o  insufficient  axial  deformation  under  applied 
spectrum  loading  to  obtain  full  design  load 
transfer  in  specimen  for  each  load  in  spectrum. 

Real  aircraft  structure  will  also  have  variable  %  load 
transfer  due  to  variable  fits  between  fasteners  and  holes. 
However,  to  properly  assess  the  effects  of  fastener  shear 
load  transfer  on  IFQ,  fatigue  cracking  in  fastener  holes 
should  be  obtained  under  controlled  %  load  transfer 
conditions . 


7 . 2  RECOMMEND AT I ONS 


1.  Further  work  is  needed  to  evaluate  the  accuracy  and 
sensitivity  of  p(i,T)  predictions  based  on  a  given  IFQ 
distribution  ^  and  service  crack  growth  master  curves 
determined  using  a  suitable  analytical  crack  growth  program 
fe.g.,  37  -  40  ].  This  research  is  necessary  to  increase  confidence 
in  the  p(i,T)  predictions  for  a  wide  range  of  "real  world" 
durability  applications. 

2.  The  IFQ  of  fastener  holes  should  be  studied  further 

using  the  IFQ  model  developed.  Existing  f ractographic 
results  from  various  sources  should  be  used  to  evaluate  IFQ 
model  parameters  (i.e.,  a  ,  0  ,  e  ,  Q  and  b  )  and  to 
establish,  on  a  statistical  basis,  suitable  IFQ 

distributions  for  durability  design  applications.  A  wealth 
of  f ractographic  results  are  available  from  this  program 
(Vol.  VIII).  Other  f ractographic  data  sources  are  also 
available  [e.g. ,  24,41,47,49,54,551  . 

3.  The  accuracy  of  crack  exceedance  predictions,  based 
on  the  same  EIFS  cumulative  distribution,  need  to  be 
evaluated  for  different  design  conditions.  Also,  IFQ  model 
parameter  sensitivity  studies  need  to  be  performed  to  better 
understand  the  average  parameter  values  and  variances  and 
the  impact  of  these  parameters  on  the  IFQ  for  different 
f ractographic  data  sets. 
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4.  Fractographic  data  from  various  sources  should  be 
pooled  and  IFQ  model  parameters  determined.  Model 
parameters  should  be  analyzed  statistically  and  a  "generic" 
EIFS  distribution  for  IFQ  should  be  determined,  if  possible, 
for  a  wide  range  of  durability  design  applications. 

5:  Further  research  is  required  to  confirm  the  IFQ 

distributions  for  different  materials,  load  spectra,  stress 
levels,  fastener  types/diameters/fit,  %  load  transfer,  etc. 

A  considerable  amount  of  fractographic  results  exists  which 
need  to  be  evaluated  using  the  IFQ  model. 

6.  The  durability  analysis  methodology  developed  under 
this  program  needs  to  be  clinically  evaluated  using  the 
wealth  of  fractographic  data  available.  Basic  assumptions, 
analysis  procedures,  design  guidelines  and  overall 
implementation  need  to  be  critically  reviewed  and  the 
methodology  appropriately  refined  if  necessary. 

7.  The  effects  of  fretting,  clamp-up,  corrosion,  size 
effect  (scale-up  from  coupon  to  component),  faying  surface 
sealant,  interference-fit  fasteners,  residual  stresses,  etc., 
on  IFQ  need  to  be  investigated  and  guidelines  should  be 
determined  for  durability  design  applications. 

.  V feasibility  of  using  no-load  transfer  specimens 
with  multiple  holes  for  quantifying  the  IFQ  should  be 
evaluated  using  spectrum  and  constant  amplitude  loading. 
This  could  provide  an  economical  way  to  generate  the 
fractographic  results  needed  to  quantify  the  IFO 
distributions. 

9.  The  effect  of  percent  fastener  load  transfer  on  IFQ 
♦should  be  investigated  further.  The  "load  transfer" 

specimens  for  this  program  were  designed  for  15%,  30%  and 
40%  load  transfer  assuming  a  perfect  fit  between  fasteners 
and  holes  as  well  as  lug  pins  and  holes.  Since  the  fit 
varies  from  hole  to  hole,  the  %  load  transfer  also  varies. 
Therefore,  to  quantitatively  evaluate  the  effects  of  load 
transfer  on  IFQ,  specimens  need  to  be  fatigue  tested  in  such 
a  way  to  assure  a  specified  %  shear  load  transfer  through 
the  fastener.  ^ 

10.  The  effect  of  scale-up  (e.g.,  coupon-component),  on 
p(i,T)^  predictions  needs  to  be  investigated  further 
including  an  assessment  of  the  accuracy  of  the  predictions 
for  different  design  conditions. 

11.  The  effect  of  fastener  size  on  IFQ  needs  to  be 
investigated  further. 
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12.  The  durability  design  handbook  [21]  should  be 
expanded  to  include  more  "durability  design  data"  for 
different  materials  and  for  the  key  structural  details 
governing  durability. 

13.  Evaluate  and  verify  the  durability  analysis 
methodology  developed  under  this  program  for  cutouts,  lugs 
and  fillets.  Refine  the  methodology  if  required. 

14.  Develop  suitable  test  specimens  and  guidelines  for 
generating  crack  initiation  and  crack  growth  data  needed  to 
quantify  the  IFQ  for  cutouts,  lugs  and  fillets. 

15.  Generate  f ractographic  results  for  cutouts,  lugs 
and  fillets  for  different  materials  and  design  conditions 
(e.g.,  load  spectra,  stress  level,  environment,  etc.).  Use 
these  data  to  determine  the  corresponding  IFQ  for  cutouts,- 
lugs  and  fillets.  Study  the  resulting  IFQ  model  parameters 
and  develop  understanding  for  design  applications. 

16.  The  durability  analysis  methodology  developed  under 
this  program  should  be  extended  to  larger  crack  sizes  (e.g., 
a  >  0.10")  and  verified  for  key  structural  details  governing 
durability  (i.e.,  fastener  holes,  cutouts,  lugs,  and 
fillets) . 


17.  Further  work  is  needed  to  determine  reasonable 
durability  crack  size  limits  and  guidelines  for  the 
economical  repair  of  cutouts  and  fillets.  Repair  crack  size 
limits  are  needed  for  applicable  structural  details  to 
quantify  the  extent  of  damage  predictions  on  a  rational 
basis . 

18.  Allowable  extent  of  damage  limits  for  airframe 

structures  and  guidelines  need  to  be  developed  so  that 

"economic  life"  can  be  defined  in  meaningful  terms  for 

analytically  demonstrating  design  compliance  with  the  U.S. 
Air  Force's  durability  design  requirements. 

19.  An  efficient,  general-purpose  plotting  capability 
should  be  developed  for  presenting  durability  analysis 
results  in  various  useful  formats.  An  automated  plotting 
capability  has  been  developed  under  this  program  but 
improvements  are  needed  to  obtain  plots  more  efficiently. 

,The  program  should  be  able  to  handle  large  amounts  of 
f ractographic  data.  This  capability  could  be  used  to 

economically  generate  useful  durability  design  data  for  the 

durability  design  handbook  and  to  evaluate  the  accuracy  of 
durability  analysis  predictions. 
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20.  Durability  analysis  workshops  should  be  conducted 
to  transfer  the  technology  developed  under  this  program  to 
various  user  levels.  For  example,  such  workshops  would  be 
useful  to  aerospace  engineers  who  are  directly  concerned 
with  aircraft  structural  durability  analyses  and  compliance 
with  the  U.S.  Air  Force's  durability  design  requirements. 
This  would  provide  an  excellent  forum  for  exchanging  ideas 
and  explaining  the  application  of  the  durability  analysis 
methodology  developed.  Such  exchanges  should  lead  to 
methodology  improvements  and  further  advance  the  durability 
analysis  state-of-the-art. 
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APPENDIX  A 


EVALUATION  OF  IFQ  MODEL  PARAMETERS  FOR 
STRAIGHT- BORE  FASTENER  HOLES 

A.l  INTRODUCTION 

The  IFQ  model  described  in  Section  3.5  is  evaluated  in  this 
Appendix  using  fractographic  results  for  straight-bore  fastener 
holes  [41] .  Similar  results  for  countersunk  fastener  holes  are 
documented  in  Appendix  B  of  this  report.  Model  parameters  were 
determined  for  both  individual  and  pooled  fractographic  data  sets 
using  the  procedures  described  in  section  4.3. 

The  basic  objective  of  this  study  was  to  compare  IFQ  model 
parameters  for  different  fractographic  data  sets  for  consistency 
and  trends  and  to  learn  more  about  the  application  of  the  model. 

A  secondary  objective  was  to  evaluate  possible  relationships  for 
Q  and  Q  and  ^  and  ^  and  The  IFQ  model  parameters  for  this 
study  were  determined  using  preliminary  calibration  procedures 
rather  than  those  recommended  in  the  Durability  Design  Handbook 
[21]  .  However,  this  study  and  others  fs.g., 11-19,  31-36]  provided 
the  "experience”  needed  to  develop  the  final  understanding  for  the 
program  [21] . 

The  fractographic  data  sets  used  for  this  study  are  described 
in  Table  A1  and  details  are  given  in  Ref.  41.  Specimen  details  for 
the  no-load  transfer  and  the  15%  load  transfer  specimen  are  shown 
in  Fig.  Al.  All  specimen  fastener  holes  were  drilled  using  a 
modified  Winslow  spacematic  drill  and  holes  were  not  deburred.  NAS 
6204 (V'  dia.)  protruding  head  bolts  (clearance-fit)  were  installed 
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TABLE  A1  FHQ  DATA  SETS  USED 


DATA 

SET 

7o  LOAD 
TRANSFER 

NO. 

SPECIMENS 

a 

(ksi) 

TEST 

SPECTRUM 

MATERIAL 

WPF 

0 

37 

34 

F-16  400-Hr 

7475- 

•T7351 

XWPF 

15 

33 

34 

Block 

HYWPF 

15 

10 

40.8 

LYWPF 

15 

10 

30.6 

1 

f 

WPB 

0 

32 

34 

B- 

■1 

XWPB 

15 

31 

34 

HYWPB 

15 

10 

40.8 

LYWPB 

15 

10 

30.6 

1 

Notes 

1.  Fractographic  results  from  Ref.  41. 

2.  NAS  6204  (k  dia)  protruding  head  bolt. 
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(b) 


Fig,  Al  Fastener  Hole  Quality  Specimens;  (a)  No-Load 
Transfer  and  (b)  15%  Load  Transfer 
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in  the  fastener  holes  as  shown  in  Fig.  Al. 

This  study  also  reflected  the  following: 
o  e  *  0  and  ao  =  (ref.  Eq.  8) 
o  b  1  form  of  Eq.  2 

o  two  fractographic  crack  size  ranges  used: 

0.005"  =  a(t)  =  0.10"  and  0.005"  =  a(t)  =  0.050". 

IFQ  model  parameter  values  are  summarized  herein.  Further 
details  of  the  study  and  conclusions  are  given  in  subsection  4.8.1. 


A. 2  IFQ  MODEL  PARAMETERS  FOR  INDIVIDUAL  DATA  SETS 

In  Table  A2,  IFQ  model  parameters  are  summarized  for 
individual  fractographic  data  sets.  These  parameters  were 
determined  using  the  fractographic  results  from  Ref.  41  and 
the  calibration  procedures  described  in  section  4.3  Only  the 
fractography  for  the  largest  fatigue  crack  per  specimen  was 
used  to  determine  the  IFQ  model  parameters  (i.e. ,  £^=1;  no  $ 
scaling  used) . 

The  notation  for  the  IFQ  Model  parameters  shown  in  Table  A2 
is  consistent  with  Ref.  21.  For  example, 
oe^  =  Weibull  shape  parameter  for  the  ith  TTCI  data  set. 


★  ic  da(t)  s  i 

b^,Q^  =  crack  growth  constants  in  =  Q^[a(t)] 

for  the  ibh  fractographic  data  set. 

=  Weibull  scale  parameter  for  the  i^h  tTCI  data  set. 
^  Note:  ,  , 

^i  ~  (Ref.  Section  4.4). 
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r 

\ - 

PQ 

o 

X 

•^H 

O' 

1.306 

1.964 

4.218 

1.924 

1.096 

1.793 

1.602 

1.937 

•H 

o? 

ca 

16992 

11352 

4230 

25386 

37092 

16956 

7672 

34428 

■K*H 

0.8399 

0.8613 

0.8277 

0.9765 

0.9211 

0.9224 

0.7166 

1.0574 

•H 
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5.349 

4.157 
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4.375 

3.934 

4.671 
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A.  3  IFQ  MODEL  PARAMETERS  FOR  POOLED  DATA  SETS 

IFQ  Model  parameters  were  determined  using  the  pooled 
fractographic  results  for  different  data  sets  and  data  set 
combinations.  Two  different  fractographic  crack  size  ranges 
were  used, i.e.,  0.005”  =  a(t)  =  0.100”  and  0.005”  =  a(t)  = 

0.050”.  The  computed  IFQ  Model  parameters  for  the  different 
fractographic crack  size  ranges  are  summarized  in  Tables  A3  and 
A4. 

In  Table  A3  and  A4,  a  and  b  were  determined  using  the 
pooled  fractographic  results  for  the  data  sets  indicated.  Given 
”ol” ,  the  corresponding  3^^  value  for  each  data  set  pooled  was 
determined  using  the  procedures  described  in  Section  4.3.  0*  for 

each  data  set  was  determined  using  the  applicable  ”b”  value  for 
the  pooled  data  sets. 

The  product  is  tabluated  in  Table  A3  and  A4  for  each 

data  set  and  data  pooling  combination.  Also,  the  average  Q|B£^ 
and  the  standard  deviation  and  coefficient  of  variation  for 
Q.*3n.  are  shown. 

ot  and  3^^  depend  on  the  TTCIs  for  a  given  a^.  However,  in 
some  cases  fractographic  results  for  selected  fatigue  cracks  had 
to  be  extrapolated  to  obtain  TTCI  values  for  aQ  =  0.03”.  Extra¬ 
polations  were  made  using  the  generalized  crack  size-time  rela¬ 
tionship  based  on  Eq.  2  and  b?^l.  Since  the  corresponding  "Q” 
and  ”b”  values  for  a  given  fatigue  crack  depends  on  the  fracto¬ 
graphic  crack  size  range  used,  this  explains  the  small  difference 
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TABLE  A3  Model  Parameters  for  Fooled  FHQ  Data  Sets  (Straight  Bore  Holes;  0.005-0.100"  Crack  Sire  Range.) 


CASE 

DATA 

SET 

NO. 

SPECIMENS 

a 

(ksi) 

Z 

LT 

a 

b 

Q*xl0^ 

AVE. 

STD. 

DEV. 

C.O.V.(Z) 

I 

XWPF 

33 

34 

15 

4.245 

0.8706 

11152 

2.045 

2.281 

2.566 

0.588 

22 

,9 

HYWPF 

10 

40.8 

15 

1 

1 

4240 

5.129 

2.175 

1 

1 

LYWPF 

10 

30.6 

15 

1 

26789 

1.211 

3.243 

1 

♦ 

* 

II 

XWPB 

31 

34 

15 

4.299 

0.9422 

16930 

1.952 

3.304 

3.499 

0,419 

12 

.0 

HYWPB 

10 

40.8 

15 

7761 

4.138 

3.212 

1 

1 

1 

LYWPB 

10 

30.6 

15 

I 

* 

34174 

1.165 

3.981 

t 

III 

WPF 

37 

34 

0 

4.309 

0.8707 

16952 

1.503 

2.549 

2.566 

0. 

484 

18 

.8 

XWPF 

33 

34 

15 

11164 

2.046 

2.284 

HWYPF 

10 

40.8 

15 

4247 

5.132 

2.179 

LYWPF 

10 

30.6 

15 

26859 

1.211 

3.253 

IV 

WPB 

32 

34 

0 

4.430 

0.9322 

36942 

1.149 

4.247 

3.583 

0. 

549 

15 

.3 

XWPB 

31 

34 

15 

1 

16974 

1,869 

3.174 

HYWPB 

10 

40.8 

15 

i 

7793 

3,968 

3.092 

LYWPB 

10 

30.6 

15 

1 

f 

34265 

1.114 

3.819 

V 

XWPF 

33 

34 

15 

4.: 

269 

0.8902 

11157 

2.228 

2.486 

2.798 

0.453 

16.2 

HYWPF 

10 

40.8 

4243 

5.611 

2.380 

LYWPF 

10 

30.'6 

26816 

1.319 

3.538 

XWPB 

31 

'34 

16920 

1.562 

2,643 

HYWPB 

10 

40.8 

7754 

3.325 

2.578 

LYWPB 

10 

30.6 

15 

34153 

0.926 

3.162 

VI 

WPF 

37 

34 

0 

4. 

364 

0.8826 

16974 

1.588 

2.695 

2, 

800 

0.449 

16.1 

XWPF 

33 

34 

15 

11174 

2.156 

2.409 

HYWPF 

10 

40.8 

15 

4253 

5.419 

2,305 

LYWPF 

10 

30.6 

15 

26918 

1.276 

3.435 

WPB 

32 

34 

0 

36889 

0.928 

3.425 

XWPB 

31 

34 

15 

16952 

1.512 

2.563, 

LYWPB 

10 

30.6 

15 

34219 

0.895 

3.064 

HYWPB 

10 

40.8 

15 

7777 

3.221 

2.505 
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TABLE  A4  Model  Parameters  for  Pooled  FHQ  Data  Sets  (Straight  Bore  Holes;  0.005-0.050"  Crack  Size  Range) 
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in  the  ot  and  values  for  comparable  cases  in  Tables  A3  and 
A4. 

should  be  a  constant  for  the  EIFS  distribution  to 

be  "generic”.  Since  Q^3£.  depends  on  "ot"  and  "b",  the  Q*go. 

values  are  comparable  only  for  the  data  sets  pooled.  For 

example,  in  Table  A3  the  Q|3£^  values  for  Case  I  are  not 

directly  comparable  with  any  of  the  other  cases  in  the  table. 

Individual  QJ3£^  values  from  Table  A3,  as  well  as  average, 

standard  deviation  and  coefficient  of  variation  values  for  0*Bo  . 

i 

are  tabulated  for  comparison  in  Tables  A5  and  A6 .  The  0*6 o 

1  ^i 

values  don't  vary  significantly  for  different  stress  levels, 

load  spectra  and  %  load  transfer.  This  suggests  that  0*6o. 

'1 

values  are  independent  of  load  spectra,  stress  level  and  % 
load  transfer  through  the  fastener. 

A. 4  STUDY  IFQ  MODEL  RELATIONSHIPS 

The  IFQ  model  parameters  from  Table  A3  for  pooled  data 
sets  were  used  to  evaluate  the  following  assumed  relationships. 

Q  =  K6^  (ai) 

6  =  ka  (A2) 

Q  =  (A3) 

In  Eqs .  Al  through  A3,  A,  K,  k,  J,  ^ and  y  are  empirical  con¬ 
stants  and<r=gross  stress.  Results  are  summarized  in  Table  A7 
and  plots  are  presented  in  Figs.  A2  through  A6 . 
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TABLE  A5  Comparison  of  Q*3n .  Values  for  Different  Spectra  and  Stress  Levels 

i 

(15%  LT,  Straight  Bore  Holes) 


% 

0 

/I  J  \ 

h 

LT 

(ksi) 

F-16  400  Hr. 

B-1 

15 

30.6 

3.538 

3.162 

4.269 

0.8902 

1 

34 

2.486 

2.643 

1 

1 

1 

40.8 

2.380 

2.578 

i 

i 

j 

dve.  Q*3£,^ 

2.801 

2.794 

SCQ*6  ) 

"a 

0.640 

0.320 

%  c.o.v. 

22.8 

11.4 

TABLE  A6 


Comparison  of  Values 


for  Different  %  LT , 


Stress  Levels 


and  Spectra  (Straight  Bore  Holes) 


% 

LT 

a 

(ksi) 

a 

b 

F-16  400  Hr. 

B-1 

0 

34 

2.695 

3.425 

4.3 

164 

0.8826 

15 

1 

30.6 

34 

40.8 

3.435 

2.409 

2.305 

3.064 

2.563 

2.505 

Av 

s( 

% 

e.  Q*3£^ 

C.O.V. 

2.711 

0.510 

18.8 

2.889 

0.437 

15.1 
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Table  A7  Siaranary  of  Model  Parameters  and  Empirical  Constants  for 
Pooled  Data  Sets 
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qTxIO 


’3  10  100 

/^(lOOO  Flight  Hours) 


Fig.  A2  LnQ^  Versus  Ln/3  for  Pooled  Data  Sets  (Case  V;  Table  A3) 

J-  A  • 


134 


<r(Ksi) 

Fig.  A3  Versus  LnO"  for  Pooled  Data  Sets  (Case  I;  Table  A3) 
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cr(Ksi) 

Fig,  A4  Ln/«j,  Versus  Ln  CT  for  Pooled  Data  Sets  (Case  V;  Table  A3) 
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cfxlO 


(TCksi) 


Fig.  A5 


Versus  LnCT  for  Pooled  Data  Sets 


(Case  I;  Table  A3) 
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Fig.  A6  I^Q*  Versus  Lucr  for  Pooled  Data  Sets  (Case  V;  Table  A3) 
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Q*3£i  values  should  be  constant  for  different  data  sets  for 
EIFS  values  to  be  "generic'^  This  means  the  constant  "A"  in  Eq. 

A1  should  be  -1.0  for  a  generic  EIFS. .  A  typical  plot  of  In  0^ 
versus  In  is  shown  in  Fig.  A2  for  Case  V.  "A”  values, 

summarized  in  Table  A7  for  various  cases,  vary  from  -0.7861  to 
-0.8571. 

QJ  versus  8^^  values  from  Table  A3  are  plotted  in  Fig.  A2 
for  selected  fighter  and  bomber  data  sets.  This  plot  suggests 
that:  (1)  Q  and  3  values  for  different  data  sets  are  directly 
comparable  when  they  are  determined  using  consistent  data  pool- 
ing  procedures,  C2)  Eq.  A1  is  a  promising  relationship  and  (3) 
is  independent  of  load  spectra  and  stress  level. 

Plots  for  ln3£^  versus  In  o (ksi)  are  shown  in  Figs.  A3  and 
A4.  The  plot  in  Fig.  A3  is  for  three  fighter  data  sets  and  the 
plot  in  Fig.  A4  is  for  pooled  fighter  and  bomber  data  sets.  Two 
parallel  straight  lines  are  shown  in  Fig.  A4  for  the  fighter  and 
bomber  data  sets.  These  lives  were  determined  using  a  least 
square  fit  with  Eq.  A2  and  a  common  slope  (i.e.,  a  constant  J  value) 
imposed.  Eq.  A2  correlates  very  well  with  the  plotted  results. 

Plots  for  In  Q*  versus  In  a  (ksi)  are  shown  in  Figs.  A5  and 
A6,  Eq.  A3,  based  on  a  least  squares  fit,  fits  the  plotted  results 
very  well.  Constants  5  and  ^  are  summarized  for  various  data  set 
combinations  in  Table  A7. 
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A  service  crack  growth  master  curve  is  needed  to  make  p(i,T) 
predictions  using  the  EIFS  ctamulative  distribution,  Eq. 

A3  could  be  particularly  useful  for  durability  analyses  if  the 
Q[a(t)]^  model  is  acceptable  for  defining  the  service  crack  growth 
master  curve  for  different  design  stress  levels.  For  example,  Eq. 
A3  might  be  used  to  help  define  y^^^  (t)  (i.e.,  the  EIFS  value 
corresponding  to  a  crack  size  at  service  time  t)  for  different 
stress  levels  than  the  f ractographic data  base  used  to  define  the 
IFQ  distribution.  Such  an  approach  has  been  used  for  the  final 
demonstration  of  the  durability  methodology  for  the  F-16  lower 
wing  skin  [36] . 
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APPENDIX  B 


EVALUATION  OF  IFQ  MODEL  PARAMETERS  FOR 
COUNTERSUNK  FASTENER  HOLES 

B.l  INTRODUCTION 

A  similar  study  to  the  one  documented  in  Appendix  A  is 
presented  in  this  appendix  for  countersunk  fastener  holes. 

In  this  study  the  IFQ  model  parameters  were  determined  for 
different  combinations  of  pooled  fractographic  data  sets. 

The  model  parameters  were  determined  using  the  preliminary 
calibration  procedures  described  in  Section  4.3,  the  fracto¬ 
graphic  data  from  Ref.  6  and  the  computer  program  from  Ref.  7. 
Also,  the  empirical  relationships  of  Eqs .  A1-A3  (Appendix  A) 
were  also  evaluated. 

The  fifteen  fractographic  data  sets  used  for  this  study 
are  described  in  Table  B1  and  details  are  given  in  Ref.6,:^ 

Details  for  the  15%,  30%  and  407o  load  transfer  specimens  are 
shown  in  Figs.  19-21.  All  fastener  holes  were  drilled  with  a 
modified  Winslow  spacematic  drill  and  the. holes  were  not  deburred. 

This  study  reflected  the  following: 

o  e=0  and  (ref.  Eq.  8) 

o  b?^!  form  of  Eq.  2 

o  a  fractographic  crack  size  range  of  0.005"=a(t)=0.10" 
was  used. 

It  should  be  noted  that  the  procedures  used  to  determine  the 
IFQ  Model  parameters  herein  deviated  slightly  from  those  recom- 
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TABLE  B1  FRACTOGRAPHIC  DATA  SETS  USED 


DATA 

DIA* 

7.  LOAD 

(T 

TEST 

SET 

■  (IN) 

TRANSFER 

(ksi) 

SPECTRUM 

MATERIAL 

AFXLR3 

3/16 

15 

32 

F-16  400 

7475- 

•T7351 

AFXMR3 

34 

HR 

BLOCK 

Aluminum 

AFXHR3 

I 

38 

AFXLR4 

1/4 

32 

AFXMR4 

34 

AFXHR4 

1 

1 

38 

ABXLR4 

1/4 

36 

B 

-1 

ABXMR4 

38 

Bomber 

ABXHR4 

40.8 

ABYLR4 

3 

0 

34 

ABYMR4 

1 

1 

36 

ABYHR4 

1 

1 

38 

ABZLR4 

40 

34 

ABZMR4 

1 

36 

ABZHR4 

j 

P 

38 

*MS90353  Countersunk  Rivet  (Blind,  Pull-through) 
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mended  in  the  Durability  Design  Handbook  [21] .  This  study  and 
others  (e.g.,  11-19,  31-36)  provided  the  ''experience''  needed  to 
develop  the  recommended  calibration  procedures  and  guidelines 
for  the  program  [21] . 

IFQ  model  parameters  are  summarized  herein.  Further  details 
about  the  study  and  conclusions  are  given  in  subsection  4.8.1. 

B.2  IFQ  MODEL  PARAMETERS  FOR  INDIVIDUAL  DATA  SETS 

IFQ  Model  parameters  are  stimmarized  in  Table  B2  for  pooled 
fractographic  data  sets.  The  parameters  were  determined  using 
the  fractography  for  the  largest  fatigue  crack  per  specimen 
(i.e.,  no  6  scaling  used). 

In  Table  B2,  the  IFQ  Model  parameter  notations  are  consis¬ 
tent  with  those  of  Ref.  21  .  Refer  to  Section  A2  of  Appendix  A 
for  a  description  of  the  subscripted/superscripted  notations. 

The  Of  and  b  values  in  Table  B2  were  determined  using  the 
pooled  fractographic  data  sets  indicated.  Once  the  pooled  "a” 
value  had  been  determined,  the  corresponding  value  for 

data  set  pooled  was  determined  using  the  computer  program 
of  Ref.  7, 

values  are  tabulated  in  Table  B2  for  various  fracto- 
S^s-phic  data  pooling  combinati'ons ,  In  Table  B2,  the  average, 
standard  deviation  and  coefficient  of  variation  for  are 

also  shown. 
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TABLE  B2  Model  Parameters  for  Pooled  Data  Sets  (Csk  Rivets:  0.005"-0. 100"  Crack  Size  Range) 
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TABLE  B2  Model  Parameters  for  Pooled  Data  Sets  (Csk  Rivets:  0.005-0.100"  Crack  Size  Range) 
(Continued) 
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One  of  the  conditions  for  a  ’’generic'"  EIFS  distribution  is 
that  should  be  a  constant  for  a  given  material,  drilling 

procedure,  fastener  type  and  fit  £21].  Since  values  depend 

on  the  applicable  and  ”b”,  the  Q*3££  values  are  comparable 
only  for  the  data  sets  pooled.  For  example,  in  Table  B2  the 
average  values  for  Cases  V  and  VI  are  not  compatible  be¬ 

cause  the  ”a"  and  ”b”  values  are  different  for  each  case.  It  has 
also  been  shown  that  Q  and  b  are  cross  correlated  constants  [33] , 
i.e.,  Q  depends  on  b  and  vice  versa, 

Q*/^  values  and  statistics  for  selected  fractographic  data 
sets  for  case  X  of  Table  B2  are  compared  in  Tables  B3-B5.  For 
example,  individual  values  are  shown  for  different  situations 

and  the  average,  standard  deviation  and  coefficient  of  variation 
-k  "k 

for  are  also  shown.  values  are  compared  for  different 

stress  levels  and  load  spectra,  for  different  hole  diameters  (D  = 
1/4”  and  3/16”)  and  for  different  %  load  transfer  and  stress  levels 
in  Tables  B3,  B4  and  B5,  respectively. 


B.3  STUDY  IFQ  MODEL  RELATIONSHIPS 

Empirical  relationships,  given  by  Eqs.  A1-A3,  were  evalu¬ 
ated  using  results  for  Case  X  from  Table  B2.  Plots  for  InQ^ 
versus  in  In  versus  In  a  and  In  versus  In  a  are 

presented  in  Figs.  Bl,  B2  and  B3,  respectively.  The  empirical 
constants  for  each  figure  were  determined  using  a  least  squares 
fit. 
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ic 

TABLE  B3  Comparison  of  Values  for  Different  Stress  Levels  and  Spectra 

(15%  LT.  1/4"  Dia.  Csk  Hole) 
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TABLE  B4 


•k 

Conparlson  of  Q^B^^Values  for  Different  Hole  Diaiaeters 

(F-16  400-Hour  Spectrum;  Csk  Hole) 
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TABLE  B5  Comparison  of  Values  for  Different  %  LT  and  Stress  Levels 

(B-1  Bomber  Spectrum;  Csk  Hole) 


% 

LT 

<r 

(ksi) 

opi, 

Av..  (f Bj. 

c.o 

(Z 

.V. 

) 

1 

5 

36 

38 

40.8 

0.725 

0.559 

0.574 

1 

0.( 

i 

j19 

1 

0.0 

1 

92 

1 

14 

' 

t 

.8 

3 

0 

34 

36 

38 

0.721 

0.521 

0.668 

0.( 

537 

0.1 

04 

16 

.3 

4 

0 

1 

34 

36 

38 

0.574 

0.566 

0.643 

0. 

\ 

594 

L 

0.0 

1 

42 

L 

7 

.1 

L 

Ave, 

1 

1.) 

1 

%  c.o.v. 

0.617 

0.075 

12.1 

a-  2.477 

b  -  0.6097 

149 


2.477 
b  =  0.609 


aQ  “  0.03" 
0.005"=a(t)=0.100" 


Fig.  Bl  Versus  Pooled  Data  Sets  (Case  X;  Table  B2) 
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xlO 


oL  =  2.477 
b  »  0.609 
ao  =  0.03” 
0.005”»a(t)=0.100” 


0^  (F-16)  “  1.579x10 

-13  5.638 

Q*  (B-1)  =  0.656x10 


Fig.  B2  LnQ*  Versus  LnCT  for  Pooled  Data  Sets  (Case  X;  Table  B2) 
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Ln/^  Versus  LnCTfor  Pooled  Data  Sets  (Case  X;  Table  B2) 
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APPENDIX  C 


RELATED  DURABILITY  STUDIES 

C.l  SCALE-UP  EFFECTS  AND  HOLE  INTERACTION 
C.1.1  Introduction 

In  general,  experimental  results  on  the  crack  growth  observed  in  coupon 
specimens,  full-scale  structural  components  and  prototype  structures  manifest 
scale-up  and  interaction  effects.  When  these  test  specimens  have  fastener 
holes,  it  is  reasonable  to  assume  that,  under  the  usual  loading  conditions, 
the  crack  growth  always  concentrates  at  the  points  of  highest  stress  concen¬ 
tration  around  the  fastener  holes. 

The  following  list  indicates  the  major  sources  of  the  scale-up  effects. 

A.  Increase  in  the  number  of  fastener  holes 

B.  Change  in  stress  field 

C.  Increased  variability  in  workmanship 

D.  Increased  variability  in  material  properties 

Preliminary  analytical  models  for  the  scale-up  effects  resulting  from 
these  sources  are  proposed  in  the  following.  It  appears,  however,  that  the 
interaction  effects  can  be  disregarded  since  the  crack  size  to  be  dealt  with 
in  the  present  durability  study  is  at  most  of  the  order  of  50  mils  while  the 
distances  between  the  two  closest  adjacent  fastener  holes  are  obviously  much 
larger  than  50  mils. 

C.1.2  Scale-Up  Effect  Hue  to  Increase  in  Number  of  Fastener  Holes 

Consider  a  set  of  nominally  identical  wing  skin  panels,  each  with  N 
fastener  holes  subjected  to  the  same  state  of  stress.  Let  p(x,t)  denote 
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the  probability  that  the  crack  size  around  a  fastener  hole  will  be  larger 
than  X  at  time  t  and  N(x,t)  the  number  of  holes  with  such  a  crack  size  at 
time  t.  Then,  the  expected  value  N^Cx.t)  of  N(x,t)  can  be  shown  to  be 

ir(x,t)  =  N  p(x,t)  (Cl) 

and  the  variance  a^(x,t)  of  N(x,t) 

a^Cx.t)  =  N  p(x,t)  [1  -  p(x.t)]  (C2) 

if  the  probability  p(x,t)  associated  with  one  fastener  hole  is  independent 
of  the  probabilities  associated  with  other  fastener  holes,  while 

o^(x,t)  =  N2p(x,t)  [1  -  p(x,t)]  (C3) 

if  the  probability  is  completely  dependent.  It  is  recognized,  however,  that 
reality  lies  between  these  two  extreme  cases  and  so  does  the  true  variance. 
The  following  observations  can  be  made  from  Eqs.Cl  -  C3:  (a)  While  the 

expected  value  of  the  number  N(x,t)  of  fastener  holes  with  cracks  larger 
than  X  at  time  t  increases  proportionally  to  the  number  N  of  the  holes  in 
the  panel,  the  variance  a^(x,t)  of  N{x,t)  is  proportional  to  N  or  contin¬ 
gent  upon  whether  the  crack  growth  probability  p(x,t)  is  independent  or  de¬ 
pendent,  and  (b)  this  makes  it  possible  to  determine  whether  the  probability 
is  closer  to  the  case  of  ideal  independence  or  to  that  of  complete  dependence 
by  estimating  the  variance  from  experiments  performed  on  a  set  of  nominally 
identical  panels. 

C.1.3  Scale-Up  Effect  Due  to  Change  in  Stress  Field 
The  probability  p(x,t)  that  each  hole  has  a  crack  larger  than  x  at 
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time  t  depends  on  the  state  of  stress  around  the  hole.  Therefore,  if  the 
i-th  fastener  hole  is  subjected  to  the  stress  state  unique  to  that  hole,  the 
probability  p(x,t)  for  this  hole  should  reflect  such  a  state  of  stress  and 
can  be  denoted  by  p^. (x,t,s)  in  which  s  symbolically  represents  the  state  of 
stress  around  the  hole.  Under  the  assumption  that  the  probability  p^(x,t,s) 
is  independent  of  p.(x,t,s)  (j7i),  it  can  be  shown  that 

J 

N 

N(x,t)  =  I  p.(x,t,s)  (C4) 

i*! 

N 

^^(x.t)  =  I  p.(x,t,s)  [1  -  p,.(x,t,s)]  (Cb) 

i  =  l  ^  ^ 

C.1.4  Scale-Up  Effect  Due  to  Variability  in  Workmanship 

The  spatial  variability  in  workmanship  is  expected  to  produce  variabil¬ 
ity  in  the  stress  state  around  each  fastener  hole.  Therefore,  Eqs.  C4  and 
C5  can  also  be  used  in  this  case  as  the  expected  value  and  variance  of  N(x,t), 
respectively. 

♦  C.1,5  Scale-Up  Effect  Due  to  Variability  in  Material  Properties 

The  existence  of  nonhomogeneity  in  the  material  properties  requires, 
in  principle,  that  the  parameters  of  the  crack  growth  equation  be  adjusted 
accordingly  and  thus  that  different  master  curves  be  used  for  individual 
fastener  holes.  Therefore,  the  probability  p(x,t)  that  the  crack  size  will 
be  larger  than  x  around  a  fastener  hole  at  time  t  varies  from  fastener  hole 
to  fastener  hole,  and  hence  it  will  be  denoted  by  p^(x,t).  Then,  Eqs.  C4 
and  C5  are  again  valid  for  the  expected  value  and  variance  of  N(x,t). 
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C.1.6  Concluding  Remarks 


The  analytical  expressions  for  the  expected  value  and  variance  of  the 
number  of  fastener  holes  with  a  crack  larger  than  x  at  time  t  are  derived 
in  terms  of  the  total  number  of  fastener  holes  in  the  specimen  and  the  prob¬ 
ability  that  such  a  crack  will  be  found  around  a  fastener  hole  at  time  t. 

This  probability  can  be  estimated  from  a  crack  growth  analysis  based  on  the 
fracture  mechanics  theory  and  experimental  evidence.  The  significance  of 
the  dependence  and  independence  of  the  same  probabilities  associated  with 
individual  fastener  holes  has  been  pointed  out.  The  way  in  which  these 
probabilities  depend  on  the  change  in  the  stress  field  and  the  variabilities 
due  to  workmanship  and  material  properties  has  also  been  indicated.  Develop¬ 
ment  of  a  data  base  for  more  detailed  theoretical  and  numerical  analysis  is 
strongly  recommended. 

C2.  DEVELOPMENT  OF  TRANSFER  FUNCTIONS 

C.2.1  Introduction 

Figured  schematically  shows  five  crack  growth  curves  for  a  crack  of 
initial  crack  size  c^  under  five  different  load  spectra  I,  II,  V. 

These  crack  growth  curves  imply  that  the  intensity  of  load  spectrum  I  is 
the  highest,  that  of  II  the  second  highest  and  so  on.  A  question  immedi¬ 
ately  arises  as  to  how  such  a  (spectral)  intensity  can  be  quantitatively 
defined.  One  possible  way  in  which  this  can  be  achieved  is  to  use  a  single 
parameter  indicator  of  the  intensity  such  as  the  root  mean  square  (rms) 
value  of  the  rise  and  fall  of  the  stress  history.  Other  indicators 

involving  two  or  more  parameters  may  be  found  to  be  desirable.  However, 
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since  the  present  study  places  its  emphasis  on  the  ej^loratory  development  of 
transfer  function  methodology,  the  rms  value  of  the  rise  and  fall  of  a  stress 
history  is  used  for  the  definition  of  the  spectral  intensity  primarily  be- 
cause  of  its  simplicity. 

F1gureC2also  schematically  shows  the  relationship  between  the  spectral 
intensity  and  the  flight  hours  needed  for  a  crack  to  grow  from  its  initial 
size  c^  to  5c„  (Curve  A),  lOc^  (Curve  B)  and  20c^  (Curve  C).  For  example. 
Curve  B  indicates  that  Tj,  T^,  ....  Tj  flight  hours  are  required  for  a  crack 
of  initial  size  c„  to  grow  to  that  of  size  lOc^  under  the  spectral  intensity 

associated  with  the  spectra  I.  II,  ...,  V,  respectively. 

Both  Figs, Cl  andC2show  the  transfer  characteristics  among  load  spectra; 
Fig.Clmay  be  used  to  trace  the  process  of  crack  growth  under  a  combined  load 
spectrum.  For  example,  the  initial  crack  of  size  c^.  subjected  to  load  spec¬ 
trum  III,  grows  to  10c„  (Point  P^  in  Fig  cl)  after  T3  flight  hours.  This 
crack  of  size  lOc^  then  grows  under  spectrum  II  to  20c„  after  flight 
hours  (P-  ♦  P3).  If  this  crack  of  size  BOc^  is  further  subjected  to  load 
spectrum  IV,  it  will  grow  to  the  critical  size  c^^  after  T‘  flight  hours 

(P  *  p,).  On  the  other  hand.  Fig  C2  can  be  used  to  establish  a  crack 
'4  5  ^ 

growth  curve  corresponding  to  a  load  spectrum,  say,  1^,  for  which  no  crack 
growth  curve  is  available.  This  can  be  done  by  assuming  that  under  spec¬ 
trum  Ip,  the  initial  crack  of  size  c^  will  grow  to  5Cp,  lOCp  and  20Cp  after, 
respectively.  1(^2  and  7^3  flight  hours  and  interpolating  curve  Ip 

through  the  points  Q^.  Q2  and  Q3  plotted  in  Fig  Cl  which  corresponds  to  the 
points  Qj,  Q2  and  Q3  in  Fig.  C2  In  the  following,  the  curves  in  Fig.C2  are 
referred  to  as  generalized  S-N  curves  for  their  obvious  similarity  to  the 
well-known  S-N  fatigue  curves. 


157 


Spectral  Intensity 


flight  hours  (T) 


Fig.  Cl  Schematic  Comparisons  of  Crack  Growth  Curves 
for  Different  Load  Spectra 


Fig.  C2  Relationship  Between  Spectral  Intensity  and  Flight  Hours 
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C.2.2  Development  of  Crack  Growth  and  Generalized  S-N  Curves 

A  computer  program  named  "CYCGRO"  was  originally  developed  by  the  Rock¬ 
well  International  Corporation  and  later  modified  by  Modern  Analysis  Inc.  to 
accomodate  a  random  stress  sequence  as  the  input  [Cl,].  The  program  computed 
the  size  of  a  crack,  as  it  grows  from  its  initial  size  c^^,  by  integrating  the 
cycle-by-cycle  (differential)  equation  representing  the  crack  growth  rate  in 
terms  of  the  range  of  the  stress  intensity  factor,  current  crack  size,  R- 
factor,  etc.  The  modified  version  of  the  program  also  includes  a  routine  in 
which  a  crack  growth  curve  is  constructed  by  evaluating  the  rate  of  crack 
growth  per  flight  (rather  than  per  cycle)  at  various  values  of  the  existing 
crack  size  with  the  aid  of  the  unit  block  approach  ['C2.]. 

The  CYCGRO  program,  together  with  the  unit  block  approach,  is  also  util¬ 
ized  in  the  present  investigation  to  construct  crack  growth  curves  and  cor¬ 
responding  generalized  S-N  curves.  Five  of  the  more  realistic  load  spectra 
are  considered  as  listed  below.  In  all  cases,  the  initial  crack  size  is  as¬ 
sumed  to  be  Cq  =  0.01". 


ID 

LOAD  SPECTRA 

RMS 

I 

F-15  (Air-Air;  =  40  ksi) 

14.8  ksi 

II 

GD  400  hr  (oj_jn^  =  34  ksi) 

9.93  ksi 

III 

GD  400  hr  =  28.9  ksi) 

8.44  ksi 

IV 

B-1  (Combat) 

5.80  ksi 

V 

B-1  (Crew  Training) 

5.18  ksi 
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Load  spectrum  I  is  defined  by  the  random  stress  sequence  listed  in  Table  Cl 
which  represents  a  block  of  fifty  flights  (totaling  46  flight  hours)  used 
in  the  unit  block  analysis  approach.  Each  stress  level  of  this  sequence 
is  in  terms  of  the  percentage  of  the  limit  stress  which,  for  numerical 
purposes,  is  taken  as  40  ksi.  Load  spectra  II  and  III  are  also  given  in 
the  form  of  a  random  stress  sequence  in  Table C2 with  each  stress  level  again 
indicating  the  percentage  of  Load  spectra  II  and  III  are  obtained  by 

setting  =  34  ksi  and  =  28.9  ksi,  respectively,  in  this  sequence. 

In  each  case,  the  size  of  one  unit  block  used  for  computing  the  flight-by¬ 
flight  crack  growth  rate  involves  400  flight  hours.  TableC3  shows  a  rep¬ 
resentative  stress  sequence  in  one  complete  flight  under  spectrum  IV  while 
Tab! eC4 lists  the  same  under  spectrum  V.  In  both  cases,  the  size  of  one 
unit  flight  consists  of  fifty  flights.  However,  one  unit  block  is  equiva¬ 
lent  to  450  and  200  flight  hours  for  spectra  IV  and  V,  respectively.  The 
rms  values  associated  with  the  load  spectra  listed  above  are  evaluated  di¬ 
rectly  from  these  stress  histories.  The  crack  growth  curves  for  all  the 
load  spectra  are  computed  using  the  CYC6R0  program  and  are  shown  in  Fig.  C3. 
From  these  crack  growth  curves,  the  time  T  (flight  hours)  needed  for  the 
initial  crack  size  c^  =  0.01"  to  grow  to  the  sizes  c  =  2Cq,  6Cq,  Sc^ 

and  10c  are  read  and  plotted  for  each  load  spectrum  which  is  now  represented 
0 

in  terms  of  the  rms  stress  intensity  value  as  shown  in  Fig.C4:  The  five  dif¬ 
ferent  values  of  T  needed  for  c^^  to  reach  2Cq,  4Cq,  Sc^,  Sc^  and  lOc^  at 
each  load  spectrum  intensity  level.  Smooth  curves  are  then  drawn  by  inter¬ 
polation  through  the  five  sets  of  five  points  in  Fig.CA,  with  each  set  indi¬ 
cating  a  generalized  S-N  curve  or  the  relationship  between  the  spectral  in¬ 
tensity  and  the  flight  hours  T  needed  to  produce  a  specific  value  of  c:  For 
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Table  Cl  F-15  Load  Spectrum  (Air -Air) 


Notes  :  Random  spectrum. 

A  block  of  fifty  (50)  flights  (46  hrs.  in  block). 
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Table  Cl  F-15  Load  Spectrum  (Air-to-Air)  (Continued) 
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Table  Cl  F-15  Load  Spectrum  (Air-to-Air)  (Continued) 
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Table  Cl  F-15  Load  Spectrum  (Air-to-Air)  (Continued) 
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Table  Cl  F-15  Load  Spectrum  (Air-to-Air)  (Continued) 
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Table  Cl 


F-15  Load  Spectrum  (Air-to-Air)  (Continued) 
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Notes  : 


Table  UZ  F-16  400  Hc3ur  Spectrum  v  ri. u^raui; 


Percentage  of  <y, 
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Table  C2  F-16  400-Hour  Spectrum  (FHQ  Program)  (Continued) 
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Table  C2  F-16  400-Hour  Spectrum  (FHQ  Program)  (Continued) 
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Table  C2  -  F-16  400-Hour  Spectrum  (FHQ  Program)  (Continued) 
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14.5 

48.2 

0.2 

15.7 

65.0 

5. 

1 3-  8 

33.D 

m 

15.3 

56.9 

3. 

_ -1.3 

6.9 

1* 

— — SSst? 

1 

6-9 

21.0 

11. 

13.8 

53.9 

7, 

- 9i6  -  ■ 

63*9  - 

-  6. 
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Table  C2  F-16  400-Hour  Spectrum  (FHQ  Program)  (Continued) 


13. A 

41.3 

8. 

15.3 

37.5 

122. 

lO.'tJ' 

6T  •  u 

224  . 

13.4 

52.2 

1. 

11.9 

57.7 

9. 

15.3 

75.8 

5. 

15.3 

51.0 

34. 

9.6 

49.3  . 

25. 

-2.7  — 

60.8^ 

-  3. 

13.4 

19.1 

1. 

15.3 

56.0 

136. 

1 

i 

N 

• 

'Urn  ^ 

1  • 

13.8 

44.4 

3. 

13.4 

34.5 

29. 

13.4  — 

-  29.3 

■  57. 

-7.6 

14.2 

t  « 

2.1 

34.3 

4. 

6T9 

■  61 ;9 — 

14. 

15.7 

51.4 

17. 

13.4 

33.7 

2. 

12-4 

-  27.4 

30. 

14.2 

40.4 

47. 

11.9 

71.6 

4. 

- 13;4 - 

-41;^-  - 

-  52. 

0.0 

15.3 

57. 

15.7 

37.5 

69. 

—  —  «  4L 

lli^  — 

30 . 7 

76. 

11.5 

44.2 

14. 

4.4 

15.9 

1. 

-1472 

— 1374- 

-  1. 

2.7 

11.9 

1. 

11.5 

38.8 

479. 

- 15.7- 

-  50.=^- 

“5. 

11.9 

44.0 

22. 

13.4 

16.9 

1. 

"12. ■4' 

39  72— 

2. 

15.3 

58.1 

6. 

■  ”11.9 

- 51.U  ” 

26. 

15.3 

58*1 

3. 

-2.7 

11.5 

1. 

1  1 

oO.'T" 

1014. 

11.9 

32.0 

109. 

14.5 

44.4 

6. 

- 1573' 

- 53^7 

26. 

3.1 

19.2 

7. 

-7.1 

43.5 

1. 

•T  •  1 

• 
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Table  C3  B-1  Load  Spectrinn  (Combat) 


Notes:  All  fractional  cycles  have  been  replaced  by  unit  cycles. 

This  is  a  complete  flight  (  *  9  hrs.). 


MISSION 

SEGMlNT 

WING 

ANue 

LOAP  I 

WtosS  STRESS LKS1>  < 

LYClCS/ 

mssicw 

typl 

MAxiTli;^  ll 

OKOUND 

IS 

20. OC 

o 

o 

i 

1.00 

POST  ta.i  orr 

IS 

MMEUVCR 

17  ee 

10.09 

0.10 

f4r.e8 

11.07 

1.00 

WING  SWEEP 

IS 

10 

19.67 

13.87 

1.00 

?s 

II. es 

11.26 

i.oe 

SUBSONIC  CLll'B 

^s 

MANt'JVEP 

ie.63 

9.77 

fi^io 

I5W6 

7.96 

1  .00 

1  5  .E.7 

1  1  .BG 

MoOO 

I  I.Z6 

O.BG 

3.00 

SUBSONIC  CRUISE 

Bb 

OUST 

1 S.  TO 

6.^0 

0.10 

m  ,39 

7. DO 

I. 00 

li .  If 

9.00 

1 .00 

maneuvtr 

i6.fo 

7.19 

O.iO 

16.20 

7*H9 

1.00 

If.BD 

1  1  .09 

19.00 

11.09 

8.39 

6.00 

WING  SWEEP 

Bb 

10 

1  1 

M  .^6 

uoo 

67. S 

7.70 

7.70 

1  .00 

SUPERSONIC  CLIMB 

67. S 

MAAdtTVtd 

l¥-.t6 

2.02 

0.10 

12.32 

M  .03 

o 

o 

9 

I0,Z5 

B  17 

$.oo 

8. 17 

6.6^ 

7.00 

SUPERSONIC  CRUISE 

L7.S 

1.92 

0.  lU 

10.96 

3.96 

1.00 

8.97 

6.33 

5.00 

6.33 

9.80 

Z.OO 

SUPERSONIC  DESCLNl 

67.S 

12.69 

1  .57 

0.10 

10.65 

3.97 

1.00 

t.»7 

7.17 

B.OO 

7. 17 

5.51 

7.00 

WING  SWEEP 

67.5 

IG 

7.70 

7.70 

1  .00 

?S 

1 1. 00 

It  09 

1  .00 

REFDCL 

Bb 

MANEUVEK 

<5  C9 

5-59 

0.10 

19.  G9 

6.89 

1 .00 
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Table  C3  B-1  Load  Spectrum  (Combat)  (Continued) 


173 


Table  C3  B-1  Load  Spectrum  (Combat)  (Continued) 
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laDie  OH  D-x  Xioau  Dpecxrum  \.orew  ±raxnxag; 


Notes:  All  fractional  cycles  have  been  replaced  by  unit  cycles. 

This  is  a  complete  flight  (  =  4  hrs.). 


MISSION 

SCGMCNT 

wlf;'!' 

angle 

LOAD  1 

GROGS  SIR£SSCk^-l>  ( 

:tcips/ 

vpsSIOfi 

TYPC  1 

MAXIMUM  i 

MINIMUM 

GPOUND 

15 

80  .00 

-  5.00 

1  .  UU 

POST  take  orr 

15 

maneuver 

17.05 

o 

o 

o.ot 

15.57 

8.H9 

0.  10 

m  .36 

9  69 

1  .00 

WING  SWEEP 

15 

IG 

12.0? 

18.08 

1.00 

?5 

10.02 

10.08 

1  .00 

SUBSONIC  CLIMB 

?5 

MANEUVER 

16.00 

3. 60 

0.01 

15.03 

H  .85 

0.  10 

13.36 

6.76 

1  .00 

18.  IM 

7.90 

8.00 

subsonic  cruise 

?5 

GUST 

15.90 

7.00 

0.01 

IH.SO 

5.81 

0. 10 

18.50 

7.  M 

1  .00 

MANEUVER 

IS. 65 

5.36 

0.01 

IS. 78 

6.83 

0  10 

12.9*1 

6.  10 

1.00 

12. 2S 

9.60 

8.00 

REFUEL 

?5 

MANEUVER 

IS. 68 

3.H6 

0.01 

IS. OB 

H.66 

0.  10 

1 

12. 9S 

6.83 

1  .00 

1 1 .91 

9.00 

7.00 

9.80 

7.98 

11.00 

SUBSONIC  CRUISE 

?5 

MANEUVER 

16.58 

H.6G 

'  0.01 

15. HO 

5.53 

O'.IO 

13. 6H 

7.00 

1  .00 

18.?^» 

9.80 

o 

o 

9.80 

7.97 

1 .00 

SUBSONIC  DESCENT 

85 

MANEUVLH 

IS.3S 

5.07 

O.Ol 

IS. 31 

5.75 

0.10 

12.52 

7.81 

1.00 

11.67 

7.99 

1.00 

WING  SWEEP 

85 

IG 

9.70 

9.78 

1  .00 

67.5 

7.17 

7,17 

1.00 
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Table  C4  B-1  Load  Spectrum  (Crew  Training)  (Continued) 


MISSION 

WING 

LOAD 

GROSS  STRESS iKSl »  i 

CYCLES/ 

5CGMCN7 

angle 

TYPE 

MAX  1  MUM 

MINIMUM  ! 

MISSION 

.CU  LtVCL  NAVIGATION 

67. b 

GU61 

8.11 

0.61 

0.01 

7.73 

0.89 

0.  ID 

6.97 

1.66 

1 .00 

6.81 

8. Ml 

10.00 

6.70 

8.98 

17.00 

n  .m 

8.98 

0.01 

10.70 

3.36 

0.10 

9.9M 

M.I8 

1.00 

9.86 

M.e8 

10.00 

6.61 

5.S6 

61 .00 

iH.eo 

6.96 

O.Ol 

m.3S 

7. HD 

0. 10 

13.60 

6.16 

1.00 

18.91 

6.66 

10.00 

18. 3H 

9. M3 

86.00 

FLY  UP 

67.6 

3.0G 

16.01 

7.03 

1 .00 

a.5G 

15. as 

7.03 

1  .00 

.D;-J  LEVLU  NAVIGATION 

67.6 

MAI'H.UVER 

13.89 

3.03 

o.ot 

13. 

3. MM 

O.ID 

18.63 

3.59 

1 .00 

18.66 

3.76 

10.00 

11.96 

M.OO 

86.00 

10.63 

M.31 

60.00 

9.80 

S.78 

60.00 

6.73 

6.6M 

60.00 

LOW  level  bombing 

67.6 

GUST 

10.98 

8.0G 

0.01 

10.5m 

8. MM 

o 

o 

9.8M 

3.13 

1.00 

9.08 

3,96 

5.00 

6.19 

M.7e 

17.00 

maneuver 

6.10 

S  .BS 

0.01 

7.98 

5.05 

0.10 

7.67 

5.30 

1 .00 

7.M6 

5.51 

a.  00 

WING  5WELP 

67.0 

IG 

6.S9 

6.M9 

1 .00 

B5 

6.60 

8.60 

1 .00 

SUBSONIC  CLIMB 

86 

MANEUVER 

IS. 03 

3. Ml 

p.oi 

13.19 

H.Ml 

0.10 

1 1 .68 

6.00 

1 .00 

10.6a 

6.96 

a. 00 

SUBSONIC  CRUISE 

as 

MAItfUVER 

m.6o 

M.80 

0.01 
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Table  C4  B-1  Load  Spectrum  (Crew  Training)  (Continued) 


hlSSlON 

wing 

LOAD  I 

C»RC’^^  SI KCSw •  ^  -  1  ' 

:  ycles/ 

SECMCNT 

ANGl  L 

TYPL  1 

MAXIMUM  1 

minimum  I 

•IISSION 

SI-IBSONIC  CRUIS-E 

?b 

MANfUVl  fy 

13.5? 

s  .9r, 

0. 10 

1?.C9 

6.85 

1  .00 

10. eo 

8.67 

5.00 

e.67 

7.09 

1 ,00 

SUBi-ONlC  DESCENT 

ib 

r.AMuvCR 

1  1  .S5 

8.15 

0.01 

10.5? 

1  .69 

0.10 

8.97 

M.09 

1  .00 

8.  I? 

^  .79 

1  .00 

WING  SWEEP 

85 

IG 

e.Bo 

8. 60 

1  .00 

15 

9.?5 

9.85 

1  .00 

PRE- LANDING 

15 

GU3T 

IG.97 

1  .53 

0.01 

15.75 

8.75 

0.  lU 

13.89 

S.68 

1  .00 

l?.S? 

6.08 

10.00 

1  1  .88 

7.88 

93.00 

MANEUVER 

17.63 

8.60 

O.Ol 

16. N9 

3.8^ 

0.10 

m.HO 

s.6e 

1.00 

l?.30 

9.85 

80.00 

9.?5 

6.38 

15.00 

GRODNO 

15 

TAXI 

-1.6! 

-0.76 

1  .00 

FLY  AROUND 

15 

maneuver 

m  .RO 

^.B8 

1  .00 

GROUND 

15 

TAXI 

-1.61 

-0.7G 

1  .00 

Fly  around 

15 

MmIJCUVI  r 

m  .RO 

X 

CD 

CD 

1.00 

GROUND 

15 

lAXl 

-1  .61 

-0.76 

1.00 

FLY  AROUND 

15 

MAtJEUVER 

IS  .so 

^  .89 

1  .00 

ground 

15 

TAXI 

-1.61 

-0.7G 

1.00 

Fly  around 

15 

MANEUVER 

IS  .so 

H.88 

1  .00 

7.S5 

9.85 

1 .00 

177 


(inch) 


O 

o 

o 

m 


o 

o 

o 

o 


o 

o 

o 

LTl 


to 

3 

O 


{q:^6U9[  ^DBJD)  D 
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Fig.  c3  Crack  Growth  Curves  for  Five  Load  Spectra 
Based  on  CYCGRO  Program 


o 

o 

o 
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Fig.  C4  Rms  Stress  Value  -  Flight  Hour  Relationship 


example,  Curve  A  for  c  =  2Cq.  These  curves  appear  to  be  reasonable  in  rep¬ 
resenting  the  T-v^aP^  relationships  for  the  various  crack  sizes,  except  in 
the  lower  range  of  the  spectral  intensities.  Indeed,  when  the  rms  stress 
values  are  small,  the  crack  grows  slowly  and  the  relationships  may  become 
highly  sensitive  to  the  definition  of  the  spectral  intensity.  In  fact, 
under  load  spectrum  IV  with  an  rms  stress  value  of  5.80  ksi,  the  crack  grows 
more  slowly  than  under  load  spectrum  V  with  a  smaller  rms  stress  value  of 
5.18  ksi.  This  particular  anomaly  can  be  explained,  however:  one  flight 
for  load  spectrum  IV  is  nine  (9)  hours  in  duration,  while  one  flight  for 
load  spectrum  V  is  four  (4)  hours.  Therefore,  although  the  rms  value  is 
smaller,  load  spectrum  V  has  a  larger  number  of  air-ground  cycles  and,  ap¬ 
parently  because  of  this,  produces  a  slightly  larger  amount  of  damage  for 
each  flight  hour  than  load  spectrum  IV  does.  Such  anomalies  may  be  recti¬ 
fied  by  redefining  (or  revising  the  current  definition  of)  the  spectral  in¬ 
tensity.  For  example,  a  spectral  intensity  definition  that  reflects  the  ef¬ 
fect  of  the  mean  stress  or  of  the  ratio  R  of  the  maximum  to  minimum  stress 
(in  addition  to  the  rise  and  fall)  may  alleviate  this  difficulty.  This  is 
indeed  an  interesting  subject  to  be  investigated.  It  is  pointed  out  again, 
however,  that  the  T-  /S?  relationships  in  Fig.C4  present  no  anomalous  fea¬ 
tures  in  the  high  intensity  range  where  the  fighter  load  spectral  intensi¬ 
ties  are  expected  to  lie,  particularly  when  each  load  spectrum  involves  ap¬ 
proximately  the  same  number  of  flight  hours  for  one  flight. 

At  any  rate,  the  curves  in  Figs. C3  and C4 represent,  in  approximation, 
the  transfer  characteristics  as  explained  earlier:  For  example,  the  crack 
growth  curve  under  a  load  spectrum  with,  say,  ~  12.5  ksi,  can  be  con¬ 
structed  approximately  on  the  basis  of  Fig.C4by  reading  the  values  of  T 
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at  c  =  2Cq,  4Cq,  etc respectively,  as  850  hrs,  1,700  hrs,  etc.  As  also 
indicated  previously,  it  is  important  to  note  that  these  curves  establish 
the  equivalent  flight  hours  for  the  various  load  spectra.  For  example, 
from  Fig.C3,  one  observes  that  1,300  hours  of  flight  under  load  spectrum  II 
is  equivalent  to  2,250  hours  of  flight  under  load  spectrum  III,  since,  in 
both  cases,  the  initial  crack  c^  =  0.01"  will  grow  into  c  =  0.02"  during 
their  respective  flight  hours. 

In  a  further  application  of  this  concept  of  equivalent  flight  hours,  a 
simplified  spectrum  consisting  only  of  a  few  stress  cycles  (including,  how¬ 
ever,  an  air-ground  stress  cycle)  can  be  constructed  for,  say,  experimental 
purposes.  Fig. C5  depicts  such  a  simplified  spectrum  whose  parameters  are 

a  ,  a  .  ,  a'^  and  n  (number  of  stress  cycles).  For  example,  adjusting 
max  mm  mm 

these  parameter  values,  one  can  construct  five  types  of  load  spectra  iden¬ 
tified  by  ID  Numbers  1  through  5  in  Table C5  The  designation  of  these  spec¬ 
tra  as  A-A  (severe),  A-A,  A-G,  Composite  and  I-N,  is  rather  arbitrary,  how¬ 
ever.  Under  the  assumption  that  one  flight  involves  1.3  flight  hours  for 
A-A,  0.92  hours  for  A-G,  1.22  hours  for  Composite  and  2.07  hours  for  I-N, 
the  crack  growth  curves  are  constructed  and  plotted  in  Fig.C3  (dashed  curves). 
Comparing  curves  1  %  5  with  curves  I  ^  V  in  Fig.C3,  one  observes  that  spectra 
1  .(simplified  A-A,  severe)  and  3  (simplified  A-G)  are  almost  equivalent  to 
spectra  I  (F-15,  A-A)  and  III  (GD  400  hours  with  =  34  ksi),  respec¬ 
tively,  at  least  up  to  c  =  0.1",  and  spectrum  5  (simplified  I-N)  somewhat 
underestimates  the  crack  growth  associated  with  spectrum  III  (GD  400  hours 
with  =  28.9  ksi).  The  other  simplified  spectra  6  15  are  also  listed 

in  Table CSwith  the  corresponding  crack  growth  curves  plotted  in  Fig.C6. 
Comparison  of  Fig.C6with  Fig. C3  shows,  among  other  things,  that  simplified 
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1 - - 


Fig.  C5  A  Simplified  Load  Spectrxim 


y 
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Table  C5  Simplified  Amplitude  FI ight-by-Fl ight  Load  Spectra  =  40  ksi) 


0 


ID 

# 

Mission 

"■min- 

"■min- 

""maK- 

(ao^)*^,  ksi 

n,  cycles 

1 

A  -  A 
(severe) 

-2 

4 

20 

16.28 

24 

2 

A  -  A 

-2 

4 

16 

12.30 

24 

3 

A  -  G 

-4 

4 

16 

12.52 

19 

4 

Composite 

-3 

4 

16 

12.40 

21 

5 

I  -  N 

-.8 

1.6 

14 

12.64 

10 

6 

Composite 

-3 

6 

16 

10.58 

21 

7 

-3 

8 

16 

8.80 

8 

-3 

12 

16 

5.63 

9 

-3 

6 

14 

8.62 

10 

-3 

6 

18 

12.55 

11 

-3 

6 

20 

14.53 

12 

-3 

6 

22 

16.52 

13 

0 

4 

16 

12.21 

14 

-.4 

4 

16 

12.23 

15 

-1.6 

4 

16 

12.31 

Note:  A-A; 

A-G; 

Composite; 

I-N; 


1  fl ight  =  1.3  hrs 
1  flight  =0.92  hrs 
1  flight  =1.22  hrs 
1  flight  =  2.07  hrs 
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crack  length 


Fig.  C6  Crack  Growth  Under  Simplified  Spectra  (Hissions  6-15) 
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spectra  7  and  9  are  equivalent  in  approximation  to  spectrum  V  (B-1,  Crew 
Training),  with  spectrum  7  overestimating  and  spectrum  9  underestimating 
the  crack  growth  associated  with  spectrum  V. 

C.2.3  Concluding  Remarks 

The  transfer  characteristics  among  load  spectra  are  described  with  the 
aid  of  crack  growth  and  generalized  S-N  curves.  The  generalized  S-N  curve 
is  defined  as  showing  the  relationship  between  the  load  spectral  intensity 
(S)  and  the  flight  hours  (N)  required  for  a  crack  of  initial  size  c^  to  grow 
to  a  specific  size  c(>  c^).  For  simplicity,  the  rms  value  of  the  rise  and 
fall  of  the  stress  history  representing  each  load  spectrum  is  used  in  this 
investigation  as  the  indicator  of  the  spectral  intensity.  For  the  five 
more  realistic  load  spectra  I  -v.  V,  crack  growth  and  generalized  S-N  curves 
are  constructed.  Similarly,  the  crack  growth  curves  for  the  simplified  load 
spectra  1  'v  15  are  constructed  for  the  purpose  of  further  comparison. 

C.3  PARAlffiTRIC  SENSITIVITY  ANALYSIS 

In  the  preceding  sections,  a  number  of  analytical  models  are  introduced 
primarily  with  respect  to  the  crack  propagation  law,  the  distribution  func¬ 
tion  of  TTCI  and  the  corresponding  distribution  function  of  EIFS  for  the  pur¬ 
pose  of  evaluating  the  probability  p(x,t)  that,  say,  a  fastener  hole  will 
have  a  dominant  crack  of  a  size  larger  than  x  at  time  t.  Obviously,  the 
value  of  p(x,t)  depends  on  these  analytical  models.  However,  these  models 
are  constructed  on  the  basis  of  the  current  (and  therefore  not  necessarily 
perfect)  state  of  engineering  knowledge  on  the  subject  matters  involved, 
with  due  consideration  not  only  for  analytical  and  numerical  tractability 
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but  also  for  data  availability.  Therefore,  it  is  highly  desirable  to  inves¬ 
tigate  the  sensitivity  of  the  probability  to  these  models.  Such  sensitivity 
studies  can  be  performed  in  terms  of  numerical  examples  for  those  items  for 
which  competing  or  alternative  models  are  proposed.  For  example,  a  signif¬ 
icant  difference  in  the  probability  may  result  depending  on  whether  one 
chooses  a  Wei  bull  compatible  distribution  or  a  log-normal  compatible  distri¬ 
bution  for  the  EIFS.  Such  results  are  highly  useful  in  numerically  assess¬ 
ing  the  significance  of  the  various  models  described  in  the  preceding  sec¬ 
tions. 

After  the  decision  is  made  to  choose  a  particular  model  among  the  al¬ 
ternatives  proposed,  it  is  not  unusual  to  face  uncertainties  of  a  statisti¬ 
cal  or  other  nature  under  which  the  values  are  to  be  assigned  to  the  param¬ 
eters  of  the  model.  Then,  sensitivity  studies  are  again  desirable  with  res¬ 
pect  to  these  parameters  for  the  purpose  of  identifying  the  more  important 
parameters  in  terms  of  their  contributions  to  the  probability.  Such  a  param¬ 
eter  sensitivity  analysis,  which  is  the  main  theme  of  this  sub-section,  is 
more  straightforward  than  a  model  sensitivity  analysis,  at  least  in  principle. 

Writing  these  parameters  as  X^,  X2,  .  X^  ,  the  probability  P  =  p(x,t) 

can  be  formally  written  as  a  function  of  X^,  X2,  .  X^  once  the  analytical 

models  to  be  used  are  decided  upon; 

P  =  P(Xj,  X2,  . X^)  (C6) 

In  actuality,  the  dependence  of  P  on  the  parameters  is  obviously  highly  com¬ 
plex  as  seen  from  the  analysis  developed  in  the  preceding  sections. 

Treating  the  uncertainties  associated  with  these  parameters  as  if  they 
were  all  from  statistical  origins  (or  treating  X^,  X2,  .  X^  as  random 
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variables),  the  sensitivity  of  the  probability  p(x,t)  to  each  of  these  param> 
eters  is  evaluated  with  the  first  order  and  second  moment  approach.  Within 
the  framework  of  this  approach,  the  most  reliable  value  of  the  parameter 
is  treated  as  its  expected  value  and  written  as  X ^ ,  whereas  its  uncertainty 
is  expressed  in  terms  of  the  coefficient  of  variation  with 

being  the  standard  deviation  of  X^.  The  expected  value  "P  and  coefficient  of 
variation  Vp  of  P  are  then  obtained  from 

P  =  P(Xj,  X2 . X^)  (C7) 


where  ot^  is  referred  to  as  the  sensitivity  index  and  is  given  by 


(C9) 


where  the  partial  derivatives  are  evaluated  at  X^  =  X^  {i=l,2,. .. ,n).  If  P 
and  Vp  are  interpreted  as  representing  the  most  reliable  value  of  P  and  its 
uncertainty,  respectively,  the  sensitivity  index  indicates  the  contribu¬ 
tion  of  the  uncertainty  associated  with  the  parameter  X^  to  the  uncertainty 
of  P  through  Eq.  C8. 

As  mentioned  earlier,  the  dependence  of  P  on  X^,  X2,  X^  is  highly 
complex  and  the  analytical  evaluations  of  the  partial  derivatives  in  Eq.  C9 
may  not  be  possible.  Therefore,  these  derivatives  are  numerically  evaluated 
as 
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Similarly,  a  parameter  sensitivity  analysis  can  be  performed  with  res¬ 
pect  to  the  mean  N(x,t)  and  variance  o^(x,t)  of  the  number  N  of  the  fastener 
holes.  In  fact,  it  may  be  more  sensible  to  perform  a  parameter  sensitivity 
analysis  for  these  quantities  since  they  are  more  directly  related  to  the 
durability  performance. 

With  respect  to  a  sensitivity  analysis,  the  present  investigation  has 
been  limited  in  its  scope  to  a  general  discussion.  In  view  of  the  fact  that 
the  present  state  of  knowledge  is  always  imperfect  and  involves  uncertain 
factors,  it  is  highly  recommended  that  a  future  durability  study  address  it¬ 
self  exclusively  to  the  problem  of  model  as  well  as  parameter  sensitivity, 
and  that  it  examine  the  range  of  the  resulting  numerical  solutions  in  order 
to  establish  engineering,  if  not  mathematical,  confidence  in  the  models  and 
their  parameter  values  to  be  used  in  the  durability  analysis  as  developed  in 
the  preceding  sections. 
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